-
Equilibrium and Non-Equilibrium diffusion approximation for the radiative transfer equation
Authors:
Elena Demattè,
Juan J. L. Velázquez
Abstract:
In this paper we study the distribution of the temperature within a body where the heat is transported only by radiation. Specifically, we consider the situation where both emission-absorption and scattering processes take place. We study the initial boundary value problem given by the coupling of the radiative transfer equation with the energy balance equation on a convex domain…
▽ More
In this paper we study the distribution of the temperature within a body where the heat is transported only by radiation. Specifically, we consider the situation where both emission-absorption and scattering processes take place. We study the initial boundary value problem given by the coupling of the radiative transfer equation with the energy balance equation on a convex domain $ Ω\subset \mathbb{R}^3 $ in the diffusion approximation regime, i.e. when the mean free path of the photons tends to zero. Using the method of matched asymptotic expansions we will derive the limit initial boundary value problems for all different possible scaling limit regimes and we will classify them as equilibrium or non-equilibrium diffusion approximation. Moreover, we will observe the formation of boundary and initial layers for which suitable equations are obtained. We will consider both stationary and time dependent problems as well as different situations in which the light is assumed to propagate either instantaneously or with finite speed.
△ Less
Submitted 16 July, 2024;
originally announced July 2024.
-
Asymptotic Analysis of a bi-monomeric nonlinear Becker-D{ö}ring system
Authors:
Marie Doumic,
Klemens Fellner,
Mathieu Mezache,
Juan J L Velázquez
Abstract:
To provide a mechanistic explanation of sustained {then} damped oscillations observed in a depolymerisation experiment, a bi-monomeric variant of the seminal Becker-Döring system has been proposed in \cite{DFMR}. When all reaction rates are constant, the equations are the following: \begin{align*}\frac{dv}{dt} & =-vw+v\sum_{j=2}^{\infty}c_{j}, \qquad \frac{dw}{dt} =vw-w\sum_{j=1}^{\infty}c_{j},…
▽ More
To provide a mechanistic explanation of sustained {then} damped oscillations observed in a depolymerisation experiment, a bi-monomeric variant of the seminal Becker-Döring system has been proposed in \cite{DFMR}. When all reaction rates are constant, the equations are the following: \begin{align*}\frac{dv}{dt} & =-vw+v\sum_{j=2}^{\infty}c_{j}, \qquad \frac{dw}{dt} =vw-w\sum_{j=1}^{\infty}c_{j}, \\ \frac{dc_{j}}{dt} & =J_{j-1}-J_{j}\ \ ,\ \ j\geq1\ \ ,\ \ \ J_{j}=wc_{j}-vc_{j+1}\ \ ,\ \ j\geq1\ \ ,\ J_{0}=0,
\end{align*} where $v$ and $w$ are two distinct unit species, and $c_i$ represents the concentration of clusters containing $i$ units.
We study in detail the mechanisms leading to such oscillations and characterise the different phases of the dynamics, from the initial high-amplitude oscillations to the progressive damping leading to the convergence towards the unique positive stationary solution. We give quantitative approximations for the main quantities of interest: period of the oscillations, size of the damping (corresponding to a loss of energy), number of oscillations characterising each phase. We illustrate these results by numerical simulation, in line with the theoretical results, and provide numerical methods to solve the system.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Multi-parameter Hopf bifurcations of rimming flows
Authors:
Illya M. Karabash,
Christina Lienstromberg,
Juan J. L. Velázquez
Abstract:
In order to investigate the emergence of periodic oscillations of rimming flows, we study analytically the stability of steady states for the model of (Benilov, Kopteva, O'Brien, 2005), which describes the dynamics of a thin fluid film coating the inner wall of a rotating cylinder and includes effects of surface tension, gravity, and hydrostatic pressure. We apply multi-parameter perturbation meth…
▽ More
In order to investigate the emergence of periodic oscillations of rimming flows, we study analytically the stability of steady states for the model of (Benilov, Kopteva, O'Brien, 2005), which describes the dynamics of a thin fluid film coating the inner wall of a rotating cylinder and includes effects of surface tension, gravity, and hydrostatic pressure. We apply multi-parameter perturbation methods to eigenvalues of Fréchet derivatives and prove the transition of a pair of conjugate eigenvalues from the stable to the unstable complex half-plane under appropriate variations of parameters. In order to prove rigorously the corresponding branching of periodic solutions from critical equilibria, we extend the multi-parameter Hopf-bifurcation theory to the case of infinite-dimensional dynamical systems.
△ Less
Submitted 17 June, 2024;
originally announced June 2024.
-
A stochastic version of the Hopfield-Ninio kinetic proofreading model
Authors:
Eugenia Franco,
Juan J. L. Velázquez
Abstract:
In this paper we study a simple stochastic version of the Hopfield-Ninio kinetic proofreading model. The model is characterized by means of two parameters, the unbinding time, which depends on the binding energy between a ligand and a receptor, and the number of times $M \geq 1$ that a ligand attaches to a receptor. We prove that, under suitable assumptions on M, our model has an extreme specifici…
▽ More
In this paper we study a simple stochastic version of the Hopfield-Ninio kinetic proofreading model. The model is characterized by means of two parameters, the unbinding time, which depends on the binding energy between a ligand and a receptor, and the number of times $M \geq 1$ that a ligand attaches to a receptor. We prove that, under suitable assumptions on M, our model has an extreme specificity, i.e. it is capable to discriminate between different ligands, and a high sensitivity, i.e. the response of the system does not change in a significant manner for ranges of ligands varying within several orders of magnitude. Additional quantities like the amount of energy used by the network or the time required to yield a response will be also computed. We also show that our results are robust, i.e., they do not depend on the specific choice of parameters that we make in this paper.
△ Less
Submitted 17 May, 2024;
originally announced May 2024.
-
On an inhomogeneous coagulation model with a differential sedimentation kernel
Authors:
Iulia Cristian,
Barbara Niethammer,
Juan J. L. Velázquez
Abstract:
We study an inhomogeneous coagulation equation that contains a transport term in the spatial variable modeling the sedimentation of clusters. We prove local existence of mass conserving solutions for a class of coagulation kernels for which in the space homogeneous case instantaneous gelation (i.e., instantaneous loss of mass) occurs. Our result holds true in particular for sum-type kernels of hom…
▽ More
We study an inhomogeneous coagulation equation that contains a transport term in the spatial variable modeling the sedimentation of clusters. We prove local existence of mass conserving solutions for a class of coagulation kernels for which in the space homogeneous case instantaneous gelation (i.e., instantaneous loss of mass) occurs. Our result holds true in particular for sum-type kernels of homogeneity greater than one, for which solutions do not exist at all in the spatially homogeneous case. Moreover, our result covers kernels that in addition vanish on the diagonal, which have been used to describe the onset of rain and the behavior of air bubbles in water.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
Characterizing the detailed balance property by means of measurements in chemical networks
Authors:
Eugenia Franco,
Bernhard Kepka,
Juan J. L. Velázquez
Abstract:
In this paper we study how to determine if a linear biochemical network satisfies the detailed balance condition, without knowing the details of all the reactions taking place in the network. To this end, we use the formalism of response functions $R_{ij} (t) $ that measure how the system reacts to the injection of the substance $j$ at time $t=0$, by measuring the concentration of the substance…
▽ More
In this paper we study how to determine if a linear biochemical network satisfies the detailed balance condition, without knowing the details of all the reactions taking place in the network. To this end, we use the formalism of response functions $R_{ij} (t) $ that measure how the system reacts to the injection of the substance $j$ at time $t=0$, by measuring the concentration of the substance $i \neq j$ for $t >0$. In particular, we obtain a condition involving two reciprocal measurements (i.e.~$R_{ij}(t), \, R_{ji}(t)$) that is necessary, but not sufficient for the detailed balance condition to hold in the network. Moreover, we prove that this necessary condition is also sufficient if a topological condition is satisfied, as well as a stability property that guarantees that the chemical rates are not fine-tuned.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
The three dimensional magneto-hydrostatic equations with Grad-Rubin boundary value
Authors:
Diego Alonso-Orán,
Daniel Sánchez-Simón del Pino,
Juan J. L. Velázquez
Abstract:
In this work, we study the well-posedness of the three dimensional magneto-hydrostatic equation under Grad-Rubin boundary value conditions. The proof relies on a fixed point argument to construct solutions to an elliptic-hyperbolic problem in a perturbative regime by means of pseudo-differential operators with symbols with limited regularity in Hölder spaces. As a byproduct, the employed technique…
▽ More
In this work, we study the well-posedness of the three dimensional magneto-hydrostatic equation under Grad-Rubin boundary value conditions. The proof relies on a fixed point argument to construct solutions to an elliptic-hyperbolic problem in a perturbative regime by means of pseudo-differential operators with symbols with limited regularity in Hölder spaces. As a byproduct, the employed technique in this work is more flexible and simplifies the arguments of the proof for the previous two-dimensional setting.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
Interface behavior for the solutions of a mass conserving free boundary problem modeling cell polarization
Authors:
Anna Logioti,
Barbara Niethammer,
Matthias Röger,
Juan J. L. Velázquez
Abstract:
We consider a parabolic non-local free boundary problem that has been derived as a limit of a bulk-surface reaction-diffusion system which models cell polarization. In previous papers, we have established well-posedness of this problem and derived conditions on the initial data that imply continuity of the free boundary as $t\to 0$. In this paper we extend the qualitative study of the free boundar…
▽ More
We consider a parabolic non-local free boundary problem that has been derived as a limit of a bulk-surface reaction-diffusion system which models cell polarization. In previous papers, we have established well-posedness of this problem and derived conditions on the initial data that imply continuity of the free boundary as $t\to 0$. In this paper we extend the qualitative study of the free boundary by considering axisymmetric data. Under additional monotonicity assumptions on the data we prove global continuity of the free boundary. On the other hand, if the initial data violate a "no-fattening" condition we show that the free boundary can oscillate as $t \to 0$.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
Compactness and existence theory for a general class of stationary radiative transfer equations
Authors:
Elena Demattè,
Jin Woo Jang,
Juan J. L. Velázquez
Abstract:
In this paper, we study the steady-states of a large class of stationary radiative transfer equations in a $C^1$ convex bounded domain. Namely, we consider the case in which both absorption-emission and scattering coefficients depend on the local temperature $T$ and the radiation frequency $ν.$ The radiative transfer equation determines the temperature of the material at each point. The main diffi…
▽ More
In this paper, we study the steady-states of a large class of stationary radiative transfer equations in a $C^1$ convex bounded domain. Namely, we consider the case in which both absorption-emission and scattering coefficients depend on the local temperature $T$ and the radiation frequency $ν.$ The radiative transfer equation determines the temperature of the material at each point. The main difficulty in proving existence of solutions is to obtain compactness of the sequence of integrals along lines that appear in several exponential terms. We prove a new compactness result suitable to deal with such a non-local operator containing integrals on a line segment. On the other hand, to obtain the existence theory of the full equation with both absorption and scattering terms we combine the compactness result with the construction of suitable Green functions for a class of non-local equations.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
On the diffusion approximation of the stationary radiative transfer equation with absorption and emission
Authors:
Elena Demattè,
Juan J. L. Velázquez
Abstract:
We study the situation in which the distribution of temperature a body is due to its interaction with radiation. We consider the boundary value problem for the stationary radiative transfer equation under the assumption of the local thermodynamic equilibrium. We study the diffusion equilibrium approximation in the absence of scattering. We consider absorption coefficient independent of the frequen…
▽ More
We study the situation in which the distribution of temperature a body is due to its interaction with radiation. We consider the boundary value problem for the stationary radiative transfer equation under the assumption of the local thermodynamic equilibrium. We study the diffusion equilibrium approximation in the absence of scattering. We consider absorption coefficient independent of the frequency $ ν$ (the so-called Grey approximation) and the limit when the photons' mean free path tends to zero, i.e. the absorption coefficient tends to infinity. We show that the densitive of radiative energy $ u $, which is proportional to the fourth power of the temperature due to the Stefan-Boltzmann law, solves in the limit an elliptic equation where the boundary value can be determined uniquely in terms of the original boundary condition. We derive formally with the method of matched asymptotic expansions the boundary condition for the limit problem and we prove rigorously the convergence to the solution of the limit problem with a careful analysis of some non-local integral operators. The method developed here allows to prove all the results using only maximum principle arguments.
△ Less
Submitted 20 September, 2023;
originally announced September 2023.
-
Description of Chemical Systems by means of Response Functions
Authors:
E. Franco,
B. Kepka,
J. J. L. Velázquez
Abstract:
In this paper we introduce a formalism that allows to describe the response of a part of a biochemical system in terms of renewal equations. In particular, we examine under which conditions the interactions between the different parts of a chemical system, described by means of linear ODEs, can be represented in terms of renewal equations. We show also how to apply the formalism developed in this…
▽ More
In this paper we introduce a formalism that allows to describe the response of a part of a biochemical system in terms of renewal equations. In particular, we examine under which conditions the interactions between the different parts of a chemical system, described by means of linear ODEs, can be represented in terms of renewal equations. We show also how to apply the formalism developed in this paper to some particular types of linear and non-linear ODEs, modelling some biochemical systems of interest in biology (for instance, some time-dependent versions of the classical Hopfield model of kinetic proofreading). We also analyse some of the properties of the renewal equations that we are interested in, as the long-time behaviour of their solution. Furthermore, we prove that the kernels characterising the renewal equations derived by biochemical system with reactions that satisfy the detail balance condition belong to the class of completely monotone functions.
△ Less
Submitted 5 September, 2023;
originally announced September 2023.
-
Coagulation equations with source leading to anomalous self-similarity
Authors:
Marina A. Ferreira,
Eugenia Franco,
Jani Lukkarinen,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
We study the long-time behaviour of the solutions to Smoluchowski coagulation equations with a source term of small clusters. The source drives the system out-of-equilibrium, leading to a rich range of different possible long-time behaviours, including anomalous self-similarity. The coagulation kernel is non-gelling, homogeneous, with homogeneity $γ\leq -1 $, and behaves like $x^{γ+λ} y^{-λ} $ whe…
▽ More
We study the long-time behaviour of the solutions to Smoluchowski coagulation equations with a source term of small clusters. The source drives the system out-of-equilibrium, leading to a rich range of different possible long-time behaviours, including anomalous self-similarity. The coagulation kernel is non-gelling, homogeneous, with homogeneity $γ\leq -1 $, and behaves like $x^{γ+λ} y^{-λ} $ when $y \ll x$ with $γ+2λ> 1 $. Our analysis shows that the long-time behaviour of the solutions depends on the parameters $γ$ and $λ$. More precisely, we argue that the long-time behaviour is self-similar, although the scaling of the self-similar solutions depends on the sign of $γ+λ$ and on whether $γ=-1$ or $γ< -1$. In all these cases, the scaling differs from the usual one that has been previously obtained when $γ+2λ<1$ or $γ+2λ\geq 1, γ>-1$. In the last part of the paper, we present some conjectures supporting the self-similar ansatz also for the critical case $γ+2λ=1, γ\leq -1 $.
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
Fast fusion in a two-dimensional coagulation model
Authors:
Iulia Cristian,
Juan J. L. Velázquez
Abstract:
In this work, we study a particular system of coagulation equations characterized by two values, namely volume $v$ and surface area $a$. Compared to the standard one-dimensional models, this model incorporates additional information about the geometry of the particles. We describe the coagulation process as a combination between collision and fusion of particles. We prove that we are able to recov…
▽ More
In this work, we study a particular system of coagulation equations characterized by two values, namely volume $v$ and surface area $a$. Compared to the standard one-dimensional models, this model incorporates additional information about the geometry of the particles. We describe the coagulation process as a combination between collision and fusion of particles. We prove that we are able to recover the standard one-dimensional coagulation model when fusion happens quickly and that we are able to recover an equation in which particles interact and form a ramified-like system in time when fusion happens slowly.
△ Less
Submitted 16 March, 2023;
originally announced March 2023.
-
Rotating solutions to the incompressible Euler-Poisson equation with external particle
Authors:
Diego Alonso-Orán,
Bernhard Kepka,
Juan J. L. Velázquez
Abstract:
We consider a two-dimensional, incompressible fluid body, together with self-induced interactions. The body is perturbed by an external particle with small mass. The whole configuration rotates uniformly around the common center of mass. We construct solutions, which are stationary in a rotating coordinate system, using perturbative methods. In addition, we consider a large class of internal motio…
▽ More
We consider a two-dimensional, incompressible fluid body, together with self-induced interactions. The body is perturbed by an external particle with small mass. The whole configuration rotates uniformly around the common center of mass. We construct solutions, which are stationary in a rotating coordinate system, using perturbative methods. In addition, we consider a large class of internal motions of the fluid. The angular velocity is related to the position of the external particle and is chosen to satisfy a non-resonance condition.
△ Less
Submitted 2 February, 2023;
originally announced February 2023.
-
A thin film model for meniscus evolution
Authors:
Amrita Ghosh,
Juan J. L. Velázquez
Abstract:
In this paper, we discuss a particular model arising from sinking of a rigid solid into a thin film of fluid, i.e. a fluid contained between two solid surfaces and part of the fluid surface is in contact with the air. The fluid is governed by Navier-Stokes equation, while the contact point, i.e. where the gas, liquid and solid meet, is assumed to be given by a constant, non-zero contact angle. We…
▽ More
In this paper, we discuss a particular model arising from sinking of a rigid solid into a thin film of fluid, i.e. a fluid contained between two solid surfaces and part of the fluid surface is in contact with the air. The fluid is governed by Navier-Stokes equation, while the contact point, i.e. where the gas, liquid and solid meet, is assumed to be given by a constant, non-zero contact angle. We consider a scaling limit of the fluid thickness (lubrication approximation) and the contact angle between the fluid-solid and the fluid-gas interfaces is close to $π$. This resulting model is a free boundary problem for the equation $h_t + (h^3h_{xxx})_x = 0$, for which we have $h>0$ at the contact point (different from the usual thin film equation with $h=0$ at the contact point). We show that this fourth order quasilinear (non-degenerate) parabolic equation, together with the so-called partial wetting condition at the contact point, is well-posed. Also the contact point in our thin film equation can actually move, contrary to the classical thin film equation for a droplet arising from no-slip condition. Furthermore, we show the global stability of steady state solutions in a periodic setting.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
Long-time asymptotics for coagulation equations with injection that do not have stationary solutions
Authors:
Iulia Cristian,
Marina A. Ferreira,
Eugenia Franco,
Juan J. L. Velázquez
Abstract:
In this paper we study a class of coagulation equations including a source term that injects in the system clusters of size of order one. The coagulation kernel is homogeneous, of homogeneity $γ< 1$, such that $K(x,y)$ is approximately $x^{γ+ λ} y^{-λ}$, when $x$ is larger than $y$. We restrict the analysis to the case $γ+ 2 λ\geq 1 $. In this range of exponents, the transport of mass toward infin…
▽ More
In this paper we study a class of coagulation equations including a source term that injects in the system clusters of size of order one. The coagulation kernel is homogeneous, of homogeneity $γ< 1$, such that $K(x,y)$ is approximately $x^{γ+ λ} y^{-λ}$, when $x$ is larger than $y$. We restrict the analysis to the case $γ+ 2 λ\geq 1 $. In this range of exponents, the transport of mass toward infinity is driven by collisions between particles of different sizes. This is in contrast with the case when $γ+ 2 λ<1$. In that case, the transport of mass toward infinity is due to the collision between particles of comparable sizes. In the case $γ+2λ\geq 1$, the interaction between particles of different sizes leads to an additional transport term in the coagulation equation that approximates the solution of the original coagulation equation with injection for large times. We prove the existence of a class of self-similar solutions for suitable choices of $γ$ and $λ$ for this class of coagulation equations with transport. We prove that for the complementary case such self-similar solutions do not exist.
△ Less
Submitted 29 November, 2022;
originally announced November 2022.
-
On the temperature distribution of a body heated by radiation
Authors:
Jin Woo Jang,
Juan J. L. Velázquez
Abstract:
In this paper, we study the temperature distribution of a body when the heat is transmitted only by radiation. The heat transmitted by convection and conduction is ignored. We consider the stationary radiative transfer equation in the local thermodynamic equilibrium. We prove that the stationary radiative transfer equation coupled with the non-local temperature equation is well-posed in a generic…
▽ More
In this paper, we study the temperature distribution of a body when the heat is transmitted only by radiation. The heat transmitted by convection and conduction is ignored. We consider the stationary radiative transfer equation in the local thermodynamic equilibrium. We prove that the stationary radiative transfer equation coupled with the non-local temperature equation is well-posed in a generic case when emission-absorption or scattering of interacting radiation is considered. The emission-absorption and the scattering coefficients are assumed to be general and they can depend on the frequency of radiation. We also establish an entropy production formula of the system, which is used to prove the uniqueness of solutions for an incoming radiation with constant temperature.
△ Less
Submitted 20 June, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
Vanishing angular singularity limit to the hard-sphere Boltzmann equation
Authors:
Jin Woo Jang,
Bernhard Kepka,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
In this note we study Boltzmann's collision kernel for inverse power law interactions $U_s(r)=1/r^{s-1}$ for $s>2$ in dimension $ d=3 $. We prove the limit of the non-cutoff kernel to the hard-sphere kernel and give precise asymptotic formulas of the singular layer near $θ\simeq 0$ in the limit $ s\to \infty $. Consequently, we show that solutions to the homogeneous Boltzmann equation converge to…
▽ More
In this note we study Boltzmann's collision kernel for inverse power law interactions $U_s(r)=1/r^{s-1}$ for $s>2$ in dimension $ d=3 $. We prove the limit of the non-cutoff kernel to the hard-sphere kernel and give precise asymptotic formulas of the singular layer near $θ\simeq 0$ in the limit $ s\to \infty $. Consequently, we show that solutions to the homogeneous Boltzmann equation converge to the respective solutions.
△ Less
Submitted 28 September, 2022;
originally announced September 2022.
-
Coagulation equations for non-spherical clusters
Authors:
Iulia Cristian,
Juan J. L. Velázquez
Abstract:
In this work, we study the long time asymptotics of a coagulation model which describes the evolution of a system of particles characterized by their volume and surface area. The aggregation mechanism takes place in two stages: collision and fusion of particles. During the collision stage, the two particles merge at a contact point. The newly formed particle has volume and area equal to the sum of…
▽ More
In this work, we study the long time asymptotics of a coagulation model which describes the evolution of a system of particles characterized by their volume and surface area. The aggregation mechanism takes place in two stages: collision and fusion of particles. During the collision stage, the two particles merge at a contact point. The newly formed particle has volume and area equal to the sum of the respective quantities of the two colliding particles. After collision, the fusion phase begins and during it the geometry of the interacting particles is modified in such a way that the volume of the total system is preserved and the surface area is reduced. During their evolution, the particles must satisfy the isoperimetric inequality. Therefore, the distribution of particles in the volume and area space is supported in the region where $\{a\geq (36π)^{\frac{1}{3}}v^{\frac{2}{3}}\}$. We assume the coagulation kernel has a weak dependence on the area variable. We prove existence of self-similar profiles for some choices of the functions describing the fusion rate for which the particles have a shape that is close to spherical. On the other hand, for other fusion mechanisms and suitable choices of initial data, we show that the particle distribution describes a system of ramified-like particles.
△ Less
Submitted 10 July, 2023; v1 submitted 1 September, 2022;
originally announced September 2022.
-
Non-power law constant flux solutions for the Smoluchowski coagulation equation
Authors:
Marina A. Ferreira,
Jani Lukkarinen,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
It is well known that for a large class of coagulation kernels, Smoluchowski coagulation equations have particular power law solutions which yield a constant flux of mass along all scales of the system. In this paper, we prove that for some choices of the coagulation kernels there are solutions with a constant flux of mass along all scales which are not power laws. The result is proved by means of…
▽ More
It is well known that for a large class of coagulation kernels, Smoluchowski coagulation equations have particular power law solutions which yield a constant flux of mass along all scales of the system. In this paper, we prove that for some choices of the coagulation kernels there are solutions with a constant flux of mass along all scales which are not power laws. The result is proved by means of a bifurcation argument.
△ Less
Submitted 24 July, 2022; v1 submitted 19 July, 2022;
originally announced July 2022.
-
Asymptotic localization in multicomponent mass conserving coagulation equations
Authors:
Marina A. Ferreira,
Jani Lukkarinen,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
In this paper we prove that the time dependent solutions of a large class of Smoluchowski coagulation equations for multicomponent systems concentrate along a particular direction of the space of cluster compositions for long times. The direction of concentration is determined by the initial distribution of clusters. These results allow to prove the uniqueness and global stability of the self-simi…
▽ More
In this paper we prove that the time dependent solutions of a large class of Smoluchowski coagulation equations for multicomponent systems concentrate along a particular direction of the space of cluster compositions for long times. The direction of concentration is determined by the initial distribution of clusters. These results allow to prove the uniqueness and global stability of the self-similar profile with finite mass in the case of coagulation kernels which are not identically constant, but are constant along any direction of the space of cluster compositions.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
On the Grad-Rubin boundary value problem for the two-dimensional magneto-hydrostatic equations
Authors:
Diego Alonso-Orán,
Juan J. L. Velázquez
Abstract:
In this work, we study the solvability of a boundary value problem for the magneto-hydrostatic equations originally proposed by Grad and Rubin in the late 50's. The proof relies on a fixed point argument which combines the so-called current transport method together with Hölder estimates for a class of non-convolution singular integral operators. The same method allows to solve an analogous bounda…
▽ More
In this work, we study the solvability of a boundary value problem for the magneto-hydrostatic equations originally proposed by Grad and Rubin in the late 50's. The proof relies on a fixed point argument which combines the so-called current transport method together with Hölder estimates for a class of non-convolution singular integral operators. The same method allows to solve an analogous boundary value problem for the steady incompressible Euler equations.
△ Less
Submitted 8 March, 2022;
originally announced March 2022.
-
Long-time asymptotics and regularity estimates for weak solutions to a doubly degenerate thin-film equation in the Taylor-Couette setting
Authors:
Christina Lienstromberg,
Juan J. L. Velázquez
Abstract:
We study the dynamic behaviour of solutions to a fourth-order quasilinear degenerate parabolic equation for large times arising in fluid dynamical applications. The degeneracy occurs both with respect to the unknown and with respect to the sum of its first and third spatial derivative. The modelling equation appears as a thin-film limit for the interface separating two immiscible viscous fluid fil…
▽ More
We study the dynamic behaviour of solutions to a fourth-order quasilinear degenerate parabolic equation for large times arising in fluid dynamical applications. The degeneracy occurs both with respect to the unknown and with respect to the sum of its first and third spatial derivative. The modelling equation appears as a thin-film limit for the interface separating two immiscible viscous fluid films confined between two cylinders rotating at small relative angular velocity. More precisely, the fluid occupying the layer next to the outer cylinder is considered to be Newtonian, i.e. it has constant viscosity, while we assume that the layer next to the inner cylinder is filled by a shear-thinning power-law fluid.
Using energy methods, Fourier analysis and suitable regularity estimates for higher-order parabolic equations, we prove global existence of positive weak solutions in the case of low initial energy. Moreover, these global solutions are polynomially stable, in the sense that interfaces which are initially close to a circle, converge at rate $1/t^{1/β}$ for some $β> 0$ to a circle, as time tends to infinity.
In addition, we provide regularity estimates for general nonlinear degenerate parabolic equations of fourth order.
△ Less
Submitted 28 February, 2022;
originally announced March 2022.
-
Qualitative properties of solutions to a non-local free boundary problem modeling cell polarization
Authors:
Anna Logioti,
Barbara Niethammer,
Matthias Röger,
Juan J. L. Velázquez
Abstract:
We consider a parabolic non-local free boundary problem that has been derived as a limit of a bulk-surface reaction-diffusion system which models cell polarization. The authors have justified the well-posedness of this problem and have further proved uniqueness of solutions and global stability of steady states. In this paper we investigate qualitative properties of the free boundary. We present n…
▽ More
We consider a parabolic non-local free boundary problem that has been derived as a limit of a bulk-surface reaction-diffusion system which models cell polarization. The authors have justified the well-posedness of this problem and have further proved uniqueness of solutions and global stability of steady states. In this paper we investigate qualitative properties of the free boundary. We present necessary and sufficient conditions for the initial data that imply continuity of the support at $t = 0$. If one of these assumptions fail, then jumps of the support take place. In addition we provide a complete characterization of the jumps for a large class of initial data.
△ Less
Submitted 21 April, 2023; v1 submitted 13 February, 2022;
originally announced February 2022.
-
Revisiting Shikhmurzaev's approach to the contact line problem
Authors:
Amrita Ghosh,
Barbara Niethammer,
Juan J. L. Velazquez
Abstract:
In this paper, we revisit a model for the contact line problem which has been proposed by Shikhmurzaev (1993). In the first part, in addition to rederiving the model, we study in detail the assumptions required to obtain the isothermal limit of the model. We also derive in this paper several lubrication approximation models, based on Shikhmurzaev's approach. The first two lubrication models descri…
▽ More
In this paper, we revisit a model for the contact line problem which has been proposed by Shikhmurzaev (1993). In the first part, in addition to rederiving the model, we study in detail the assumptions required to obtain the isothermal limit of the model. We also derive in this paper several lubrication approximation models, based on Shikhmurzaev's approach. The first two lubrication models describe thin film flow of incompressible fluids on solid substrates, based on different orders of magnitude of the slip length parameter. The third lubrication model describes a meniscus formation where a wedge-shaped solid immerses in a thin film of fluid.
△ Less
Submitted 11 August, 2022; v1 submitted 23 December, 2021;
originally announced December 2021.
-
Some thin film models based on Shikhmurzaev's approach
Authors:
Amrita Ghosh,
Barbara Niethammer,
Juan J. L. Velazquez
Abstract:
We derive several lubrication approximation models, using Shikhmurzaev's approach to the contact line problem, obtained in \cite{GNV}. The first two lubrication models describe thin film flow of incompressible fluids on solid substrates, based on different orders of magnitude of the slip length parameter. The third lubrication model describes a meniscus formation where a wedge-shaped solid immerse…
▽ More
We derive several lubrication approximation models, using Shikhmurzaev's approach to the contact line problem, obtained in \cite{GNV}. The first two lubrication models describe thin film flow of incompressible fluids on solid substrates, based on different orders of magnitude of the slip length parameter. The third lubrication model describes a meniscus formation where a wedge-shaped solid immerses in a thin film of fluid.
△ Less
Submitted 29 December, 2021; v1 submitted 22 December, 2021;
originally announced December 2021.
-
LTE and Non-LTE Solutions in Gases Interacting with Radiation
Authors:
Jin Woo Jang,
Juan J. L. Velázquez
Abstract:
In this paper, we study a class of kinetic equations describing radiative transfer in gases which include also the interaction of the molecules of the gas with themselves. We discuss several scaling limits and introduce some Euler-like systems coupled with radiation as an aftermath of specific scaling limits. We consider scaling limits in which local thermodynamic equilibrium (LTE) holds, as well…
▽ More
In this paper, we study a class of kinetic equations describing radiative transfer in gases which include also the interaction of the molecules of the gas with themselves. We discuss several scaling limits and introduce some Euler-like systems coupled with radiation as an aftermath of specific scaling limits. We consider scaling limits in which local thermodynamic equilibrium (LTE) holds, as well as situations in which this assumption fails (non-LTE). The structure of the equations describing the gas-radiation system is very different in the LTE and non-LTE cases. We prove the existence of stationary solutions with zero velocities to the resulting limit models in the LTE case. We also prove the non-existence of a stationary state with zero velocities in a non-LTE case.
△ Less
Submitted 7 December, 2021; v1 submitted 21 September, 2021;
originally announced September 2021.
-
Edge States for generalised Iwatsuka models: Magnetic fields having a fast transition across a curve
Authors:
Arianna Giunti,
Juan J. L. Velázquez
Abstract:
In this paper, we study the localization and propagation properties of the edge states associated to a class of magnetic laplacians in $\mathbb{R}^2$. We assume that the intensity of the magnetic field has a fast transition along a regular and compact curve $Γ$. Our main results extend to a general regular curve the study of the localised eigenfunction obtained when $Γ$ is a straight line (i.e. Iw…
▽ More
In this paper, we study the localization and propagation properties of the edge states associated to a class of magnetic laplacians in $\mathbb{R}^2$. We assume that the intensity of the magnetic field has a fast transition along a regular and compact curve $Γ$. Our main results extend to a general regular curve the study of the localised eigenfunction obtained when $Γ$ is a straight line (i.e. Iwatsuka models). Furthermore, we include in our analysis the case of magnetic fields that slowly change along the curve $Γ$ and we obtain a rigorous and explicit characterization of the asymptotic mass distribution of the edge state along $Γ$.
△ Less
Submitted 20 September, 2021;
originally announced September 2021.
-
On the self-similar behaviour of coagulation systems with injection
Authors:
Marina A. Ferreira,
Eugenia Franco,
Juan J. L. Velázquez
Abstract:
In this paper we prove the existence of a family of self-similar solutions for a class of coagulation equations with a constant flux of particles from the origin. These solutions are expected to describe the longtime asymptotics of Smoluchowski's coagulation equations with a time independent source of clusters concentrated in small sizes. The self-similar profiles are shown to be smooth, provided…
▽ More
In this paper we prove the existence of a family of self-similar solutions for a class of coagulation equations with a constant flux of particles from the origin. These solutions are expected to describe the longtime asymptotics of Smoluchowski's coagulation equations with a time independent source of clusters concentrated in small sizes. The self-similar profiles are shown to be smooth, provided the coagulation kernel is also smooth. Moreover, the self-similar profiles are estimated from above and from below by $x^{-(γ+3)/2}$ as $x \to 0$, where $γ<1 $ is the homogeneity of the kernel, and are proven to decay at least exponentially as $x \to \infty$.
△ Less
Submitted 31 May, 2022; v1 submitted 23 June, 2021;
originally announced June 2021.
-
Multicomponent coagulation systems: existence and non-existence of stationary non-equilibrium solutions
Authors:
Marina A. Ferreira,
Jani Lukkarinen,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
We study multicomponent coagulation via the Smoluchowski coagulation equation under non-equilibrium stationary conditions induced by a source of small clusters. The coagulation kernel can be very general, merely satisfying certain power law asymptotic bounds in terms of the total number of monomers in a cluster. The bounds are characterized by two parameters and we extend previous results for one-…
▽ More
We study multicomponent coagulation via the Smoluchowski coagulation equation under non-equilibrium stationary conditions induced by a source of small clusters. The coagulation kernel can be very general, merely satisfying certain power law asymptotic bounds in terms of the total number of monomers in a cluster. The bounds are characterized by two parameters and we extend previous results for one-component systems to classify the parameter values for which the above stationary solutions do or do not exist. Moreover, we also obtain criteria for the existence or non-existence of solutions which yield a constant flux of mass towards large clusters.
△ Less
Submitted 18 March, 2021;
originally announced March 2021.
-
Oscillations in a Becker-Döring model with injection and depletion
Authors:
Barbara Niethammer,
Robert L. Pego,
André Schlichting,
Juan J. L. Velázquez
Abstract:
We study the Becker-Döring bubblelator, a variant of the Becker-Döring coagulation-fragmentation system that models the growth of clusters by gain or loss of monomers. Motivated by models of gas evolution oscillators from physical chemistry, we incorporate injection of monomers and depletion of large clusters. For a wide range of physical rates, the Becker-Döring system itself exhibits a dynamic p…
▽ More
We study the Becker-Döring bubblelator, a variant of the Becker-Döring coagulation-fragmentation system that models the growth of clusters by gain or loss of monomers. Motivated by models of gas evolution oscillators from physical chemistry, we incorporate injection of monomers and depletion of large clusters. For a wide range of physical rates, the Becker-Döring system itself exhibits a dynamic phase transition as mass density increases past a critical value. We connect the Becker-Döring bubblelator to a transport equation coupled with an integrodifferential equation for excess monomer density by formal asymptotics in the near-critical regime. For suitable injection/depletion rates, we argue that time-periodic solutions appear via a Hopf bifurcation. Numerics confirm that the generation and removal of large clusters can become desynchronized, leading to temporal oscillations associated with bursts of large-cluster nucleation.
△ Less
Submitted 27 September, 2021; v1 submitted 12 February, 2021;
originally announced February 2021.
-
Boundary value problems for two dimensional steady incompressible fluids
Authors:
Diego Alonso-Orán,
Juan Juan J. L. Velázquez
Abstract:
In this paper we study the solvability of different boundary value problems for the two dimensional steady incompressible Euler equation. Two main methods are currently available to study those problems, namely the Grad-Shafranov method and the vorticity transport method. We describe for which boundary value problems these methods can be applied. The obtained solutions have non-vanishing vorticity…
▽ More
In this paper we study the solvability of different boundary value problems for the two dimensional steady incompressible Euler equation. Two main methods are currently available to study those problems, namely the Grad-Shafranov method and the vorticity transport method. We describe for which boundary value problems these methods can be applied. The obtained solutions have non-vanishing vorticity.
△ Less
Submitted 18 January, 2021;
originally announced January 2021.
-
Analysis of a two-fluid Taylor-Couette flow with one non-Newtonian fluid
Authors:
Christina Lienstromberg,
Tania Pernas-Castaño,
Juan J. L. Velázquez
Abstract:
We study the dynamic behaviour of two viscous fluid films confined between two concentric cylinders rotating at a small relative velocity. It is assumed that the fluids are immiscible and that the volume of the outer fluid film is large compared to the volume of the inner one. Moreover, while the outer fluid is considered to have constant viscosity, the rheological behaviour of the inner thin film…
▽ More
We study the dynamic behaviour of two viscous fluid films confined between two concentric cylinders rotating at a small relative velocity. It is assumed that the fluids are immiscible and that the volume of the outer fluid film is large compared to the volume of the inner one. Moreover, while the outer fluid is considered to have constant viscosity, the rheological behaviour of the inner thin film is determined by a strain-dependent power-law. Starting from a Navier--Stokes system, we formally derive evolution equations for the interface separating the two fluids. Two competing effects drive the dynamics of the interface, namely, the surface tension and the shear stresses induced by the rotation of the cylinders. When the two effects are comparable, the solutions behave, for large times, as in the Newtonian regime. We also study the regime in which the surface tension effects dominate the stresses induced by the rotation of the cylinders. In this case, we prove local existence of positive weak solutions both for shear-thinning and shear-thickening fluids. In the latter case, we show that interfaces which are initially close to a circle converge to a circle in finite time and keep that shape for later times.
△ Less
Submitted 30 September, 2021; v1 submitted 19 December, 2020;
originally announced December 2020.
-
Kinetic Models for Semiflexible Polymers in a Half-plane
Authors:
Jin Woo Jang,
Juan J. L. Velázquez
Abstract:
Based on a general discrete model for a semiflexible polymer chain, we introduce a formal derivation of a kinetic equation for semiflexible polymers in the half-plane via a continuum limit. It turns out that the resulting equation is the kinetic Fokker-Planck-type equation with the Laplace-Beltrami operator under a non-local trapping boundary condition. We then study the well-posedness and the lon…
▽ More
Based on a general discrete model for a semiflexible polymer chain, we introduce a formal derivation of a kinetic equation for semiflexible polymers in the half-plane via a continuum limit. It turns out that the resulting equation is the kinetic Fokker-Planck-type equation with the Laplace-Beltrami operator under a non-local trapping boundary condition. We then study the well-posedness and the long-chain asymptotics of the solutions of the resulting equation. In particular, we prove that there exists a unique measure-valued solution for the corresponding boundary value problem. In addition, we prove that the equation is hypoelliptic and the solutions are locally Hölder continuous near the singular boundary. Finally, we provide the asymptotic behaviors of the solutions for large polymer chains.
△ Less
Submitted 28 July, 2022; v1 submitted 3 November, 2020;
originally announced November 2020.
-
On the dynamics of thin layers of viscous flows inside another viscous fluid
Authors:
Tania Pernas-Castaño,
Juan J. L. Velázquez
Abstract:
In this work we will study the dynamics of a thin layer of a viscous fluid which is embedded in the interior of another viscous fluid. The resulting flow can be approximated by means of the solutions of a free boundary problem for the Stokes equation in which one of the unknowns is the shape of a curve which approximates the geometry of the thin layer of fluid. We also derive the equation yielding…
▽ More
In this work we will study the dynamics of a thin layer of a viscous fluid which is embedded in the interior of another viscous fluid. The resulting flow can be approximated by means of the solutions of a free boundary problem for the Stokes equation in which one of the unknowns is the shape of a curve which approximates the geometry of the thin layer of fluid. We also derive the equation yielding the thickness of this fluid. This model, that will be termed as the Geometric Free Boundary Problem, will be derived using matched asymptotic expansions. We will prove that the Geometric Free Boundary Problem is well posed and the solutions of the thickness equation are well defined (in particular they do not yield breaking of fluid layers) as long as the solutions of the Geometric Free Boundary Problem exist.
△ Less
Submitted 28 October, 2020;
originally announced October 2020.
-
A parabolic free boundary problem arising in a model of cell polarization
Authors:
Anna Logioti,
Barbara Niethammer,
Matthias Röger,
Juan J. L. Velázquez
Abstract:
The amplification of an external signal is a key step in direction sensing of biological cells. We consider a simple model for the response to a time-depending signal, which was previously proposed by the last three authors. The model consists of a bulk-surface reaction-diffusion model. We prove that in a suitable asymptotic limit the system converges to a bulk-surface parabolic obstacle type prob…
▽ More
The amplification of an external signal is a key step in direction sensing of biological cells. We consider a simple model for the response to a time-depending signal, which was previously proposed by the last three authors. The model consists of a bulk-surface reaction-diffusion model. We prove that in a suitable asymptotic limit the system converges to a bulk-surface parabolic obstacle type problem. For this model and a reduction to a nonlocal surface equation we show an L1-contraction property and, in the case of time-constant signals, the stability of stationary states.
△ Less
Submitted 29 June, 2020;
originally announced June 2020.
-
Localization in stationary non-equilibrium solutions for multicomponent coagulation systems
Authors:
Marina A. Ferreira,
Jani Lukkarinen,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
We consider the multicomponent Smoluchowski coagulation equation under non-equilibrium conditions induced either by a source term or via a constant flux constraint. We prove that the corresponding stationary non-equilibrium solutions have a universal localization property. More precisely, we show that these solutions asymptotically localize into a direction determined by the source or by a flux co…
▽ More
We consider the multicomponent Smoluchowski coagulation equation under non-equilibrium conditions induced either by a source term or via a constant flux constraint. We prove that the corresponding stationary non-equilibrium solutions have a universal localization property. More precisely, we show that these solutions asymptotically localize into a direction determined by the source or by a flux constraint: the ratio between monomers of a given type to the total number of monomers in the cluster becomes ever closer to a predetermined ratio as the cluster size is increased. The assumptions on the coagulation kernel are quite general, with isotropic power law bounds. The proof relies on a particular measure concentration estimate and on the control of asymptotic scaling of the solutions which is allowed by previously derived estimates on the mass current observable of the system.
△ Less
Submitted 23 March, 2021; v1 submitted 26 June, 2020;
originally announced June 2020.
-
Self-similar asymptotics for a modified Maxwell-Boltzmann equation in systems subject to deformations
Authors:
Alexander Bobylev,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
In this paper we study a generalized class of Maxwell-Boltzmann equations which in addition to the usual collision term contains a linear deformation term described by a matrix A. This class of equations arises, for instance, from the analysis of homoenergetic solutions for the Boltzmann equation considered by many authors since 1950s. Our main goal is to study a large time asymptotics of solution…
▽ More
In this paper we study a generalized class of Maxwell-Boltzmann equations which in addition to the usual collision term contains a linear deformation term described by a matrix A. This class of equations arises, for instance, from the analysis of homoenergetic solutions for the Boltzmann equation considered by many authors since 1950s. Our main goal is to study a large time asymptotics of solutions under assumption of smallness of the matrix A. The main result of the paper is formulated in Theorem 2.1. Informally stated, this Theorem says that, for sufficiently small norm of A, any non-negative solution with finite second moment tends to a self-similar solution of relatively simple form for large values of time. This is what we call "the self-similar asymptotics". We also prove that the higher order moments of the self-similar profile are finite under further smallness condition on the matrix A.
△ Less
Submitted 28 December, 2019;
originally announced December 2019.
-
Edge States for the magnetic Laplacian in domains with smooth boundary
Authors:
Arianna Giunti,
Juan J. L. Velázquez
Abstract:
We are interested in the spectral properties of the magnetic Schrödinger operator $H_\varepsilon$ in a domain $Ω\subset \mathbb{R}^2$ with compact boundary and with magnetic field of intensity $\varepsilon^{-2}$. We impose Dirichlet boundary conditions on $\partialΩ$. Our main focus is the existence and description of the so-called \textit{edge states}, namely eigenfunctions for $H_{\varepsilon}$…
▽ More
We are interested in the spectral properties of the magnetic Schrödinger operator $H_\varepsilon$ in a domain $Ω\subset \mathbb{R}^2$ with compact boundary and with magnetic field of intensity $\varepsilon^{-2}$. We impose Dirichlet boundary conditions on $\partialΩ$. Our main focus is the existence and description of the so-called \textit{edge states}, namely eigenfunctions for $H_{\varepsilon}$ whose mass is localized at scale $\varepsilon$ along the boundary $\partialΩ$. When the intensity of the magnetic field is large (i.e. $\varepsilon <<1$), we show that such edge states exist. Furthermore, we give a detailed description of their localization close to the boundary $\partialΩ$, as well as how their mass is distributed along it. From this result, we also infer asymptotic formulas for the eigenvalues of $H_\varepsilon$.
△ Less
Submitted 4 September, 2020; v1 submitted 16 December, 2019;
originally announced December 2019.
-
Stationary non-equilibrium solutions for coagulation systems
Authors:
Marina A. Ferreira,
Jani Lukkarinen,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
We study coagulation equations under non-equilibrium conditions which are induced by the addition of a source term for small cluster sizes. We consider both discrete and continuous coagulation equations, and allow for a large class of coagulation rate kernels, with the main restriction being boundedness from above and below by certain weight functions. The weight functions depend on two power law…
▽ More
We study coagulation equations under non-equilibrium conditions which are induced by the addition of a source term for small cluster sizes. We consider both discrete and continuous coagulation equations, and allow for a large class of coagulation rate kernels, with the main restriction being boundedness from above and below by certain weight functions. The weight functions depend on two power law parameters, and the assumptions cover, in particular, the commonly used free molecular and diffusion limited aggregation coagulation kernels. Our main result shows that the two weight function parameters already determine whether there exists a stationary solution under the presence of a source term. In particular, we find that the diffusive kernel allows for the existence of stationary solutions while there cannot be any such solutions for the free molecular kernel. The argument to prove the non-existence of solutions relies on a novel power law lower bound, valid in the appropriate parameter regime, for the decay of stationary solutions with a constant flux. We obtain optimal lower and upper estimates of the solutions for large cluster sizes, and prove that the solutions of the discrete model behave asymptotically as solutions of the continuous model.
△ Less
Submitted 25 September, 2020; v1 submitted 23 September, 2019;
originally announced September 2019.
-
Two-speed solutions to non-convex rate-independent systems
Authors:
Filip Rindler,
Sebastian Schwarzacher,
Juan J. L. Velazquez
Abstract:
We consider evolutionary PDE inclusions of the form \[
-λ\dot{u}_λ+ Δu - \mathrm{D} W_0(u) + f \ni \partial \mathcal{R}_1(\dot{u}) \qquad\text{in $(0,T) \times Ω$,} \] where $\mathcal{R}_1$ is a positively $1$-homogeneous rate-independent dissipation potential and $W_0$ is a (generally) non-convex energy density. This work constructs solutions to the above system in the slow-loading limit…
▽ More
We consider evolutionary PDE inclusions of the form \[
-λ\dot{u}_λ+ Δu - \mathrm{D} W_0(u) + f \ni \partial \mathcal{R}_1(\dot{u}) \qquad\text{in $(0,T) \times Ω$,} \] where $\mathcal{R}_1$ is a positively $1$-homogeneous rate-independent dissipation potential and $W_0$ is a (generally) non-convex energy density. This work constructs solutions to the above system in the slow-loading limit $λ\downarrow 0$. Our solutions have more regularity both in space and time than those that have been obtained with other approaches. On the "slow" time scale we see strong solutions to a purely rate-independent evolution. Over the jumps, we obtain a detailed description of the behavior of the solution and we resolve the jump transients at a "fast" time scale, where the original rate-dependent evolution is still visible. Crucially, every jump transient splits into a (possibly countable) number of rate-dependent evolutions, for which the energy dissipation can be explicitly computed. This in particular yields a global energy equality for the whole evolution process. It also turns out that there is a canonical slow time scale that avoids intermediate-scale effects, where movement occurs in a mixed rate-dependent / rate-independent way. In this way, we obtain precise information on the impact of the approximation on the constructed solution. Our results are illustrated by examples, which elucidate the effects that can occur.
△ Less
Submitted 21 December, 2019; v1 submitted 11 July, 2019;
originally announced July 2019.
-
Long time asymptotics for homoenergetic solutions of the Boltzmann equation. Hyperbolic-dominated case
Authors:
Richard D. James,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
In this paper we continue the formal analysis of the long-time asymptotics of the homoenergetic solutions for the Boltzmann equation that we began in [18]. They have the form $f\left( x,v,t\right) =g\left(v-L\left( t\right) x,t\right) $ where $L\left( t\right) =A\left(I+tA\right) ^{-1}$ where $A$ is a constant matrix. Homoenergetic solutions satisfy an integro-differential equation which contains,…
▽ More
In this paper we continue the formal analysis of the long-time asymptotics of the homoenergetic solutions for the Boltzmann equation that we began in [18]. They have the form $f\left( x,v,t\right) =g\left(v-L\left( t\right) x,t\right) $ where $L\left( t\right) =A\left(I+tA\right) ^{-1}$ where $A$ is a constant matrix. Homoenergetic solutions satisfy an integro-differential equation which contains, in addition to the classical Boltzmann collision operator, a linear hyperbolic term. Depending on the properties of the collision kernel the collision and the hyperbolic terms might be of the same order of magnitude as $t\to\infty$, or the collision term could be the dominant one for large times, or the hyperbolic term could be the largest. The first case has been rigorously studied in [17]. Formal asymptotic expansions in the second case have been obtained in [18]. All the solutions obtained in this case can be approximated by Maxwellian distributions with changing temperature. In this paper we focus in the case where the hyperbolic terms are much larger than the collision term for large times (hyperbolic-dominated behavior). In the hyperbolic-dominated case it does not seem to be possible to describe in a simple way all the long time asymptotics of the solutions, but we discuss several physical situations and formulate precise conjectures. We give explicit formulas for the relationship between density, temperature and entropy for these solutions. These formulas differ greatly from the ones at equilibrium.
△ Less
Submitted 20 June, 2019;
originally announced June 2019.
-
Analysis of a thin film approximation for two-fluid Taylord-Couette flows
Authors:
Tania Pernas-Castaño,
Juan J. L. Velázquez
Abstract:
In this work we study the evolution of the interface between two different fluids in two concentric cylinders when the velocity is given by the Navier-Stokes equation and one of the fluids is thin. We present a formal asymptotic derivation of the evolution equation for the interface under different scaling assumptions for the surface tension. We then study the different types of the stationary sol…
▽ More
In this work we study the evolution of the interface between two different fluids in two concentric cylinders when the velocity is given by the Navier-Stokes equation and one of the fluids is thin. We present a formal asymptotic derivation of the evolution equation for the interface under different scaling assumptions for the surface tension. We then study the different types of the stationary solutions and travelling waves for the resulting equation. In particular, we state a global well posedness result and using Center Manifold Theory, we obtain detailed information about the long time asymptotics of the solutions of the problem.
△ Less
Submitted 31 May, 2019;
originally announced May 2019.
-
Temporal oscillations in Becker-Doering equations with atomization
Authors:
Robert L. Pego,
Juan J. L. Velázquez
Abstract:
We prove that time-periodic solutions arise via Hopf bifurcation in a finite closed system of coagulation-fragmentation equations. The system we treat is a variant of the Becker-Doering equations, in which clusters grow or shrink by addition or deletion of monomers. To this is added a linear atomization reaction for clusters of maximum size. The structure of the system is motivated by models of ga…
▽ More
We prove that time-periodic solutions arise via Hopf bifurcation in a finite closed system of coagulation-fragmentation equations. The system we treat is a variant of the Becker-Doering equations, in which clusters grow or shrink by addition or deletion of monomers. To this is added a linear atomization reaction for clusters of maximum size. The structure of the system is motivated by models of gas evolution oscillators in physical chemistry, which exhibit temporal oscillations under certain input/output conditions.
△ Less
Submitted 9 December, 2019; v1 submitted 7 May, 2019;
originally announced May 2019.
-
A bulk-surface reaction-diffusion system for cell polarization
Authors:
Barbara Niethammer,
Matthias Röger,
Juan J. L. Velázquez
Abstract:
We propose a model for cell polarization as a response to an external signal which results in a system of PDEs for different variants of a protein on the cell surface and interior respectively. We study stationary states of this model in certain parameter regimes in which several reaction rates on the membrane as well as the diffusion coefficient within the cell are large. It turns out that in sui…
▽ More
We propose a model for cell polarization as a response to an external signal which results in a system of PDEs for different variants of a protein on the cell surface and interior respectively. We study stationary states of this model in certain parameter regimes in which several reaction rates on the membrane as well as the diffusion coefficient within the cell are large. It turns out that in suitable scaling limits steady states converge to solutions of some obstacle type problems. For these limiting problems we prove the onset of polarization if the total mass of protein is sufficiently small. For some variants we can even characterize precisely the critical mass for which polarization occurs.
△ Less
Submitted 15 February, 2019;
originally announced February 2019.
-
On the homogenization of random stationary elliptic operators in divergence form
Authors:
Arianna Giunti,
Juan J. L. Velázquez
Abstract:
In this note we comment on the homogenization of a random elliptic operator in divergence form $-\nabla \cdot a\nabla$, where the coefficient field $a$ is distributed according to a stationary, but not necessarily ergodic, probability measure $P$. We generalize the well-known case for $P$ stationary and ergodic by showing that the operator $-\nabla \cdot a(\frac{\cdot}{\varepsilon})\nabla$ almost…
▽ More
In this note we comment on the homogenization of a random elliptic operator in divergence form $-\nabla \cdot a\nabla$, where the coefficient field $a$ is distributed according to a stationary, but not necessarily ergodic, probability measure $P$. We generalize the well-known case for $P$ stationary and ergodic by showing that the operator $-\nabla \cdot a(\frac{\cdot}{\varepsilon})\nabla$ almost surely homogenizes to a constant-coefficient, random operator $-\nabla \cdot A_h\nabla$. Furthermore, we use a disintegration formula for $P$ with respect to a family of ergodic and stationary probability measures to show that the law of $A_h$ may be obtained by using the standard homogenization results on each probability measure of the previous family. We finally provide a more explicit formula for $A_h$ in the case of coefficient fields which are a function of a stationary Gaussian field.
△ Less
Submitted 17 September, 2018;
originally announced September 2018.
-
Long time asymptotics for homoenergetic solutions of the Boltzmann equation. Collision-dominated case
Authors:
Richard D. James,
Alessia Nota,
Juan J. L. Velázquez
Abstract:
In this paper we present a formal analysis of the long-time asymptotics of a particular class of solutions of the Boltzmann equation, known as homoenergetic solutions, which have the form $f\left( x,v,t\right)=g\left( v-L\left( t\right) x,t\right)$ where $L\left( t\right) =A\left(I+tA\right) ^{-1}$ with the matrix $A$ describing a shear flow or a dilatation or a combination of both. We began this…
▽ More
In this paper we present a formal analysis of the long-time asymptotics of a particular class of solutions of the Boltzmann equation, known as homoenergetic solutions, which have the form $f\left( x,v,t\right)=g\left( v-L\left( t\right) x,t\right)$ where $L\left( t\right) =A\left(I+tA\right) ^{-1}$ with the matrix $A$ describing a shear flow or a dilatation or a combination of both. We began this study in \cite{JNV1}. Homoenergetic solutions satisfy an integro-differential equation which contains, in addition to the classical Boltzmann collision operator, a linear hyperbolic term. In \cite{JNV1} it has been proved rigorously the existence of self-similar solutions which describe the change of the average energy of the particles of the system in the case in which there is a balance between the hyperbolic and the collision term. In this paper we focus in homoenergetic solutions for which the collision term is much larger than the hyperbolic term (collision-dominated behavior). In this case the long time asymptotics for the distribution of velocities is given by a time dependent Maxwellian distribution with changing temperature.
△ Less
Submitted 20 August, 2018;
originally announced August 2018.
-
Hopf Bifurcation in Structural Population Models
Authors:
Narek Hovsepyan,
Juan J. L. Velázquez
Abstract:
We study a nonlinear PDE problem motivated by the peculiar patterns arising in myxobacteria, namely counter-migrating cell density waves. We rigorously prove the existence of Hopf bifurcations for some specific values of the parameters of the system. This shows the existence of periodic solutions for the systems under consideration.
We study a nonlinear PDE problem motivated by the peculiar patterns arising in myxobacteria, namely counter-migrating cell density waves. We rigorously prove the existence of Hopf bifurcations for some specific values of the parameters of the system. This shows the existence of periodic solutions for the systems under consideration.
△ Less
Submitted 24 July, 2018;
originally announced July 2018.
-
Self-similar asymptotic behavior for the solutions of a linear coagulation equation
Authors:
Barbara Niethammer,
Alessia Nota,
Sebastian Throm,
Juan J. L. Velázquez
Abstract:
In this paper we consider the long time asymptotics of a linear version of the Smoluchowski equation which describes the evolution of a tagged particle moving at constant speed in a random distribution of fixed particles. The volumes $v$ of the particles are independently distributed according to a probability distribution which decays asymptotically as a power law $v^{-σ}$. The validity of the eq…
▽ More
In this paper we consider the long time asymptotics of a linear version of the Smoluchowski equation which describes the evolution of a tagged particle moving at constant speed in a random distribution of fixed particles. The volumes $v$ of the particles are independently distributed according to a probability distribution which decays asymptotically as a power law $v^{-σ}$. The validity of the equation has been rigorously proved in \cite{NoV} for values of the exponent $σ>3$. The solutions of this equation display a rich structure of different asymptotic behaviours according to the different values of the exponent $σ$. Here we show that for $\frac{5}{3}<σ<2$ the linear Smoluchowski equation is well posed and that there exists a unique self-similar profile which is asymptotically stable.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.
-
Homogenization for the Poisson equation in randomly perforated domains under minimal assumptions on the size of the holes
Authors:
Arianna Giunti,
Richard Höfer,
Juan J. L. Velázquez
Abstract:
This paper deals with the homogenization of the Poisson equation in a bounded domain of $\mathbb{R}^d$, $d>2$, which is perforated by a random number of small spherical holes with random radii and positions. We show that for a class of stationary short-range correlated measures for the centres and radii of the holes, we recover in the homogenized limit an averaged analogue of the "strange term" ob…
▽ More
This paper deals with the homogenization of the Poisson equation in a bounded domain of $\mathbb{R}^d$, $d>2$, which is perforated by a random number of small spherical holes with random radii and positions. We show that for a class of stationary short-range correlated measures for the centres and radii of the holes, we recover in the homogenized limit an averaged analogue of the "strange term" obtained by Cioranescu and Murat in the periodic case [D. Cioranescu and F. Murat, \textit{Un term étrange venu d'ailleurs} (1986)]. We stress that we only require that the random radii have finite $(d-2)$-moment, which is the minimal assumption in order to ensure that the average of the capacity of the balls is finite. Under this assumption, there are holes which overlap with probability one. However, we show that homogenization occurs and that the clustering holes do not have any effect in the resulting homogenized equation.
△ Less
Submitted 27 March, 2018;
originally announced March 2018.
