Researcher profile

Stig Larsson

Stig Larsson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

11 published item(s)

preprint2020arXiv

Mittag-Leffler Euler integrator for a stochastic fractional order equation with additive noise

Motivated by fractional derivative models in viscoelasticity, a class of semilinear stochastic Volterra integro-differential equations, and their deterministic counterparts, are considered. A generalized exponential Euler method, named here as the Mittag-Leffler Euler integrator, is used for the temporal discretization, while the spatial discretization is performed by the spectral Galerkin method. The temporal rate of strong convergence is found to be (almost) twice compared to when the backward Euler method is used together with a convolution quadrature for time discretization. Numerical experiments that validate the theory are presented.

preprint2019arXiv

A greedy algorithm for optimal heating in powder-bed-based additive manufacturing

Powder-bed-based additive manufacturing involves melting of a powder bed using a moving laser or electron beam as a heat source. In this paper, we formulate an optimization scheme that aims to control this type of melting. The goal consists of tracking maximum temperatures on lines that run along the beam path. Time-dependent beam parameters (more specifically, beam power, spot size, and speed) act as control functions. The scheme is greedy in the sense that it exploits local properties of the melt pool in order to divide a large optimization problem into several small ones. As illustrated by numerical examples, the scheme can resolve heat conduction issues such as concentrated heat accumulation at turning points and non-uniform melt depths.

preprint2019arXiv

Error estimates of the backward Euler-Maruyama method for multi-valued stochastic differential equations

In this paper, we derive error estimates of the backward Euler-Maruyama method applied to multi-valued stochastic differential equations. An important example of such an equation is a stochastic gradient flow whose associated potential is not continuously differentiable, but assumed to be convex. We show that the backward Euler-Maruyama method is well-defined and convergent of order at least $1/4$ with respect to the root-mean-square norm. Our error analysis relies on techniques for deterministic problems developed in [Nochetto, Savaré, and Verdi, Comm.\ Pure Appl.\ Math., 2000]. We verify that our setting applies to an overdamped Langevin equation with a discontinuous gradient and to a spatially semi-discrete approximation of the stochastic $p$-Laplace equation.

preprint2018arXiv

On a randomized backward Euler method for nonlinear evolution equations with time-irregular coefficients

In this paper we introduce a randomized version of the backward Euler method, that is applicable to stiff ordinary differential equations and nonlinear evolution equations with time-irregular coefficients. In the finite-dimensional case, we consider Carathéodory type functions satisfying a one-sided Lipschitz condition. After investigating the well-posedness and the stability properties of the randomized scheme, we prove the convergence to the exact solution with a rate of $0.5$ in the root-mean-square norm assuming only that the coefficient function is square integrable with respect to the temporal parameter. These results are then extended to the numerical solution of infinite-dimensional evolution equations under monotonicity and Lipschitz conditions. Here we consider a combination of the randomized backward Euler scheme with a Galerkin finite element method. We obtain error estimates that correspond to the regularity of the exact solution. The practicability of the randomized scheme is also illustrated through several numerical experiments.

preprint2015arXiv

A weak space-time formulation for the linear stochastic heat equation

We apply the well-known Banach-Necas-Babuska inf-sup theory in a stochastic setting to introduce a weak space-time formulation of the linear stochastic heat equation with additive noise. We give sufficient conditions on the the data and on the covariance operator associated to the driving Wiener process, in order to have existence and uniqueness of the solution. We show the relation of the obtained solution to the so-called mild solution and to the variational solution of the same problem. The spatial regularity of the solution is also discussed. Finally, an extension to the case of linear multiplicative noise is presented.

preprint2015arXiv

Compressive Space-Time Galerkin Discretizations of Parabolic Partial Differential Equations

We study linear parabolic initial-value problems in a space-time variational formulation based on fractional calculus. This formulation uses "time derivatives of order one half" on the bi-infinite time axis. We show that for linear, parabolic initial-boundary value problems on $(0,\infty)$, the corresponding bilinear form admits an inf-sup condition with sparse tensor product trial and test function spaces. We deduce optimality of compressive, space-time Galerkin discretizations, where stability of Galerkin approximations is implied by the well-posedness of the parabolic operator equation. The variational setting adopted here admits more general Riesz bases than previous work; in particular, no stability in negative order Sobolev spaces on the spatial or temporal domains is required of the Riesz bases accommodated by the present formulation. The trial and test spaces are based on Sobolev spaces of equal order $1/2$ with respect to the temporal variable. Sparse tensor products of multi-level decompositions of the spatial and temporal spaces in Galerkin discretizations lead to large, non-symmetric linear systems of equations. We prove that their condition numbers are uniformly bounded with respect to the discretization level. In terms of the total number of degrees of freedom, the convergence orders equal, up to logarithmic terms, those of best $N$-term approximations of solutions of the corresponding elliptic problems.

preprint2015arXiv

Discontinuous Galerkin method for an integro-differential equation modeling dynamic fractional order viscoelasticity

An integro-differential equation, modeling dynamic fractional order viscoelasticity, with a Mittag-Leffler type convolution kernel is considered. A discontinuous Galerkin method, based on piecewise constant polynomials is formulated for temporal semidiscretization of the problem. Stability estimates of the discrete problem are proved, that are used to prove optimal order a priori error estimates. The theory is illustrated by a numerical example.

preprint2013arXiv

Covariance structure of parabolic stochastic partial differential equations

In this paper parabolic random partial differential equations and parabolic stochastic partial differential equations driven by a Wiener process are considered. A deterministic, tensorized evolution equation for the second moment and the covariance of the solutions of the parabolic stochastic partial differential equations is derived. Well-posedness of a space-time weak variational formulation of this tensorized equation is established.

preprint2013arXiv

Posterior Contraction Rates for the Bayesian Approach to Linear Ill-Posed Inverse Problems

We consider a Bayesian nonparametric approach to a family of linear inverse problems in a separable Hilbert space setting with Gaussian noise. We assume Gaussian priors, which are conjugate to the model, and present a method of identifying the posterior using its precision operator. Working with the unbounded precision operator enables us to use partial differential equations (PDE) methodology to obtain rates of contraction of the posterior distribution to a Dirac measure centered on the true solution. Our methods assume a relatively weak relation between the prior covariance, noise covariance and forward operator, allowing for a wide range of applications.

preprint2011arXiv

Optimal Regularity for Semilinear Stochastic Partial Differential Equations with Multiplicative Noise

This paper deals with the spatial and temporal regularity of the unique Hilbert space valued mild solution to a semilinear stochastic partial differential equation with nonlinear terms that satisfy global Lipschitz conditions. It is shown that the mild solution has the same optimal regularity properties as the stochastic convolution. The proof is elementary and makes use of existing results on the regularity of the solution, in particular, the Hölder continuity with a non-optimal exponent.

preprint2010arXiv

Optimal closing of a pair trade with a model containing jumps

A pair trade is a portfolio consisting of a long position in one asset and a short position in another, and it is a widely applied investment strategy in the financial industry. Recently, Ekström, Lindberg and Tysk studied the problem of optimally closing a pair trading strategy when the difference of the two assets is modelled by an Ornstein-Uhlenbeck process. In this paper we study the same problem, but the model is generalized to also include jumps. More precisely we assume that the above difference is an Ornstein-Uhlenbeck type process, driven by a Lévy process of finite activity. We prove a verification theorem and analyze a numerical method for the associated free boundary problem. We prove rigorous error estimates, which are used to draw some conclusions from numerical simulations.