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

Assyr Abdulle

Assyr Abdulle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.APmath.PR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
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

8 published item(s)

preprint2022arXiv

Eigenfunction martingale estimating functions and filtered data for drift estimation of discretely observed multiscale diffusions

We propose a novel method for drift estimation of multiscale diffusion processes when a sequence of discrete observations is given. For the Langevin dynamics in a two-scale potential, our approach relies on the eigenvalues and the eigenfunctions of the homogenized dynamics. Our first estimator is derived from a martingale estimating function of the generator of the homogenized diffusion process. However, the unbiasedness of the estimator depends on the rate with which the observations are sampled. We therefore introduce a second estimator which relies also on filtering the data and we prove that it is asymptotically unbiased independently of the sampling rate. A series of numerical experiments illustrate the reliability and efficiency of our different estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Explicit stabilized multirate method for stiff differential equations

Stabilized Runge-Kutta methods are especially efficient for the numerical solution of large systems of stiff nonlinear differential equations because they are fully explicit. For semi-discrete parabolic problems, for instance, stabilized Runge-Kutta methods overcome the stringent stability condition of standard methods without sacrificing explicitness. However, when stiffness is only induced by a few components, as in the presence of spatially local mesh refinement, their efficiency deteriorates. To remove the crippling effect of a few severely stiff components on the entire system of differential equations, we derive a modified equation, whose stiffness solely depend on the remaining mildly stiff components. By applying stabilized Runge-Kutta methods to this modified equation, we then devise an explicit multirate Runge-Kutta-Chebyshev (mRKC) method whose stability conditions are independent of a few severely stiff components. Stability of the mRKC method is proved for a model problem, whereas its efficiency and usefulness are demonstrated through a series of numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A parabolic local problem with exponential decay of the resonance error for numerical homogenization

This paper aims at an accurate and efficient computation of effective quantities, e.g., the homogenized coefficients for approximating the solutions to partial differential equations with oscillatory coefficients. Typical multiscale methods are based on a micro-macro coupling, where the macro model describes the coarse scale behaviour, and the micro model is solved only locally to upscale the effective quantities, which are missing in the macro model. The fact that the micro problems are solved over small domains within the entire macroscopic domain, implies imposing artificial boundary conditions on the boundary of the microscopic domains. A naive treatment of these artificial boundary conditions leads to a first order error in $\varepsilon/δ$, where $\varepsilon < δ$ represents the characteristic length of the small scale oscillations and $δ^d$ is the size of micro domain. This error dominates all other errors originating from the discretization of the macro and the micro problems, and its reduction is a main issue in today&#39;s engineering multiscale computations. The objective of the present work is to analyze a parabolic approach, first announced in [A. Abdulle, D. Arjmand, E. Paganoni, C. R. Acad. Sci. Paris, Ser. I, 2019], for computing the homogenized coefficients with arbitrarily high convergence rates in $\varepsilon/δ$. The analysis covers the setting of periodic micro structure, and numerical simulations are provided to verify the theoretical findings for more general settings, e.g. random stationary micro structures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A posteriori error analysis of a local adaptive discontinuous Galerkin method for convection-diffusion-reaction equations

We introduce a local adaptive discontinuous Galerkin method for convection-diffusion-reaction equations. The proposed method is based on a coarse grid and iteratively improves the solution&#39;s accuracy by solving local elliptic problems in refined subdomains. For purely diffusion problems, we already proved that this scheme converges under minimal regularity assumptions [A. Abdulle and G.Rosilho de Souza, ESAIM: M2AN, 53(4):1269--1303, 2019]. In this paper, we provide an algorithm for the automatic identification of the local elliptic problems&#39; subdomains employing a flux reconstruction strategy. Reliable error estimators are derived for the local adaptive method. Numerical comparisons with a classical nonlocal adaptive algorithm illustrate the efficiency of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analytical and numerical study of a modified cell problem for the numerical homogenization of multiscale random fields

A central question in numerical homogenization of partial differential equations with multiscale coefficients is the accurate computation of effective quantities, such as the homogenized coefficients. Computing homogenized coefficients requires solving local corrector problems followed by upscaling relevant local data. The most naive way of computing homogenized coefficients is by solving a local elliptic problem, which is known to suffer from the so-called resonance error dominating all other errors inherent in multiscale computations. A far more efficient modelling strategy, based on adding an exponential correction term to the standard local elliptic problem, has recently been proved to result in exponentially decaying error bounds with respect to the size of the local geometry. The questions in relation with the accuracy and computational efficiency of this approach has been previously addressed in the context of periodic homogenization. The present article concerns the extension of mathematical and numerical study of this modified elliptic corrector problem to stochastic homogenization problems. In particular, we assume a stationary, ergodic micro-structure and i) establish the well-posedness of the corrector equation, ii) analyse the bias (or the systematic error) originating from additional exponential correction term in the model. Numerical results corroborating our theoretical findings are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective models and numerical homogenization for wave propagation in heterogeneous media on arbitrary timescales

A family of effective equations for wave propagation in periodic media for arbitrary timescales $\mathcal{O}(\varepsilon^{-α})$, where $\varepsilon\ll1$ is the period of the tensor describing the medium, is proposed. The well-posedness of the effective equations of the family is ensured without requiring a regularization process as in previous models [A. Benoit and A. Gloria, arXiv:1701.08600, 2017], [G. Allaire, A. Lamacz, and J. Rauch, arXiv:1803.09455, 2018]. The effective solutions in the family are proved to be $\varepsilon$ close to the original wave in a norm equivalent to the $L^\infty(0,\varepsilon^{-α}T;L^2(Ω))$ norm. In addition, a numerical procedure for the computation of the effective tensors of arbitrary order is provided. In particular, we present a new relation between the correctors of arbitrary order, which allows to substantially reduce the computational cost of the effective tensors of arbitrary order. This relation is not limited to the effective equations presented in this paper and can be used to compute the effective tensors of alternative effective models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Instabilities and order reduction phenomenon of an interpolation based multirate Runge-Kutta-Chebyshev method

An explicit stabilized additive Runge-Kutta scheme is proposed. The method is based on a splitting of the problem in severely stiff and mildly stiff subproblems, which are then independently solved using a Runge-Kutta-Chebyshev scheme. The number of stages is adapted according to the subproblem&#39;s stiffness and leads to asynchronous integration needing ghost values. Whenever ghost values are needed, linear interpolation in time between stages is employed. One important application of the scheme is for parabolic partial differential equations discretized on a nonuniform grid. The goal of this paper is to introduce the scheme and prove on a model problem that linear interpolations trigger instabilities into the method. Furthermore, we show that it suffers from an order reduction phenomenon. The theoretical results are confirmed numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Random time step probabilistic methods for uncertainty quantification in chaotic and geometric numerical integration

A novel probabilistic numerical method for quantifying the uncertainty induced by the time integration of ordinary differential equations (ODEs) is introduced. Departing from the classical strategy to randomize ODE solvers by adding a random forcing term, we show that a probability measure over the numerical solution of ODEs can be obtained by introducing suitable random time-steps in a classical time integrator. This intrinsic randomization allows for the conservation of geometric properties of the underlying deterministic integrator such as mass conservation, symplecticity or conservation of first integrals. Weak and mean-square convergence analysis are derived. We also analyse the convergence of the Monte Carlo estimator for the proposed random time step method and show that the measure obtained with repeated sampling converges in the mean-square sense independently of the number of samples. Numerical examples including chaotic Hamiltonian systems, chemical reactions and Bayesian inferential problems illustrate the accuracy, robustness and versatility of our probabilistic numerical method.

0 citations0 reviewsNo code linkedOpen access