-
Approximate solutions to a nonlinear functional differential equation
Authors:
Nicholas Hale,
Enrique Thomann,
JAC Weideman
Abstract:
A functional differential equation related to the logistic equation is studied by a combination of numerical and perturbation methods. Parameter regions are identified where the solution to the nonlinear problem is approximated well by known series solutions of the linear version of the equation. The solution space for a certain class of functions is then mapped out using a continuation approach.
A functional differential equation related to the logistic equation is studied by a combination of numerical and perturbation methods. Parameter regions are identified where the solution to the nonlinear problem is approximated well by known series solutions of the linear version of the equation. The solution space for a certain class of functions is then mapped out using a continuation approach.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
On the impact of insurance on households susceptible to random proportional losses: An analysis of poverty trapping
Authors:
Kira Henshaw,
Jorge M. Ramirez,
José M. Flores-Contró,
Enrique A. Thomann,
Sooie-Hoe Loke,
Corina Constantinescu
Abstract:
In this paper, we consider a risk process with deterministic growth and multiplicative jumps to model the capital of a low-income household. Reflecting the high-risk nature of the low-income environment, capital losses are assumed to be proportional to the level of accumulated capital at the jump time. Our aim is to derive the probability that a household falls below the poverty line, i.e. the tra…
▽ More
In this paper, we consider a risk process with deterministic growth and multiplicative jumps to model the capital of a low-income household. Reflecting the high-risk nature of the low-income environment, capital losses are assumed to be proportional to the level of accumulated capital at the jump time. Our aim is to derive the probability that a household falls below the poverty line, i.e. the trapping probability, where ``trapping" occurs when the level of capital of a household holds falls below the poverty line, to an area from which it is difficult to escape without external help. Considering the remaining proportion of capital to be distributed as a special case of the beta distribution, closed-form expressions for the trapping probability are obtained via analysis of the Laplace transform of the infinitesimal generator of the process. To study the impact of insurance on this probability, introduction of an insurance product offering proportional coverage is presented. The infinitesimal generator of the insured process gives rise to non-local differential equations. To overcome this, we propose a recursive method for deriving a closed-form solution of the integro-differential equation associated with the infinitesimal generator of the insured process and provide a numerical estimation method for obtaining the trapping probability. Constraints on the rate parameters of the process that prevent certain trapping are derived in both the uninsured and insured cases using classical results from risk theory.
△ Less
Submitted 22 September, 2023;
originally announced October 2023.
-
Errata to Stochastic explosion and non-uniqueness for $α$-Riccati equation
Authors:
Radu Dascaliuc,
Tuan N. Pham,
Enrique A. Thomann,
Edward C. Waymire
Abstract:
An error occurs in a part of the statement and proof of Proposition 2.2 in Jour. Math. Anal. and Appl., 476, (2019), 53-85 that is corrected in this erratum. The revised result reveals a new and unexpected critical phenomenon, having further implications for non-uniqueness of solutions to a nonlinear differential equation of the Riccati type.
An error occurs in a part of the statement and proof of Proposition 2.2 in Jour. Math. Anal. and Appl., 476, (2019), 53-85 that is corrected in this erratum. The revised result reveals a new and unexpected critical phenomenon, having further implications for non-uniqueness of solutions to a nonlinear differential equation of the Riccati type.
△ Less
Submitted 5 August, 2023; v1 submitted 2 March, 2023;
originally announced March 2023.
-
Doubly stochastic Yule cascades (Part II): The explosion problem in the non-reversible case
Authors:
Radu Dascaliuc,
Tuan N. Pham,
Enrique Thomann,
Edward C. Waymire
Abstract:
We analyze the explosion problem for a class of stochastic models introduced in Part I (arXiv:2103.06912), referred to as doubly stochastic Yule cascades. These models arise naturally in the construction of solutions to evolutionary PDEs as well as in purely probabilistic first passage percolation phenomena having a Markov-type statistical dependence, new for this context. Using cut-set arguments…
▽ More
We analyze the explosion problem for a class of stochastic models introduced in Part I (arXiv:2103.06912), referred to as doubly stochastic Yule cascades. These models arise naturally in the construction of solutions to evolutionary PDEs as well as in purely probabilistic first passage percolation phenomena having a Markov-type statistical dependence, new for this context. Using cut-set arguments and a greedy algorithm, we respectively establish criteria for non-explosion and explosion without requiring the time-reversibility of the underlying branching Markov chain (a condition required in Part I). Notable applications include the explosion of the self-similar cascade of the Navier-Stokes equations in dimension $d=3$ and non-explosion in dimensions $d\ge 12$.
△ Less
Submitted 2 December, 2021; v1 submitted 28 July, 2021;
originally announced July 2021.
-
Doubly Stochastic Yule Cascades (Part I): The explosion problem in the time-reversible case
Authors:
Radu Dascaliuc,
Tuan N. Pham,
Enrique Thomann,
Edward C. Waymire
Abstract:
Motivated by the probabilistic methods for nonlinear differential equations introduced by McKean (1975) for the Kolmogorov-Petrovski-Piskunov (KPP) equation, and by Le Jan and Sznitman (1997) for the incompressible Navier-Stokes equations, we identify a new class of stochastic cascade models, referred to as Doubly Stochastic Yule cascades. We establish non-explosion criteria under the assumption t…
▽ More
Motivated by the probabilistic methods for nonlinear differential equations introduced by McKean (1975) for the Kolmogorov-Petrovski-Piskunov (KPP) equation, and by Le Jan and Sznitman (1997) for the incompressible Navier-Stokes equations, we identify a new class of stochastic cascade models, referred to as Doubly Stochastic Yule cascades. We establish non-explosion criteria under the assumption that the randomization of Yule intensities from generation to generation is by an ergodic time-reversible Markov process. In addition to the cascade models that arise in the analysis of certain deterministic nonlinear differential equations, this model includes the multiplicative branching random walks, the branching Markov processes, and the stochastic generalizations of the percolation and/or cell aging models introduced by Aldous and Shields (1988) and independently by Athreya (1985).
△ Less
Submitted 30 September, 2022; v1 submitted 11 March, 2021;
originally announced March 2021.
-
On Le Jan-Sznitman's stochastic approach to the Navier-Stokes equations
Authors:
Radu Dascaliuc,
Tuan N. Pham,
Enrique Thomann
Abstract:
The paper explores the symbiotic relation between the Navier-Stokes equations and the associated stochastic cascades. Specifically, we examine how some well-known existence and uniqueness results for the Navier-Stokes equations can inform about the probabilistic features of the associated stochastic cascades, and how some probabilistic features of the stochastic cascades can, in turn, inform about…
▽ More
The paper explores the symbiotic relation between the Navier-Stokes equations and the associated stochastic cascades. Specifically, we examine how some well-known existence and uniqueness results for the Navier-Stokes equations can inform about the probabilistic features of the associated stochastic cascades, and how some probabilistic features of the stochastic cascades can, in turn, inform about the existence and uniqueness (or the lack thereof) of solutions. Our method of incorporating the stochastic explosion gives a simpler and more natural method to construct the solution compared to the original construction by Le Jan and Sznitman. This new stochastic construction is then used to show the finite-time blowup and non-uniqueness of the initial value problem for the Montgomery-Smith equation. We exploit symmetry properties inherent in our construction to give a simple proof of the global well-posedness results for small initial data in scale-critical Fourier-Besov spaces. We also obtain the pointwise convergence of the Picard's iteration associated with the Fourier-transformed Navier-Stokes equations.
△ Less
Submitted 7 December, 2021; v1 submitted 12 October, 2019;
originally announced October 2019.
-
Global analysis of a predator-prey model with variable predator search rate
Authors:
Ben Dalziel,
Enrique Thomann,
Jan Medlock,
Patrick De Leenheer
Abstract:
We consider a modified Rosenzweig-MacArthur predator-prey model, based on the premise that the search rate of predators is dependent on the prey density, rather than constant. A complete analysis of the global behavior of the model is presented, and shows that the model exhibits a dichotomy similar to the classical Rosenzweig-MacArthur model: either the coexistence steady state is globally stable;…
▽ More
We consider a modified Rosenzweig-MacArthur predator-prey model, based on the premise that the search rate of predators is dependent on the prey density, rather than constant. A complete analysis of the global behavior of the model is presented, and shows that the model exhibits a dichotomy similar to the classical Rosenzweig-MacArthur model: either the coexistence steady state is globally stable; or it is unstable, and then a unique, globally stable limit cycle exists. We discuss the similarities, but also important differences between our model and the Rosenweig-MacArthur model. The main differences are that: 1. The paradox of enrichment which always occurs in the Rosenzweig-MacArthur model, does not always occur here, and 2. Even when the paradox of enrichment occurs, predators can adapt by lowering their search rate, and effectively stabilize the system.
△ Less
Submitted 13 June, 2019;
originally announced June 2019.
-
A Delayed Yule Process
Authors:
Radu Dascaliuc,
Nicholas Michalowski,
Enrique Thomann,
Edward C. Waymire
Abstract:
In now classic work, David Kendall (1966) recognized that the Yule process and Poisson process could be related by a (random) time change. Furthermore, he showed that the Yule population size rescaled by its mean has an almost sure exponentially distributed limit as $t\to \infty$. In this note we introduce a class of coupled delayed Yule processes parameterized by $0 < α\le 1$ that includes the Po…
▽ More
In now classic work, David Kendall (1966) recognized that the Yule process and Poisson process could be related by a (random) time change. Furthermore, he showed that the Yule population size rescaled by its mean has an almost sure exponentially distributed limit as $t\to \infty$. In this note we introduce a class of coupled delayed Yule processes parameterized by $0 < α\le 1$ that includes the Poisson process at $α= {1/2}$. Moreover we extend Kendall's limit theorem to include a larger class of positive martingales derived from functionals that gauge the population genealogy. A somewhat surprising connection with the Holley-Liggett smoothing transformation also emerges in this context. Specifically, the latter is exploited to uniquely characterize the moment generating functions of distributions of the limit martingales, generalizing Kendall's mean one exponential limit.
△ Less
Submitted 7 July, 2016;
originally announced July 2016.
-
Continuity of Local Time: An applied perspective
Authors:
Jorge M. Ramirez,
Edward C. Waymire,
Enrique A. Thomann
Abstract:
Continuity of local time for Brownian motion ranks among the most notable mathematical results in the theory of stochastic processes. This article addresses its implications from the point of view of applications. In particular an extension of previous results on an explicit role of continuity of (natural) local time is obtained for applications to recent classes of problems in physics, biology an…
▽ More
Continuity of local time for Brownian motion ranks among the most notable mathematical results in the theory of stochastic processes. This article addresses its implications from the point of view of applications. In particular an extension of previous results on an explicit role of continuity of (natural) local time is obtained for applications to recent classes of problems in physics, biology and finance involving discontinuities in a dispersion coefficient. The main theorem and its corollary provide physical principles that relate macro scale continuity of deterministic quantities to micro scale continuity of the (stochastic) local time.
△ Less
Submitted 16 March, 2015;
originally announced March 2015.
-
Symmetry Breaking and Uniqueness for the Incompressible Navier-Stokes Equations
Authors:
Radu Dascaliuc,
Nicholas Michalowski,
Enrique Thomann,
Edward C. Waymire
Abstract:
The present article establishes connections between the structure of the deterministic Navier-Stokes equations and the structure of (similarity) equations that govern self-similar solutions as expected values of certain naturally associated stochastic cascades. A principle result is that explosion criteria for the stochastic cascades involved in the probabilistic representations of solutions to th…
▽ More
The present article establishes connections between the structure of the deterministic Navier-Stokes equations and the structure of (similarity) equations that govern self-similar solutions as expected values of certain naturally associated stochastic cascades. A principle result is that explosion criteria for the stochastic cascades involved in the probabilistic representations of solutions to the respective equations coincide. While the uniqueness problem itself remains unresolved, these connections provide interesting problems and possible methods for investigating symmetry breaking and the uniqueness problem for Navier-Stokes equations. In particular, new branching Markov chains, including a dilogarithmic branching random walk on the multiplicative group $(0,\infty)$, naturally arise as a result of this investigation.
△ Less
Submitted 24 February, 2015;
originally announced February 2015.
-
Ergodic and Mixing Properties of the Boussinesq Equations with a Degenerate Random Forcing
Authors:
Juraj Földes,
Nathan Glatt-Holtz,
Geordie Richards,
Enrique Thomann
Abstract:
We establish the existence, uniqueness and attraction properties of an ergodic invariant measure for the Boussinesq Equations in the presence of a degenerate stochastic forcing acting only in the temperature equation and only at the largest spatial scales. The central challenge is to establish time asymptotic smoothing properties of the Markovian dynamics corresponding to this system. Towards this…
▽ More
We establish the existence, uniqueness and attraction properties of an ergodic invariant measure for the Boussinesq Equations in the presence of a degenerate stochastic forcing acting only in the temperature equation and only at the largest spatial scales. The central challenge is to establish time asymptotic smoothing properties of the Markovian dynamics corresponding to this system. Towards this aim we encounter a Lie bracket structure in the associated vector fields with a complicated dependence on solutions. This leads us to develop a novel Hörmander-type condition for infinite-dimensional systems. Demonstrating the sufficiency of this condition requires new techniques for the spectral analysis of the Malliavin covariance matrix.
△ Less
Submitted 14 November, 2013;
originally announced November 2013.
-
Numerical Methods for Linear Diffusion Equations in the Presence of an Interface
Authors:
V. A. Bokil,
N. L. Gibson,
S. L. Nguyen,
E. A. Thomann,
E. Waymire
Abstract:
We consider numerical methods for linear parabolic equations in one spatial dimension having piecewise constant diffusion coefficients defined by a one parameter family of interface conditions at the discontinuity. We construct immersed interface finite element methods for an alternative formulation of the original deterministic diffusion problem in which the interface condition is recast as a nat…
▽ More
We consider numerical methods for linear parabolic equations in one spatial dimension having piecewise constant diffusion coefficients defined by a one parameter family of interface conditions at the discontinuity. We construct immersed interface finite element methods for an alternative formulation of the original deterministic diffusion problem in which the interface condition is recast as a natural condition on the interfacial flux for which the given operator is self adjoint. An Euler-Maruyama method is developed for the stochastic differential equation corresponding to the alternative divergence formulation of the equation having a discontinuous coefficient and a one-parameter family of interface conditions. We then prove convergence estimates for the Euler scheme. The main goal is to develop numerical schemes that can accommodate specification of any one of the possible interface conditions, and to illustrate the implementation for each of the deterministic and stochastic formulations, respectively. The issues pertaining to speed-ups of the numerical schemes are left to future work.
△ Less
Submitted 30 October, 2013;
originally announced October 2013.
-
Advection-Dispersion Across Interfaces
Authors:
Jorge M. Ramirez,
Enrique A. Thomann,
Edward C. Waymire
Abstract:
This article concerns a systemic manifestation of small scale interfacial heterogeneities in large scale quantities of interest to a variety of diverse applications spanning the earth, biological and ecological sciences. Beginning with formulations in terms of partial differential equations governing the conservative, advective-dispersive transport of mass concentrations in divergence form, the sp…
▽ More
This article concerns a systemic manifestation of small scale interfacial heterogeneities in large scale quantities of interest to a variety of diverse applications spanning the earth, biological and ecological sciences. Beginning with formulations in terms of partial differential equations governing the conservative, advective-dispersive transport of mass concentrations in divergence form, the specific interfacial heterogeneities are introduced in terms of (spatial) discontinuities in the diffusion coefficient across a lower-dimensional hypersurface. A pathway to an equivalent stochastic formulation is then developed with special attention to the interfacial effects in various functionals such as first passage times, occupation times and local times. That an appreciable theory is achievable within a framework of applications involving one-dimensional models having piecewise constant coefficients greatly facilitates our goal of a gentle introduction to some rather dramatic mathematical consequences of interfacial effects that can be used to predict structure and to inform modeling.
△ Less
Submitted 23 December, 2013; v1 submitted 28 October, 2013;
originally announced October 2013.
-
Interfacial Phenomena and Natural Local Time
Authors:
Thilanka Appuhamillage,
Vrushali Bokil,
Enrique Thomann,
Edward Waymire,
Brian Wood
Abstract:
This article addresses a modification of local time for stochastic processes, to be referred to as `natural local time'. It is prompted by theoretical developments arising in mathematical treatments of recent experiments and observations of phenomena in the geophysical and biological sciences pertaining to dispersion in the presence of an interface of discontinuity in dispersion coefficients. The…
▽ More
This article addresses a modification of local time for stochastic processes, to be referred to as `natural local time'. It is prompted by theoretical developments arising in mathematical treatments of recent experiments and observations of phenomena in the geophysical and biological sciences pertaining to dispersion in the presence of an interface of discontinuity in dispersion coefficients. The results illustrate new ways in which to use the theory of stochastic processes to infer macro scale parameters and behavior from micro scale observations in particular heterogeneous environments.
△ Less
Submitted 1 April, 2012;
originally announced April 2012.
-
First Passage Times and Breakthrough Curves Associated with Interfacial Phenomena
Authors:
Edward C. Waymire,
Thilanka A. Appuhamillage,
Vrushali A. Bokil,
Enrique A. Thomann,
Brian D. Wood
Abstract:
Advection and dispersion in highly heterogeneous environments involving interfacial discontinuities in the corresponding drift and dispersion rates are described through disparate examples from the physical and biological sciences. A mathematical framework is formulated to address specific empirical phenomena involving first passage time and occupation time functionals observed in relation to the…
▽ More
Advection and dispersion in highly heterogeneous environments involving interfacial discontinuities in the corresponding drift and dispersion rates are described through disparate examples from the physical and biological sciences. A mathematical framework is formulated to address specific empirical phenomena involving first passage time and occupation time functionals observed in relation to the interfacial parameters.
△ Less
Submitted 21 June, 2011;
originally announced June 2011.
-
Corrections for "Occupation and local times for skew Brownian motion with applications to dispersion across an interface"
Authors:
Thilanka Appuhamillage,
Vrushali Bokil,
Enrique Thomann,
Edward Waymire,
Brian Wood
Abstract:
We are making corrections and acknowledging colleagues that pointed out mistakes in our recent paper titled "Occupation and local times for skew Brownian motion with applications to dispersion across an interface" which was published in Annals of Applied Probability (2011) 21(1) 183-214. Specifically the corrections are: 1. The restriction of $γ$ to non negative values in Theorem 1.3 is not needed…
▽ More
We are making corrections and acknowledging colleagues that pointed out mistakes in our recent paper titled "Occupation and local times for skew Brownian motion with applications to dispersion across an interface" which was published in Annals of Applied Probability (2011) 21(1) 183-214. Specifically the corrections are: 1. The restriction of $γ$ to non negative values in Theorem 1.3 is not needed. But one has probabilistic interpretation only when $γ$ is non negative. 2. State the correct formulas in Corollary 3.3 as their were computational errors in the original formulas. We thank Pierre Etoiré and Miguel Martinez for pointing out these errors.
△ Less
Submitted 3 January, 2012; v1 submitted 27 September, 2010;
originally announced September 2010.