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

Yves Bourgault

Yves Bourgault contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.AP

Trust snapshot

Quick read

Trust 19 - Baseline
5works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

5 published item(s)

preprint2020arXiv

Analysis of time-stepping methods for the monodomain model

To a large extent, the stiffness of the bidomain and monodomain models depends on the choice of the ionic model, which varies in terms of complexity and realism. In this paper, we compare and analyze a variety of time-stepping methods: explicit or semi-implicit, operator splitting, exponential, and deferred correction methods. We compare these methods for solving the bidomain model coupled with three ionic models of varying complexity and stiffness: the phenomenological Mitchell-Schaeffer model, the more realistic Beeler-Reuter model, and the stiff and very complex ten Tuscher-Noble-Noble-Panfilov (TNNP) model. For each method, we derive absolute stability criteria of the spatially discretized monodomain model and verify that the theoretical critical time-steps obtained closely match the ones in numerical experiments. We also verify that the numerical methods achieve an optimal order of convergence on the model variables and derived quantities (such as speed of the wave, depolarization time), and this in spite of the local non-differentiability of some of the ionic models. The efficiency of the different methods is also considered by comparing computational times for similar accuracy. Conclusions are drawn on the methods to be used to solve the monodomain model based on the model stiffness and complexity, measured respectively by the eigenvalues of the model's Jacobian and the number of variables, and based on strict stability and accuracy criteria.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient second-order semi-implicit finite element method for fourth-order nonlinear diffusion equations

We focus here on a class of fourth-order parabolic equations that can be written as a system of second-order equations by introducing an auxiliary variable. We design a novel second-order fully discrete mixed finite element method to approximate these equations. In our approach, we propose new techniques using the second-order backward differentiation formula for the time derivative and a special technique for the approximation of nonlinear terms. The use of the proposed technique for nonlinear terms makes the developed numerical scheme efficient in terms of computational cost since the proposed method only deals with a linear system at each time step and no iterative resolution is needed. A numerical convergence study is performed using the method of manufactured and analytical solutions of the system where we investigate different boundary conditions. With respect to the spatial discretization, convergence rates are found to at least match a priori error estimates available for linear problems. The convergence analysis is completed with an investigation of the temporal discretization where we numerically demonstrate the second-order time-accuracy of the proposed scheme using the method of reference solution. We present a series of numerical tests to demonstrate the efficiency and robustness of the proposed scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite difference and numerical differentiation: General formulae from deferred corrections

This paper provides a new approach to derive various arbitrary high order finite difference formulae for the numerical differentiation of analytic functions. In this approach, various first and second order formulae for the numerical approximation of analytic functions are given with error terms explicitly expanded as Taylor series of the analytic function. These lower order approximations are successively improved by one or two (two order improvement for centered formulae) to give finite difference formulae of arbitrary high order. The new approach allows to recover the standard backward, forward, and centered finite difference formulae which are given in terms of formal power series of finite difference operators. Examples of new formulae suited for deferred correction methods are given.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Anisotropic residual based a posteriori mesh adaptation in 2D: element based approach

An element based adaptation method is developed for an anisotropic a posteriori error estimator. The adaptation does not make use of a metric, but instead equidistributes the error over elements using local mesh modifications. Numerical results are reported, comparing with three popular anisotropic adaptation methods currently in use. It was found that the new method gives favourable results for controlling the energy norm of the error in terms of degrees of freedom at the cost of increased CPU usage. Additionally, we considered a new $L^2$ variant of the estimator. The estimator is shown to be conditionally equivalent to the exact $L^2$ error. We provide examples of adapted meshes with the $L^2$ estimator, and show that it gives greater control of the $L^2$ error compared with the original estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence rate for the method of moments with linear closure relations

We study linear closure relations for the moments' method applied to simple kinetic equations. The equations are linear collisional models (velocity jump processes) which are well suited to this type of approximation. In this simplified, 1 dimensional setting, we are able to prove stability estimates for the method (with a kinetic interpretation by a BGK model). Moreover we are also able to obtain convergence rates which automatically increase with the smoothness of the initial data.

0 citations0 reviewsNo code linkedOpen access