BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Antoine Tambue

Antoine Tambue contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.PRBiological PhysicsComputational Engineering, Finance, and Sciencemath.APmath.FAnlin.PSq-fin.PRTissues and Organs

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2020arXiv

A fitted finite volume method for stochastic optimal control Problems

In this article, we provide a numerical method based on fitted finite volume method to approximate the Hamilton-Jacobi-Bellman (HJB) equation coming from stochastic optimal control problems. The computational challenge is due to the nature of the HJB equation, which may be a second-order degenerated partial differential equation coupled with optimization. In the work, we discretize the HJB equation using the fitted finite volume method and show that matrix resulting from spatial discretization is an M-matrix. The optimization problem is solved at every time step using iterative method. Numerical results are presented to show the robustness of the fitted finite volume numerical method comparing to the standard finite difference method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal error estimate of the finite element approximation of second order semilinear non-autonomous parabolic PDEs

In this work, we investigate the numerical approximation of the second order non-autonomous semilnear parabolic partial differential equation (PDE) using the finite element method. To the best of our knowledge, only the linear case is investigated in the literature. Using an approach based on evolution operator depending on two parameters, we obtain the error estimate of the scheme toward the mild solution of the PDE under polynomial growth condition of the nonlinearity. Our convergence rate are obtain for smooth and non-smooth initial data and is similar to that of the autonomous case. Our convergence result for smooth initial data is very important in numerical analysis. For instance, it is one step forward in approximating non-autonomous stochastic partial differential equations by the finite element method. In addition, we provide realistic conditions on the nonlinearity, appropriated to achieve optimal convergence rate without logarithmic reduction by exploiting the smooth properties of the two parameters evolution operator.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A fitted L-Multi-point Flux Approximation method for pricing options

In this paper, we introduce a special kind of finite volume method called Multi-Point Flux Approximation method (MPFA) to price European and American options in two dimensional domain. We focus on the L-MPFA method for space discretization of the diffusion term of Black-Scholes operator. The degeneracy of the Black Scholes operator is tackled using the fitted finite volume method. This combination of fitted finite volume method and L-MPFA method coupled to upwind methods gives us a novel scheme called the fitted L-MPFA method. Numerical experiments show the accuracy of the novel fitted L-MPFA method comparing to well known schemes for pricing options.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Convergence of the Two Point Flux Approximation and a novel fitted Two-Point Flux Approximation method for pricing options

In this paper, we deal with numerical approximations for solving the Black-Scholes Partial Differential Equation (PDE). This PDE is well known to be degenerated. The space discretization is performed using the classical finite volume method with Two-Point Flux Approximation (TPFA) and a novel scheme called fitted Two-Point Flux Approximation (TPFA).The fitted Two-Point Flux Approximation (TPFA) combines the fitted finite volume method and the standard TPFA method. More precisely the fitted finite volume method is used when the stock price approaches zero in order to handle the degeneracy of the PDE and the TPFA method is used on the rest of space domain. This combination yields our novel fitted TPFA scheme. The Euler method is used for the time discretization. We provide the rigorous convergence proofs of the two fully discretized schemes. Numerical experiments to support our theoretical results are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimal strong convergence rates of some Euler-type timestepping schemes for the finite element discretization SPDEs driven by additive fractional Brownian motion and Poisson random measure

In this paper, we study the numerical approximation of a general second order semilinear stochastic partial differential equation (SPDE) driven by a additive fractional Brownian motion (fBm) with Hurst parameter $H>\frac 12$ and Poisson random measure, more realistic in modelling real world phenomena. To the best of our knowledge, numerical schemes for such SPDE have been lacked in scientific literature. The approximation is done with the standard finite element method in space and three Euler-type timestepping methods in time, more precisely linear implicit method, exponential integrator and exponential Rosenbrock scheme are used for time discretisation. In contract to the current literature in the field for SPDE driven only by fBm, our linear operator is not necessary self-adjoint and optimal strong convergence rates have been achieved for SPDE driven only by fBm and SPDE driven by fBm and Poisson measure. The results examine how the convergence orders depend on the regularity of the noise and the initial data and reveal that the full discretization attains an optimal convergence rate of order $\mathcal{O}(h^2+Δt)$ for the exponential integrator and implicit schemes (linear operator $A$ self-adjoint for implicit). Numerical experiments are provided to illustrate our theoretical results for the case of SPDE driven with fBm noise.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Strong convergence of the backward Euler approximation for the finite element discretization of semilinear parabolic SPDEs with non-global Lipschitz drift driven by additive noise

This paper deals with the backward Euler method applied to semilinear parabolic stochastic partial differential equations (SPDEs) driven by additive noise. The SPDE is discretized in space by the finite element method and in time by the backward Euler. We consider a larger class of nonlinear drift functions, which are of Nemytskii type and polynomial of any odd degree with negative leading term, instead of only dealing with the special case of stochastic Allen-Chan equation as in the up to date literature.Moreover our linear operator is of second order and not necessarily self-adjoint, therefore makes estimates more challenging than in the case of self-adjoint operator. We prove the strong convergence of our fully discrete schemes toward the mild solution and results indicate how the convergence rates depend on the regularities of the initial data and the noise. In particular, for trace class noise, we achieve convergence order $\mathcal{O} (h^{2}+Δt^{1-ε})$, where $ε>0$ is positive number, small enough. We also provide numerical experiments to illustrate our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Strong Convergence of the Linear Implicit Euler Method for the Finite Element Discretization of Semilinear SPDEs Driven by Multiplicative and Additive Noise

This paper aims to investigate the numerical approximation of a general second order parabolic stochastic partial differential equation(SPDE) driven by multiplicative and additive noise under more relaxed conditions. The SPDE is discretized in space by the finite element method and in time by the linear implicit Euler method. This extends the current results in the literature to not necessary self-adjoint operator with more general boundary conditions. As a consequence key part of the proof does not rely on the spectral decomposition of the linear operator. We achieve optimal convergence orders which depend on the regularity of the noise and the initial data. In particular, for multiplicative noise we achieve optimal order $\mathcal{O}(h^2+Δt^{1/2})$ and for additive noise, we achieve optimal order $\mathcal{O}(h^2+Δt)$. In contrast to current work in the literature, where the optimal convergence orders are achieved for additive noise by incorporating further regularity assumptions on the nonlinear drift function, our optimal convergence orders are obtained under only the standard Lipschitz condition of the nonlinear drift term. Numerical experiments to sustain our theoretical results are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A modified semi--implict Euler-Maruyama Scheme for finite element discretization of SPDEs with additive noise

We consider the numerical approximation of a general second order semi--linear parabolic stochastic partial differential equation (SPDE) driven by additive space-time noise. We introduce a new modified scheme using a linear functional of the noise with a semi--implicit Euler--Maruyama method in time and in space we analyse a finite element method (although extension to finite differences or finite volumes would be possible). We prove convergence in the root mean square $L^{2}$ norm for a diffusion reaction equation and diffusion advection reaction equation. We present numerical results for a linear reaction diffusion equation in two dimensions as well as a nonlinear example of two-dimensional stochastic advection diffusion reaction equation. We see from both the analysis and numerics that the proposed scheme has better convergence properties than the standard semi--implicit Euler--Maruyama method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Localized modulated wave solutions in diffusive glucose-insulin systems

We investigate intercellular insulin dynamics in an array of diffusively coupled pancreatic islet \b{eta}-cells. The cells are connected via gap junction coupling, where nearest neighbor interactions are included. Through the multiple scale expansion in the semi-discrete approximation, we show that the insulin dynamics can be governed by the complex Ginzburg-Landau equation. The localized solutions of this equation are reported. The results suggest from the biophysical point of view that the insulin propagates in pancreatic islet \b{eta}-cells using both temporal and spatial dimensions in the form of localized modulated waves.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Localized numerical impulses solutions in diffuse neural networks modeled by the complex fractional Ginzburg-Landau equation

We investigate localized wave solutions in a network of Hindmarsh-Rose neural model taking into account the long-range diffusive couplings. We show by a specific analytical technique that the model equations in the infrared limit (wave number $k\rightarrow 0$) can be governed by the complex fractional Ginzburg-Landau (CFGL) equation. According to the stiffness of the system, we propose both the semi and the linearly implicit Riesz fractional finite-difference schemes to solve efficiently the CFGL equation. The obtained fractional numerical solutions for the nerve impulse reveal localized short impulse properties. We also show the equivalence between the continuous CFGL and the discrete Hindmarsh-Rose models for relatively large network.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An exponential Integrator for finite volume discretization of nonlinear parabolic differential equation

We consider the numerical approximation of a general second order semi--linear parabolic partial differential equation. Equations of this type arise in many contexts, such as transport in porous media which is fundamental in many geo-engineering applications, including oil and gas recovery from subsurface. Using the finite volume with two-point flux approximation on regular mesh combined with exponential time differencing of order one (ETD1) for temporal discretization, we derive the $L^{2}$ estimate under the assumption that the non linear term is locally Lipschitz. Numerical simulations to sustain the theoretical results are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of the semi-tamed and tamed Euler schemes for stochastic differential equations with jumps under non-global Lipschitz condition

Under non-global Lipschitz condition, Euler Explicit method fails to converge strongly to the exact solution, while Euler implicit method converges but requires much computational efforts. Tamed scheme was first introduced in [2] to overcome this failure of the standard explicit method. This technique is extended to SDEs driven by Poisson jump in [3] where several schemes were analyzed. In this work, we investigate their nonlinear stability under non-global Lipschitz and their linear stability. Numerical simulations to sustain the theoretical results are also provided.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weak convergence for a stochastic exponential integrator and finite element discretization of SPDE for multiplicative \& additive noise

We consider a finite element approximation of a general semi-linear stochastic partial differential equation (SPDE) driven by space-time multiplicative and additive noise. We examine the full weak convergence rate of the exponential Euler scheme when the linear operator is self adjoint and provide preliminaries results toward the full weak convergence rate for non-self-adjoint linear operator. Key part of the proof does not rely on Malliavin calculus. Depending of the regularity of the noise and the initial solution, we found that in some cases the rate of weak convergence is twice the rate of the strong convergence. Our convergence rate is in agreement with some numerical results in two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Stochastic Delay Model for Pricing Debt and Equity: Numerical Techniques and Applications

In the accompanied paper [14], a delayed nonlinear model for pricing corporate liabilities was developed. Using self-financed strategy and duplication we were able to derive two Random Partial Differential Equations (RPDEs) describing the evolution of debt and equity values of the corporate in the last delay period interval. In this paper, we provide numerical techniques to solve our delayed nonlinear model along with the corresponding RPDEs modeling the debt and equity values of the corporate. Using financial data from some firms, we compare numerical solutions from both our nonlinear model and classical Merton model [7] to the real corporate data. From this comparison, it comes up that in corporate finance the past dependence of the firm value process may be an important feature and therefore should not be ignored.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Stochastic Delay Model For Pricing Debt And Loan Guarantees: Theoretical Results

We consider that the price of a firm follows a non linear stochastic delay differential equation. We also assume that any claim value whose value depends on firm value and time follows a non linear stochastic delay differential equation. Using self-financed strategy and replication we are able to derive a Random Partial Differential Equation (RPDE) satisfied by any corporate claim whose value is a function of firm value and time. Under specific final and boundary conditions, we solve the RPDE for the debt value and loan guarantees within a single period and homogeneous class of debt.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient Simulation of Geothermal Processes in Heterogeneous Porous Media based on the Exponential Rosenbrock-Euler and Rosenbrock-type Methods

Simulation of geothermal systems is challenging due to coupled physical processes in highly heterogeneous media. Combining the exponential Rosenbrock--Euler and Rosenbrock-type methods with control-volume (two-point flux approximation) space discretizations leads to efficient numerical techniques for simulating geothermal systems. In terms of efficiency and accuracy, the exponential Rosenbrock--Euler time integrator has advantages over standard time-dicretization schemes, which suffer from time-step restrictions or excessive numerical diffusion when advection processes are dominating. Based on linearization of the equation at each time step, we make use of matrix exponentials of the Jacobian from the spatial discretization, which provide the exact solution in time for the linearized equations. This is at the expense of computing the matrix exponentials of the stiff Jacobian matrix, together with propagating a linearized system. However, using a Krylov subspace or Leja points techniques make these computations efficient. The Rosenbrock-type methods use the appropriate rational functions of the Jacobian from the spatial discretization. The parameters in these schemes are found in consistency with the required order of convergence in time. As a result, these schemes are A-stable and only a few linear systems are solved at each time step. The efficiency of the methods compared to standard time-discretization techniques are demonstrated in numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stochastic exponential integrators for finite element discretization of SPDEs for multiplicative and additive noise

We consider the numerical approximation of a general second order semi--linear parabolic stochastic partial differential equation (SPDEs) driven by space-time noise, for multiplicative and additive noise. We examine convergence of exponential integrators for multiplicative and additive noise. We consider noise that is in trace class and give a convergence proof in the mean square $L^{2}$ norm. We discretize in space with the finite element method and in our implementation we examine both the finite element and the finite volume methods. We present results for a linear reaction diffusion equation in two dimensions as well as a nonlinear example of two-dimensional stochastic advection diffusion reaction equation motivated from realistic porous media flow.

0 citations0 reviewsNo code linkedOpen access