liu.seSearch for publications in DiVA
Change search
Refine search result
12345 1 - 50 of 243
CiteExportLink to result list
Permanent link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • oxford
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Rows per page
  • 5
  • 10
  • 20
  • 50
  • 100
  • 250
Sort
  • Standard (Relevance)
  • Author A-Ö
  • Author Ö-A
  • Title A-Ö
  • Title Ö-A
  • Publication type A-Ö
  • Publication type Ö-A
  • Issued (Oldest first)
  • Issued (Newest first)
  • Created (Oldest first)
  • Created (Newest first)
  • Last updated (Oldest first)
  • Last updated (Newest first)
  • Disputation date (earliest first)
  • Disputation date (latest first)
  • Standard (Relevance)
  • Author A-Ö
  • Author Ö-A
  • Title A-Ö
  • Title Ö-A
  • Publication type A-Ö
  • Publication type Ö-A
  • Issued (Oldest first)
  • Issued (Newest first)
  • Created (Oldest first)
  • Created (Newest first)
  • Last updated (Oldest first)
  • Last updated (Newest first)
  • Disputation date (earliest first)
  • Disputation date (latest first)
Select
The maximal number of hits you can export is 250. When you want to export more records please use the Create feeds function.
  • 1.
    Abgrall, Remi
    et al.
    Institute of Mathematics, University of Zurich, Switzerland.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering. Department of Mathematics and Applied Mathematics, University of Johannesburg, South Africa.
    Öffner, Philipp
    Institute of Mathematics, University of Zurich, Switzerland. Institute of Mathematics, Johannes Gutenberg-Universtiy, Germany.
    Tokareva, Svetlana
    Theoretical Division, Applied Mathematics and Plasma Physics Group (T-5), Los Alamos National Laboratory, USA.
    Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems2020In: Journal of Scientific Computing, ISSN 0885-7474, E-ISSN 1573-7691, Vol. 85, no 2, article id 43Article in journal (Refereed)
    Abstract [en]

    In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

    Download full text (pdf)
    fulltext
  • 2.
    Achieng, Pauline
    et al.
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Berntsson, Fredrik
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Chepkorir, Jennifer
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Analysis of Dirichlet–Robin Iterations for Solving the Cauchy Problem for Elliptic Equations2021In: Bulletin of the Iranian Mathematical Society, ISSN 1735-8515, Vol. 47, p. 1681-1699Article in journal (Refereed)
    Abstract [en]

    The Cauchy problem for general elliptic equations of second order is considered. In a previous paper (Berntsson et al. in Inverse Probl Sci Eng 26(7):1062–1078, 2018), it was suggested that the alternating iterative algorithm suggested by Kozlov and Maz’ya can be convergent, even for large wavenumbers k2, in the Helmholtz equation, if the Neumann boundary conditions are replaced by Robin conditions. In this paper, we provide a proof that shows that the Dirichlet–Robin alternating algorithm is indeed convergent for general elliptic operators provided that the parameters in the Robin conditions are chosen appropriately. We also give numerical experiments intended to investigate the precise behaviour of the algorithm for different values of k2 in the Helmholtz equation. In particular, we show how the speed of the convergence depends on the choice of Robin parameters.

    Download full text (pdf)
    fulltext
  • 3.
    Amsallem, David
    et al.
    Department of Aeronautics and Astronautics, Stanford University, Stanford, CA 94305-4035, USA.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Energy Stable Model Reduction of Neurons by Non-negative Discrete Empirical Interpolation2016In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 38, no 2, p. B297-B326Article in journal (Refereed)
    Abstract [en]

    The accurate and fast prediction of potential propagation in neuronal networks is of prime importance in neurosciences. This work develops a novel structure-preserving model reduction technique to address this problem based on Galerkin projection and nonnegative operator approximation. It is first shown that the corresponding reduced-order model is guaranteed to be energy stable, thanks to both the structure-preserving approach that constructs a distinct reduced-order basis for each cable in the network and the preservation of nonnegativity. Furthermore, a posteriori error estimates are provided, showing that the model reduction error can be bounded and controlled. Finally, the application to the model reduction of a large-scale neuronal network underlines the capability of the proposed approach to accurately predict the potential propagation in such networks while leading to important speedups.

    Download full text (pdf)
    fulltext
  • 4.
    Amsallem, David
    et al.
    Department of Aeronautics and Astronautics, Stanford University, Stanford, USA.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    High-order accurate difference schemes for the Hodgkin-Huxley equations2012Report (Other academic)
    Abstract [en]

    A novel approach for simulating potential propagation in neuronal branches with high accuracy is developed. The method relies on high-order accurate dierence schemes using the Summation-By-Parts operators with weak boundary and interface conditions applied to the Hodgkin-Huxley equations. This work is the rst demonstrating high accuracy for that equation. Several boundary conditions are considered including the non-standard one accounting for the soma presence, which is characterized by its own partial dierential equation. Well-posedness for the continuous problem as well as stability of the discrete approximation is proved for all the boundary conditions. Gains in terms of CPU times are observed when high-order operators are used, demonstrating the advantage of the high-order schemes for simulating potential propagation in large neuronal trees.

    Download full text (pdf)
    High-order accurate difference schmes for the Hodgkin-Huxley equations
  • 5.
    Amsallem, David
    et al.
    Department of Aeronautics and Astronautics, Stanford University, USA.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    High-order accurate difference schemes for the Hodgkin-Huxley equations2013In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 252, p. 573-590Article in journal (Refereed)
    Abstract [en]

    A novel approach for simulating potential propagation in neuronal branches with high accuracy is developed. The method relies on high-order accurate difference schemes using the Summation-By-Parts operators with weak boundary and interface conditions applied to the Hodgkin–Huxley equations. This work is the first demonstrating high accuracy for that equation. Several boundary conditions are considered including the non-standard one accounting for the soma presence, which is characterized by its own partial differential equation. Well-posedness for the continuous problem as well as stability of the discrete approximation is proved for all the boundary conditions. Gains in terms of CPU times are observed when high-order operators are used, demonstrating the advantage of the high-order schemes for simulating potential propagation in large neuronal trees.

    Download full text (pdf)
    fulltext
  • 6.
    Berg, Jens
    et al.
    Uppsala University, Department of Information Technology.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    A stable and dual consistent boundary treatment using finite differences on summation-by-parts form2012In: European Congress on Computational Methods in Applied Sciences and Engineering, Vienna University of Technology , 2012Conference paper (Other academic)
    Abstract [en]

    This paper is concerned with computing very high order accurate linear functionals from a numerical solution of a time-dependent partial differential equation (PDE). Based on finite differences on summation-by-parts form, together with a weak implementation of the boundary conditions, we show how to construct suitable boundary conditions for the PDE such that the continuous problem is well-posed and the discrete problem is stable and spatially dual consistent. These two features result in a superconvergent functional, in the sense that the order of accuracy of the functional is provably higher than that of the solution.

    Download full text (pdf)
    fulltext
  • 7.
    Berg, Jens
    et al.
    Uppsala University, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Duality based boundary conditions and dual consistent finite difference discretizations of the Navier–Stokes and Euler equations2014In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 259, p. 135-153Article in journal (Refereed)
    Abstract [en]

    In this paper we derive new farfield boundary conditions for the time-dependent Navier–Stokes and Euler equations in two space dimensions. The new boundary conditions are derived by simultaneously considering well-posedess of both the primal and dual problems. We moreover require that the boundary conditions for the primal and dual Navier–Stokes equations converge to well-posed boundary conditions for the primal and dual Euler equations.

    We perform computations with a high-order finite difference scheme on summation-by-parts form with the new boundary conditions imposed weakly by the simultaneous approximation term. We prove that the scheme is both energy stable and dual consistent and show numerically that both linear and non-linear integral functionals become superconvergent.

    Download full text (pdf)
    fulltext
  • 8.
    Berg, Jens
    et al.
    Uppsala University, Department of Information Technology.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Duality based boundary treatment for the Euler and Navier-Stokes equations2013In: AIAA Aerospace Sciences - Fluid Sciences Event, 2013, p. 1-19Conference paper (Other academic)
    Abstract [en]

    In this paper we construct well-posed boundary conditions for the compressible Euler and Navier-Stokes equations in two space dimensions. When also considering the dual equations, we show how to construct the boundary conditions so that both the primal and dual problems are well-posed. By considering the primal and dual problems simultaneously, we construct energy stable and dual consistent finite difference schemes on summation-by-  parts form with weak imposition of the boundary conditions.

    According to linear theory, the stable and dual consistent discretization can be used to compute linear integral functionals from the solution at a superconvergent rate. Here we evaluate numerically the superconvergence property for the non-linear Euler and Navier{ Stokes equations with linear and non-linear integral functionals.

    Download full text (pdf)
    fulltext
  • 9.
    Berg, Jens
    et al.
    Uppsala University, Department of Information Technology, SE-751 05, Uppsala, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    On the impact of boundary conditions on dual consistent finite difference discretizations2013In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 236, p. 41-55Article in journal (Refereed)
    Abstract [en]

    In this paper we derive well-posed boundary conditions for a linear incompletely parabolic system of equations, which can be viewed as a model problem for the compressible Navier{Stokes equations. We show a general procedure for the construction of the boundary conditions such that both the primal and dual equations are wellposed.

    The form of the boundary conditions is chosen such that reduction to rst order form with its complications can be avoided.

    The primal equation is discretized using finite difference operators on summation-by-parts form with weak boundary conditions. It is shown that the discretization can be made energy stable, and that energy stability is sufficient for dual consistency.

    Since reduction to rst order form can be avoided, the discretization is significantly simpler compared to a discretization using Dirichlet boundary conditions.

    We compare the new boundary conditions with standard Dirichlet boundary conditions in terms of rate of convergence, errors and discrete spectra. It is shown that the scheme with the new boundary conditions is not only far simpler, but also has smaller errors, error bounded properties, and highly optimizable eigenvalues, while maintaining all desirable properties of a dual consistent discretization.

    Download full text (pdf)
    fulltext
  • 10.
    Berg, Jens
    et al.
    Division of Scientific Computing, Department of Information Technology, Uppsala University, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Spectral analysis of the continuous and discretized heat and advection equation on single and multiple domains2012In: Applied Numerical Mathematics, ISSN 0168-9274, E-ISSN 1873-5460, Vol. 62, no 11, p. 1620-1638Article in journal (Refereed)
    Abstract [en]

    In this paper we study the heat and advectionequation in single and multipledomains. The equations are discretized using a second order accurate finite difference method on Summation-By-Parts form with weak boundary and interface conditions. We derive analytic expressions for the spectrum of the continuous problem and for their corresponding discretization matrices.

    It is shown how the spectrum of the singledomain operator is contained in the multi domain operator spectrum when artificial interfaces are introduced. The interface treatments are posed as a function of one parameter, and the impact on the spectrum and discretization error is investigated as a function of this parameter. Finally we briefly discuss the generalization to higher order accurate schemes.

    Download full text (pdf)
    fulltext
  • 11.
    Berg, Jens
    et al.
    Uppsala University, Department of Information Technology, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Superconvergent functional output for time-dependent problems using finite differences on summation-by-parts form2012In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 231, no 20, p. 6846-6860Article in journal (Refereed)
    Abstract [en]

    Finitedifference operators satisfying the summation-by-parts (SBP) rules can be used to obtain high order accurate, energy stable schemes for time-dependent partial differential equations, when the boundary conditions are imposed weakly by the simultaneous approximation term (SAT).

    In general, an SBP-SAT discretization is accurate of order p + 1 with an internal accuracy of 2p and a boundary accuracy of p. Despite this, it is shown in this paper that any linear functional computed from the time-dependent solution, will be accurate of order 2p when the boundary terms are imposed in a stable and dual consistent way.

    The method does not involve the solution of the dual equations, and superconvergent functionals are obtained at no extra computational cost. Four representative model problems are analyzed in terms of convergence and errors, and it is shown in a systematic way how to derive schemes which gives superconvergentfunctionaloutputs.

    Download full text (pdf)
    fulltext
  • 12.
    Berntsson, Fredrik
    Linköping University, Department of Mathematics, Computationel Mathematics. Linköping University, The Institute of Technology.
    An Inverse Heat Conduction Problem and Improving Shielded Thermocouple Accuracy2012In: Numerical Heat Transfer, Part A Applications, ISSN 1040-7782, E-ISSN 1521-0634, Vol. 61, no 10, p. 754-763Article in journal (Refereed)
    Abstract [en]

    A shielded thermocouple is a measurement device used for monitoring the temperature in chemically, or mechanically, hostile environments. The sensitive parts of the thermocouple are protected by a shielding layer. In order to improve the accuracy of the measurement device, we study an inverse heat conduction problem where the temperature on the surface of the shielding layer is sought, given measured temperatures in the interior of the thermocouple. The procedure is well suited for real-time applications where newly collected data is continuously used to compute current estimates of the surface temperature. Mathematically we can formulate the problem as a Cauchy problem for the heat equation, in cylindrical coordinates, where data is given along the line r = r 1 and the solution is sought at r 1 < r ≤ r 2. The problem is ill-posed, in the sense that the solution (if it exists) does not depend continuously on the data. Thus, regularization techniques are needed. The ill–posedness of the problem is analyzed and a numerical method is proposed. Numerical experiments demonstrate that the proposed method works well.

  • 13.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Eldén, Lars
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Numerical Solution of Cauchy Problems for Elliptic Equations in ``Rectangle-like'' Geometries2005In: Proceedings for the FEMLAB Conference 2005, 2005Conference paper (Other academic)
    Abstract [en]

    We consider two dimensional inverse steady state heat conductionproblems in complex geometries. The coefficients of the elliptic equation are assumed to be non-constant. Cauchy data are given on onepart of the boundary and we want to find the solution in the wholedomain. The problem is ill--posed in the sense that the solution doesnot depend continuously on the data.

    Using an orthogonal coordinate transformation the domain is mappedonto a rectangle. The Cauchy problem can then be solved by replacing one derivative by a bounded approximation. The resulting well--posed problem can then be solved by a method of lines. A bounded approximation of the derivative can be obtained by differentiating a cubic spline, that approximate the function in theleast squares sense. This particular approximation of the derivativeis computationally efficient and flexible in the sense that its easy to handle different kinds of boundary conditions.This inverse problem arises in iron production, where the walls of amelting furnace are subject to physical and chemical wear. Temperature and heat--flux data are collected by several thermocouples locatedinside the walls. The shape of the interface between the molten ironand the walls can then be determined by solving an inverse heatconduction problem.  In our work we make extensive use of Femlab for creating testproblems. By using FEMLAB we solve relatively complex model problems for the purpose of creating numerical test data used for validating our methods. For the types of problems we are intressted in numerical artefacts appear, near corners in the domain, in the gradients that Femlab calculates. We demonstrate why this happen and also how we deal with the problem.

    Download full text (pdf)
    fulltext
  • 14.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Ghosh, Arpan
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Nazarov, S. A.
    St Petersburg State Univ, Russia; Inst Problems Mech Engn RAS, Russia.
    A one dimensional model of blood flow through a curvilinear artery2018In: Applied Mathematical Modelling, ISSN 0307-904X, E-ISSN 1872-8480, Vol. 63, p. 633-643Article in journal (Refereed)
    Abstract [en]

    We present a one-dimensional model describing the blood flow through a moderately curved and elastic blood vessel. We use an existing two dimensional model of the vessel wall along with Navier-Stokes equations to model the flow through the channel while taking factors, namely, surrounding muscle tissue and presence of external forces other than gravity into account. Our model is obtained via a dimension reduction procedure based on the assumption of thinness of the vessel relative to its length. Results of numerical simulations are presented to highlight the influence of different factors on the blood flow. (C) 2018 Elsevier Inc. All rights reserved.

    Download full text (pdf)
    fulltext
  • 15.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Karlsson, Matts
    Linköping University, Department of Management and Engineering, Applied Thermodynamics and Fluid Mechanics. Linköping University, Faculty of Science & Engineering. Linköping University, Center for Medical Image Science and Visualization (CMIV).
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Nazarov, Sergei A.
    St Petersburg State University, Russia; St Petersburg State Polytech University, Russia.
    A Modification to the Kirchhoff Conditions at a Bifurcation and Loss Coefficients2018Report (Other academic)
    Abstract [en]

    One dimensional models for fluid flow in tubes are frequently used tomodel complex systems, such as the arterial tree where a large numberof vessels are linked together at bifurcations. At the junctions transmission conditions are needed. One popular option is the classic Kirchhoffconditions which means conservation of mass at the bifurcation andprescribes a continuous pressure at the joint.

    In reality the boundary layer phenomena predicts fast local changesto both velocity and pressure inside the bifurcation. Thus it is not appropriate for a one dimensional model to assume a continuous pressure. In this work we present a modification to the classic Kirchhoff condi-tions, with a symmetric pressure drop matrix, that is more suitable forone dimensional flow models. An asymptotic analysis, that has beencarried out previously shows that the new transmission conditions hasen exponentially small error.

    The modified transmission conditions take the geometry of the bifurcation into account and can treat two outlets differently. The conditions can also be written in a form that is suitable for implementationin a finite difference solver. Also, by appropriate choice of the pressuredrop matrix we show that the new transmission conditions can producehead loss coefficients similar to experimentally obtained ones.

    Download full text (pdf)
    fulltext
  • 16.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Karlsson, Matts
    Linköping University, Department of Management and Engineering, Applied Thermodynamics and Fluid Mechanics. Linköping University, Faculty of Science & Engineering.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Nazarov, Sergey A.
    St Petersburg State University, St Petersburg State Polytechnical University, and Institute of Problems of Mechanical Engineering RAS, Russia..
    A one-dimensional model of a false aneurysm2017In: International Journal of Research in Engineering and Science (IJRES), ISSN 2320-9356, Vol. 5, no 6, p. 61-73Article in journal (Refereed)
    Abstract [en]

     A false aneurysm is a hematoma, i.e. collection ofblood outside of a blood vessel, that forms due to a hole  in the wall of an artery . This represents a serious medical condition that needs to be monitored and, under certain conditions, treatedurgently. In this work a one-dimensional model of a false aneurysm isproposed. The new model is based on a one-dimensional model of anartery previously presented by the authors and it takes into accountthe interaction between the hematoma  and the surrounding musclematerial. The model equations are derived  using rigorous asymptoticanalysis for the case of a simplified geometry.   Even though the model is simple it still supports a realisticbehavior for the system consisting of the vessel and the  hematoma. Using numerical simulations we illustrate the behavior ofthe model. We also investigate the effect  of changing the size of the hematoma. The simulations show that ourmodel can reproduce realistic solutions. For instance we show thetypical strong pulsation of an aneurysm by blood entering the hematoma during the work phase of the cardiac cycle, and the blood returning tothe vessel during the resting phase. Also we show that the aneurysmgrows  if the pulse rate is increased due to, e.g., a higher work load. 

    Download full text (pdf)
    A one-dimensional model of a false aneurysm
  • 17.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Kozlov, Vladimir A.
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, The Institute of Technology.
    Mpinganzima, Lydie
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, The Institute of Technology. University of Rwanda.
    Turesson, Bengt-Ove
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, The Institute of Technology.
    Numerical Solution of the Cauchy Problem for the Helmholtz Equation2014Report (Other academic)
    Abstract [en]

    The Cauchy problem for the Helmholtz equation appears in applications related to acoustic or electromagnetic wave phenomena. The problem is ill–posed in the sense that the solution does not depend on the data in a stable way. In this paper we give a detailed study of the problem. Specifically we investigate how the ill–posedness depends on the shape of the computational domain and also on the wave number. Furthermore, we give an overview over standard techniques for dealing with ill–posed problems and apply them to the problem.

    Download full text (pdf)
    Numerical Solution of the Cauchy Problem for the Helmholtz Equation
  • 18.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Mpinganzima, L.
    University of Rwanda, Rwanda.
    Turesson, Bengt-Ove
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Iterative Tikhonov regularization for the Cauchy problem for the Helmholtz equation2017In: Computers and Mathematics with Applications, ISSN 0898-1221, E-ISSN 1873-7668, Vol. 73, no 1, p. 163-172Article in journal (Refereed)
    Abstract [en]

    The Cauchy problem for the Helmholtz equation appears in various applications. The problem is severely ill-posed and regularization is needed to obtain accurate solutions. We start from a formulation of this problem as an operator equation on the boundary of the domain and consider the equation in (H-1/2)* spaces. By introducing an artificial boundary in the interior of the domain we obtain an inner product for this Hilbert space in terms of a quadratic form associated with the Helmholtz equation; perturbed by an integral over the artificial boundary. The perturbation guarantees positivity property of the quadratic form. This inner product allows an efficient evaluation of the adjoint operator in terms of solution of a well-posed boundary value problem for the Helmholtz equation with transmission boundary conditions on the artificial boundary. In an earlier paper we showed how to take advantage of this framework to implement the conjugate gradient method for solving the Cauchy problem. In this work we instead use the Conjugate gradient method for minimizing a Tikhonov functional. The added penalty term regularizes the problem and gives us a regularization parameter that can be used to easily control the stability of the numerical solution with respect to measurement errors in the data. Numerical tests show that the proposed algorithm works well. (C) 2016 Elsevier Ltd. All rights reserved.

    Download full text (pdf)
    fulltext
  • 19.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, The Institute of Technology.
    Mpinganzima, Lydie
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, The Institute of Technology.
    Turesson, Bengt-Ove
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, The Institute of Technology.
    An alternating iterative procedure for the Cauchy problem for the Helmholtz equation2014In: Inverse Problems in Science and Engineering, ISSN 1741-5977, E-ISSN 1741-5985, Vol. 22, no 1, p. 45-62Article in journal (Refereed)
    Abstract [en]

    We present a modification of the alternating iterative method, which was introduced by V.A. Kozlov and V. Maz’ya in for solving the Cauchy problem for the Helmholtz equation in a Lipschitz domain. The method is implemented numerically using the finite difference method.

    Download full text (pdf)
    fulltext
  • 20.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Mpinganzima, Lydie
    Univ Rwanda, Rwanda.
    Turesson, Bengt-Ove
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Robin-Dirichlet algorithms for the Cauchy problem for the Helmholtz equation2018In: Inverse Problems in Science and Engineering, ISSN 1741-5977, E-ISSN 1741-5985, Vol. 26, no 7, p. 1062-1078Article in journal (Refereed)
    Abstract [en]

    The Cauchy problem for the Helmholtz equation is considered. It was demonstrated in a previous paper by the authors that the alternating algorithm suggested by V.A. Kozlov and V.G. Mazya does not converge for large wavenumbers k in the Helmholtz equation. Here, we present some simple modifications of the algorithm which may restore the convergence. They consist of the replacement of the Neumann-Dirichlet iterations by the Robin-Dirichlet ones which repairs the convergence for less than the first Dirichlet-Laplacian eigenvalue. In order to treat large wavenumbers, we present an algorithm based on iterative solution of Robin-Dirichlet boundary value problems in a sufficiently narrow border strip. Numerical implementations obtained using the finite difference method are presented. The numerical results illustrate that the algorithms suggested in this paper, produce convergent iterative sequences.

  • 21.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Kozlov, Vladimir
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Wokiyi, Dennis
    Makerere Univ, Uganda.
    Solvability of a non-linear Cauchy problem for an elliptic equation2019In: International Journal of Computer Mathematics, ISSN 0020-7160, E-ISSN 1029-0265, Vol. 96, no 12, p. 2317-2333Article in journal (Refereed)
    Abstract [en]

    We study a non-linear operator equation originating from a Cauchy problem for an elliptic equation. The problem appears in applications where surface measurements are used to calculate the temperature below the earth surface. The Cauchy problem is ill-posed and small perturbations to the used data can result in large changes in the solution. Since the problem is non-linear certain assumptions on the coefficients are needed. We reformulate the problem as an non-linear operator equation and show that under suitable assumptions the operator is well-defined. The proof is based on making a change of variables and removing the non-linearity from the leading term of the equation. As a part of the proof we obtain an iterative procedure that is convergent and can be used for evaluating the operator. Numerical results show that the proposed procedure converges faster than a simple fixed point iteration for the equation in the the original variables.

    Download full text (pdf)
    fulltext
  • 22.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Ohlson, Martin
    Linköping University, Department of Mathematics, Mathematical Statistics . Linköping University, Faculty of Science & Engineering.
    More on Estimation of Banded and Banded Toeplitz Covariance Matrices2017Report (Other academic)
    Abstract [en]

    In this paper we consider two different linear covariance structures, e.g., banded and bended Toeplitz, and how to estimate them using different methods, e.g., by minimizing different norms.

    One way to estimate the parameters in a linear covariance structure is to use tapering, which has been shown to be the solution to a universal least squares problem. We know that tapering not always guarantee the positive definite constraints on the estimated covariance matrix and may not be a suitable method. We propose some new methods which preserves the positive definiteness and still give the correct structure.

    More specific we consider the problem of estimating parameters of a multivariate normal p–dimensional random vector for (i) a banded covariance structure reflecting m–dependence, and (ii) a banded Toeplitz covariance structure.

    Download full text (pdf)
    fulltext
  • 23.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Orlof, Anna
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Thim, Johan
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Error Estimation for Eigenvalues of Unbounded Linear Operators and an Application to Energy Levels in Graphene Quantum Dots2017In: Numerical Functional Analysis and Optimization, ISSN 0163-0563, E-ISSN 1532-2467, Vol. 38, no 3, p. 293-305Article in journal (Refereed)
    Abstract [en]

    The eigenvalue problem for linear differential operators is important since eigenvalues correspond to the possible energy levels of a physical system. It is also important to have good estimates of the error in the computed eigenvalues. In this work, we use spline interpolation to construct approximate eigenfunctions of a linear operator using the corresponding eigenvectors of a discretized approximation of the operator. We show that an error estimate for the approximate eigenvalues can be obtained by evaluating the residual for an approximate eigenpair. The interpolation scheme is selected in such a way that the residual can be evaluated analytically. To demonstrate that the method gives useful error bounds, we apply it to a problem originating from the study of graphene quantum dots where the goal was to investigate the change in the spectrum from incorporating electron–electron interactions in the potential.

    Download full text (pdf)
    fulltext
  • 24.
    Berntsson, Fredrik
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Chen, Lin
    State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China.
    Xu, Tao
    State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China; CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing, China.
    Wokiyi, Dennis
    Linköping University, Department of Mathematics, Mathematics and Applied Mathematics. Linköping University, Faculty of Science & Engineering. Department of Mathematics, Makerere University, Kampala, Uganda.
    An efficient regularization method for a large scale ill-posed geothermal problem2017In: Computers & Geosciences, ISSN 0098-3004, E-ISSN 1873-7803, Vol. 105, p. 1-9Article in journal (Refereed)
    Abstract [en]

    The inverse geothermal problem consists of estimating the temperature distribution below the earth's surface using measurements on the surface. The problem is important since temperature governs a variety of geologic processes, including the generation of magmas and the deformation style of rocks. Since the thermal properties of rocks depend strongly on temperature the problem is non-linear.

    The problem is formulated as an ill-posed operator equation, where the righthand side is the heat-flux at the surface level. Since the problem is ill-posed regularization is needed. In this study we demonstrate that Tikhonov regularization can be implemented efficiently for solving the operator equation. The algorithm is based on having a code for solving a well-posed problem related to the above mentioned operator. The algorithm is designed in such a way that it can deal with both 2D and 3D calculations.

    Numerical results, for 2D domains, show that the algorithm works well and the inverse problem can be solved accurately with a realistic noise level in the surface data.

    Download full text (pdf)
    fulltext
  • 25.
    Björck, Åke
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Numerical methods in matrix computations2014 (ed. 1)Book (Refereed)
  • 26.
    Björck, Åke
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Stability of Two Direct Methods for Bidiagonalization and Partial Least Squares2014In: SIAM Journal on Matrix Analysis and Applications, ISSN 0895-4798, E-ISSN 1095-7162, Vol. 35, no 1, p. 279-291Article in journal (Refereed)
    Abstract [en]

    The partial least squares (PLS) method computes a sequence of approximate solutions x(k) is an element of K-k (A(T) A, A(T) b), k = 1, 2, ..., to the least squares problem min(x) parallel to Ax - b parallel to(2). If carried out to completion, the method always terminates with the pseudoinverse solution x(dagger) = A(dagger)b. Two direct PLS algorithms are analyzed. The first uses the Golub-Kahan Householder algorithm for reducing A to upper bidiagonal form. The second is the NIPALS PLS algorithm, due to Wold et al., which is based on rank-reducing orthogonal projections. The Householder algorithm is known to be mixed forward-backward stable. Numerical results are given, that support the conjecture that the NIPALS PLS algorithm shares this stability property. We draw attention to a flaw in some descriptions and implementations of this algorithm, related to a similar problem in Gram-Schmidt orthogonalization, that spoils its otherwise excellent stability. For large-scale sparse or structured problems, the iterative algorithm LSQR is an attractive alternative, provided an implementation with reorthogonalization is used.

    Download full text (pdf)
    fulltext
  • 27.
    Björck, Åke
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Eldén, Lars
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Axel Ruhe 1942-20152015In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 55, no 3, p. 621-623Article in journal (Other academic)
    Abstract [en]

    Axel Ruhe passed away April 4, 2015. He was cross-country-skiing with friends in the Swedish mountains when after 21 km he suddenly died. He is survived by his wife Gunlaug and three children from his first marriage....

  • 28.
    Björck, Åke
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Indahl, Ulf G.
    Norwegian University of Life Science, Norway.
    Fast and stable partial least squares modelling: A benchmark study with theoretical comments2017In: Journal of Chemometrics, ISSN 0886-9383, E-ISSN 1099-128X, Vol. 31, no 8, article id e2898Article in journal (Refereed)
    Abstract [en]

    Algorithms for partial least squares (PLS) modelling are placed into a sound theoretical context focusing on numerical precision and computational efficiency. NIPALS and other PLS algorithms that perform deflation steps of the predictors (X) may be slow or even computationally infeasible for sparse and/or large-scale data sets. As alternatives, we develop new versions of the Bidiag1 and Bidiag2 algorithms. These include full reorthogonalization of both score and loading vectors, which we consider to be both necessary and sufficient for numerical precision. Using a collection of benchmark data sets, these 2 new algorithms are compared to the NIPALS PLS and 4 other PLS algorithms acknowledged in the chemometrics literature. The provably stable Householder algorithm for PLS regression is taken as the reference method for numerical precision. Our conclusion is that our new Bidiag1 and Bidiag2 algorithms are themethods of choice for problems where both efficiency and numerical precision are important. When efficiency is not urgent, the NIPALS PLS and the Householder PLS are also good choices. The benchmark study shows that SIMPLS gives poor numerical precision even for a small number of factors. Further, the nonorthogonal scores PLS, direct scores PLS, and the improved kernel PLS are demonstrated to be numerically less stable than the best algorithms. PrototypeMATLAB codes are included for the 5 PLS algorithms concluded to be numerically stable on our benchmark data sets. Other aspects of PLS modelling, such as the evaluation of the regression coefficients, are also analyzed using techniques from numerical linear algebra.

  • 29.
    Bohm, Marvin
    et al.
    University of Cologne, Cologne, Germany.
    Schermeng, Sven
    University of Cologne, Cologne, Germany.
    Winters, Andrew Ross
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Gassner, Gregor J
    University of Cologne, Cologne, Germany.
    Jacobs, Gustaaf B
    San Diego State University, San Diego, USA.
    Multi-element SIAC Filter for Shock Capturing Applied to High-Order Discontinuous Galerkin Spectral Element Methods2019In: Journal of Scientific Computing, ISSN 0885-7474, E-ISSN 1573-7691, ISSN 0885-7474, Vol. 81, no 2, p. 820-844Article in journal (Refereed)
    Abstract [en]

    We build a multi-element variant of the smoothness increasing accuracy conserving (SIAC) shock capturing technique proposed for single element spectral methods by Wissink et al. (J Sci Comput 77:579–596, 2018). In particular, the baseline scheme of our method is the nodal discontinuous Galerkin spectral element method (DGSEM) for approximating the solution of systems of conservation laws. It is well known that high-order methods generate spurious oscillations near discontinuities which can develop in the solution for nonlinear problems, even when the initial data is smooth. We propose a novel multi-element SIAC filtering technique applied to the DGSEM as a shock capturing method. We design the SIAC filtering such that the numerical scheme remains high-order accurate and that the shock capturing is applied adaptively throughout the domain. The shock capturing method is derived for general systems of conservation laws. We apply the novel SIAC filter to the two-dimensional Euler and ideal magnetohydrodynamics equations to several standard test problems with a variety of boundary conditions.

    Download full text (pdf)
    fulltext
  • 30.
    Carpenter, Mark H
    et al.
    Computational Aerosciences Branch, NASA Langley Research Center, Hampton, VA 23681, USA.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Gottlieb, David
    cDivision of Applied Mathematics, Brown University, Providence, RI 02912, USA.
    Corrigendum to “A stable and conservative interface treatment of arbitrary spatial accuracy” [J.Comput.Phys.148(1999)341–365]2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 351, p. 534-Article in journal (Other academic)
    Download full text (pdf)
    fulltext
  • 31.
    Changfoot, Donovan M.
    et al.
    University of Cape Town, Cape Town, South Africa.
    Malan, Arnaud G
    University of Cape Town, Cape Town, South Africa.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Hybrid Computational-Fluid-Dynamics Platform to Investigate Aircraft Trailing Vortices2019In: Journal of Aircraft, ISSN 0021-8669, E-ISSN 1533-3868, Vol. 56, no 1, p. 344-355Article in journal (Refereed)
    Abstract [en]

    This paper outlines the development of a parallel three-dimensional hybrid finite volume finite difference capability. The specific application area under consideration is modeling the trailing vortices shed from the wings of aircraft under transonic flight conditions. For this purpose, the Elemental finite volume code is employed in the vicinity of the aircraft, whereas the ESSENSE finite difference software is employed to accurately resolve the trailing vortices. The former method is spatially formally second-order, and the latter is set to sixth-order accuracy. The coupling of the two methods is achieved in a stable manner through the use of summation-by-parts operators and weak imposition of boundary conditions using simultaneous approximation terms. The developed hybrid solver is successfully validated against an analytical test case. This is followed by demonstrating the ability to model the flowfield, including trailing vortex structures, around the NASA Common Research Model under transonic flow conditions. The interface treatment is shown to describe the intersecting vortices in a smooth manner. In addition, insights gained in resolving the vortices include violation of underlying assumptions of analytical vortex modeling methods.

    Download full text (pdf)
    fulltext
  • 32.
    Compagnoni, Marco
    et al.
    Politecn Milan, Italy.
    Notari, Roberto
    Politecn Milan, Italy.
    Ruggiu, Andrea Alessandro
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Antonacci, Fabio
    Politecn Milan, Italy.
    Sarti, Augusto
    Politecn Milan, Italy.
    The Algebro-geometric Study of Range Maps2017In: Journal of nonlinear science, ISSN 0938-8974, E-ISSN 1432-1467, Vol. 27, no 1, p. 99-157Article in journal (Refereed)
    Abstract [en]

    Localizing a radiant source is a problem of great interest to many scientific and technological research areas. Localization based on range measurements is at the core of technologies such as radar, sonar and wireless sensor networks. In this manuscript, we offer an in-depth study of the model for source localization based on range measurements obtained from the source signal, from the point of view of algebraic geometry. In the case of three receivers, we find unexpected connections between this problem and the geometry of Kummers and Cayleys surfaces. Our work also gives new insights into the localization based on range differences.

  • 33.
    Delorme, Yann T.
    et al.
    Faculty of Mechanical Engineering, Technion-Israel Institute of Technology, Haifa, Israel.
    Puri, Kunal
    Faculty of Mechanical Engineering, Technion-Israel Institute of Technology, Haifa, Israel.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Linders, Viktor
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Dong, Suchuan
    Department of Mathematics, Purdue University, West Lafayette, IN, USA.
    Frankel, Steven H.
    Faculty of Mechanical Engineering, Technion-Israel Institute of Technology, Haifa, Israel.
    A simple and efficient incompressible Navier-Stokes solver for unsteady complex geometry flows on truncated domains2017In: Computers & Fluids, ISSN 0045-7930, E-ISSN 1879-0747, Vol. 150, p. 84-94Article in journal (Refereed)
    Abstract [en]

    Incompressible Navier-Stokes solvers based on the projection method often require an expensive numerical solution of a Poisson equation for a pressure-like variable. This often involves linear system solvers based on iterative and multigrid methods which may limit the ability to scale to large numbers of processors. The artificial compressibility method (ACM) [6] introduces a time derivative of the pressure into the incompressible form of the continuity equation creating a coupled closed hyperbolic system that does not require a Poisson equation solution and allows for explicit time-marching and localized stencil numerical methods. Such a scheme should theoretically scale well on large numbers of CPUs, GPU'€™s, or hybrid CPU-GPU architectures. The original ACM was only valid for steady flows and dual-time stepping was often used for time-accurate simulations. Recently, Clausen [7] has proposed the entropically damped artificial compressibility (EDAC) method which is applicable to both steady and unsteady flows without the need for dual-time stepping. The EDAC scheme was successfully tested with both a finite-difference MacCormack'€™s method for the two-dimensional lid driven cavity and periodic double shear layer problem and a finite-element method for flow over a square cylinder, with scaling studies on the latter to large numbers of processors. In this study, we discretize the EDAC formulation with a new optimized high-order centered finite-difference scheme and an explicit fourth-order Runge-€“Kutta method. This is combined with an immersed boundary method to efficiently treat complex geometries and a new robust outflow boundary condition to enable higher Reynolds number simulations on truncated domains. Validation studies for the Taylor-€“Green Vortex problem and the lid driven cavity problem in both 2D and 3D are presented. An eddy viscosity subgrid-scale model is used to enable large eddy simulations for the 3D cases. Finally, an application to flow over a sphere is presented to highlight the boundary condition and performance comparisons to a traditional incompressible Navier-€“Stokes solver is shown for the 3D lid driven cavity. Overall, the combined EDAC formulation and discretization is shown to be both effective and affordable.

    Download full text (pdf)
    A simple and efficient incompressible Navier-Stokes solver for unsteady complex geometry flows on truncated domains
  • 34.
    Eldén, Lars
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Computing Frechet derivatives in partial least squares regression2015In: Linear Algebra and its Applications, ISSN 0024-3795, E-ISSN 1873-1856, Vol. 473, p. 316-338Article in journal (Refereed)
    Abstract [en]

    Partial least squares is a common technique for multivariate regression. The pro- cedure is recursive and in each step basis vectors are computed for the explaining variables and the solution vectors. A linear model is fitted by projection onto the span of the basis vectors. The procedure is mathematically equivalent to Golub-Kahan bidiagonalization, which is a Krylov method, and which is equiv- alent to a pair of matrix factorizations. The vectors of regression coefficients and prediction are non-linear functions of the right hand side. An algorithm for computing the Frechet derivatives of these functions is derived, based on perturbation theory for the matrix factorizations. From the Frechet derivative of the prediction vector one can compute the number of degrees of freedom, which can be used as a stopping criterion for the recursion. A few numerical examples are given.

    Download full text (pdf)
    fulltext
  • 35.
    Eldén, Lars
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Editorial Material: Preface in BIT NUMERICAL MATHEMATICS, vol 57, issue 1, pp2017In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 57, no 1Article in journal (Other academic)
    Abstract [en]

    n/a

  • 36.
    Eldén, Lars
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Editorial Material: Preface to BIT 55:4 in BIT NUMERICAL MATHEMATICS, vol 55, issue 4, pp 897-8992015In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 55, no 4, p. 897-899Article in journal (Other academic)
    Abstract [en]

    n/a

  • 37.
    Eldén, Lars
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Preface to BIT 56:4 in BIT NUMERICAL MATHEMATICS2016In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 56, no 4, p. 1163-1164Article in journal (Other academic)
    Abstract [en]

    n/a

  • 38.
    Eldén, Lars
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Ahmadi-Asl, Salman
    Skolkovo Inst Sci and Technol Skoltech, Russia.
    Solving bilinear tensor least squares problems and application to Hammerstein identification2019In: Numerical Linear Algebra with Applications, ISSN 1070-5325, E-ISSN 1099-1506, Vol. 26, no 2, article id e2226Article in journal (Refereed)
    Abstract [en]

    Bilinear tensor least squares problems occur in applications such as Hammerstein system identification and social network analysis. A linearly constrained problem of medium size is considered, and nonlinear least squares solvers of Gauss-Newton-type are applied to numerically solve it. The problem is separable, and the variable projection method can be used. Perturbation theory is presented and used to motivate the choice of constraint. Numerical experiments with Hammerstein models and random tensors are performed, comparing the different methods and showing that a variable projection method performs best.

    Download full text (pdf)
    fulltext
  • 39.
    Eldén, Lars
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Merkel, Magnus
    Linköping University, Department of Computer and Information Science, Human-Centered systems. Linköping University, The Institute of Technology.
    Ahrenberg, Lars
    Linköping University, Department of Computer and Information Science, Human-Centered systems. Linköping University, The Institute of Technology.
    Fagerlund, Martin
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Computing Semantic Clusters by Semantic Mirroring and Spectral Graph Partitioning2013In: Mathematics in Computer Science, ISSN 1661-8270, Vol. 7, p. 293-313Article in journal (Refereed)
    Abstract [en]

    Using the technique of semantic mirroring a graph is obtained that represents words and their translationsfrom a parallel corpus or a bilingual lexicon. The connectedness of the graph holds information about the semanticrelations of words that occur in the translations. Spectral graph theory is used to partition the graph, which leadsto a grouping of the words in different clusters. We illustrate the method using a small sample of seed words froma lexicon of Swedish and English adjectives and discuss its application to computational lexical semantics andlexicography.

  • 40.
    Eldén, Lars
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Trendafilov, Nickolay
    Open Univ, England.
    Semi-sparse PCA2019In: Psychometrika, ISSN 0033-3123, E-ISSN 1860-0980, Vol. 84, no 1, p. 164-185Article in journal (Refereed)
    Abstract [en]

    It is well known that the classical exploratory factor analysis (EFA) of data with more observations than variables has several types of indeterminacy. We study the factor indeterminacy and show some new aspects of this problem by considering EFA as a specific data matrix decomposition. We adopt a new approach to the EFA estimation and achieve a new characterization of the factor indeterminacy problem. A new alternative model is proposed, which gives determinate factors and can be seen as a semi-sparse principal component analysis (PCA). An alternating algorithm is developed, where in each step a Procrustes problem is solved. It is demonstrated that the new model/algorithm can act as a specific sparse PCA and as a low-rank-plus-sparse matrix decomposition. Numerical examples with several large data sets illustrate the versatility of the new model, and the performance and behaviour of its algorithmic implementation.

    Download full text (pdf)
    fulltext
  • 41.
    Elfving, Tommy
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Row and column based Iterations2018In: Applied Analysis and Optimization, ISSN 2432-1656, Vol. 2, no 2, p. 219-236Article in journal (Other academic)
    Abstract [en]

    We survey old and more recent results on row- and column-action iterative methods for solving ill-conditioned linear systems. Our main application is in X-ray tomography problems with missing and/or noisy data. We consider the stationary case with cyclic control. A unified framework is presented the use of which allows deriving both necessary and sufficient convergence conditions for many of the methods presented.

  • 42.
    Elfving, Tommy
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Christian Hansen, Per
    Technical University of Denmark, Denmark.
    Nikazad, Touraj
    Iran University of Science and Technology, Iran.
    Correction: Convergence analysis for column-action methods in image reconstruction (vol 74, pg 905, 2017)2017In: Numerical Algorithms, ISSN 1017-1398, E-ISSN 1572-9265, Vol. 74, no 3, p. 925-925Article in journal (Other academic)
    Abstract [en]

    n/a

  • 43.
    Elfving, Tommy
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Christian Hansen, Per
    Technical University of Denmark, Denmark .
    Nikazad, Touraj
    Iran University of Science and Technology, Iran .
    Semi-convergence properties of Kaczmarzs method2014In: Inverse Problems, ISSN 0266-5611, E-ISSN 1361-6420, Vol. 30, no 5, p. 055007-Article in journal (Refereed)
    Abstract [en]

    Kaczmarzs method-sometimes referred to as the algebraic reconstruction technique-is an iterative method that is widely used in tomographic imaging due to its favorable semi-convergence properties. Specifically, when applied to a problem with noisy data, during the early iterations it converges very quickly toward a good approximation of the exact solution, and thus produces a regularized solution. While this property is generally accepted and utilized, there is surprisingly little theoretical justification for it. The purpose of this paper is to present insight into the semi-convergence of Kaczmarzs method as well as its projected counterpart (and their block versions). To do this we study how the data errors propagate into the iteration vectors and we derive upper bounds for this noise propagation. Our bounds are compared with numerical results obtained from tomographic imaging.

  • 44.
    Elfving, Tommy
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Hansen, Per Christian
    Tech Univ Denmark, Denmark.
    UNMATCHED PROJECTOR/BACKPROJECTOR PAIRS: PERTURBATION AND CONVERGENCE ANALYSIS2018In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 40, no 1, p. A573-A591Article in journal (Refereed)
    Abstract [en]

    In tomographic reconstruction problems it is not uncommon that there are errors in the implementation of the forward projector and/or the backprojector, and hence we encounter a so-called unmatched projektor/backprojector pair. Consequently, the matrices that represent the two projectors are not each others transpose. Surprisingly, the influence of such errors in algebraic iterative reconstruction methods has received little attention in the literature. The goal of this paper is to perform a rigorous first-order perturbation analysis of the minimization problems underlying the algebraic methods in order to understand the role played by the nonmatch of the matrices. We also study the convergence properties of linear stationary iterations based on unmatched matrix pairs, leading to insight into the behavior of some important row-and column-oriented algebraic iterative methods. We conclude with numerical examples that illustrate the perturbation and convergence results.

  • 45.
    Elfving, Tommy
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Hansen, Per Christian
    Technical University of Denmark, Denmark.
    Nikazad, Touraj
    Iran University of Science and Technology, Iran.
    Convergence analysis for column-action methods in image reconstruction2017In: Numerical Algorithms, ISSN 1017-1398, E-ISSN 1572-9265, Vol. 74, no 3, p. 905-924Article in journal (Refereed)
    Abstract [en]

    Column-oriented versions of algebraic iterative methods are interesting alternatives to their row-version counterparts: they converge to a least squares solution, and they provide a basis for saving computational work by skipping small updates. In this paper we consider the case of noise-free data. We present a convergence analysis of the column algorithms, we discuss two techniques (loping and flagging) for reducing the work, and we establish some convergence results for methods that utilize these techniques. The performance of the algorithms is illustrated with numerical examples from computed tomography.

  • 46.
    Elfving, Tommy
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Hansen, Per Christian
    Technical University of Denmark, Lyngby, Denmark.
    Nikazad, Touraj
    Iran University of Science and Technology, Narmak, Tehran, Iran .
    Semiconvergence and Relaxation Parameters for Projected SIRT Algorithms2012In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 34, no 4, p. A2000-A2017Article in journal (Refereed)
    Abstract [en]

    We give a detailed study of the semiconvergence behavior of projected nonstationary simultaneous iterative reconstruction technique (SIRT) algorithms, including the projected Landweber algorithm. We also consider the use of a relaxation parameter strategy, proposed recently for the standard algorithms, for controlling the semiconvergence of the projected algorithms. We demonstrate the semiconvergence and the performance of our strategies by examples taken from tomographic imaging.

    Download full text (pdf)
    fulltext
  • 47.
    Eliasson, Peter
    et al.
    Department of Aeronautics and Autonomous Systems, FOI, Swedish Defense Research Agency, SE-164 90, Stockholm, Sweden .
    Kupiainen, Marco
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Higher Order Accurate Solutions for Flow in a Cavity: Experiences and Lessons Learned2015In: Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2014 / [ed] Mejdi Azaïez, Henda El Fekih, Jan S. Hesthaven, Springer, 2015, p. 189-196Chapter in book (Refereed)
    Abstract [en]

    Experiences from using a higher order accurate finite difference multiblock solver to compute the time dependent flow over a cavity is summarized. The work has been carried out as part of a work in a European project called IDIHOM in a collaboration between the Swedish Defense Research Agency (FOI) and University of Linköping (LiU). The higher order code is based on Summation By Parts operators combined with the Simultaneous Approximation Term approach for boundary and interface conditions. The spatial accuracy of the code is verified by calculations over a cyclinder by monitoring the decay of the errors of known wall quantities as the grid is refined. The focus is on the validation for a test case of transonic flow over a rectangular cavity with hybrid RANS/LES calculations. The results are compared to reference numerical results from a second order finite volume code as well as with experimental results with a good overall agreement between the results.

    Download full text (pdf)
    Higher Order Accurate Solutions for Flow in a Cavity: Experiences and Lessons Learned
  • 48.
    Eliasson, Peter
    et al.
    FOI, Swedish Defence Research Agency, SE-16490 Stockholm, Sweden.
    Lundquist, Tomas
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    A global time integration approach for realistic unsteady flow computations2016Conference paper (Refereed)
    Abstract [en]

    A novel time integration approach is explored for unsteady flow computations. It is a multi-block formulation in time where one solves for all time levels within a block simultaneously. The time discretization within a block is based on the summation-by-parts (SBP) technique in time combined with the simultaneous-approximation-term (SAT) technique for imposing the initial condition. The approach is implicit, unconditionally stable and can be made high order accurate in time. The implicit system is solved by a dual time stepping technique. The technique has been implemented in a flow solver for unstructured grids and applied to an unsteady flow problem with vortex shedding over a cylinder. Four time integration approaches being 2nd to 5th order accurate in time are evaluated and compared to the conventional 2nd order backward difference (BDF2) method and a 4th order diagonally implicit Runge-Kutta scheme (ESDIRK64). The obtained orders of accuracy are higher than expected and correspond to the accuracy in the interior of the blocks, up to 8th order accuracy is obtained. The influence on the accuracy from the size of the time blocks is small. Smaller blocks are computationally more efficient though, and the efficiency increases with increased accuracy of the SBP operator and reduced size of time steps. The most accurate scheme, with a small time step and block size, is approximately as efficient as the ESDIRK64 scheme. There is a significant potential for improvements ranging from convergence acceleration techniques in dual time, alternative initialization of time blocks, and by introducing smaller time blocks based on alternative SBP operators.

    Download full text (pdf)
    A global time integration approach for realistic unsteady flow computations
  • 49.
    Eliasson, Peter
    et al.
    Swedish Defence Research Agency (FOI), Stockholm, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    The Influence of Viscous Operator and Wall Boundary Conditions on the Accuracy of the Navier-Stokes Equations2013In: AIAA Aerospace Sciences - Fluid Sciences Event, 2013, p. 1-15Conference paper (Other academic)
    Abstract [en]

    The discretization of the viscous operator in an edge-based flow solver for unstructured grids has been investigated. A compact discretization of the Laplace and thin-layer operators in the viscous terms is used with two different wall boundary conditions. Furthermore, a wide discretization of the same operators is investigated. The resulting numerical operators are all formally second order accurate in space; the wide operator has higher truncation errors. The operators are implemented weakly using a penalty formulation and are shown to be stable for a scalar model problem with given constraints on the penalty coefficients. The different operators are applied to a set of grid convergence test cases for laminar flow in two dimensions up to a large-scale three dimensional turbulent flow problem. The operators converge to the same solutions as the grids are refined with one exception where the wide operator converges to a solution with higher drag. The 2nd compact discretization, being locally more accurate at a wall boundary than the original 1st compact operator, reduces the grid dependency somewhat for most test cases. The wide operator performs very well and leads for most test cases to results with minimum spread between coarsest and finest grids. For one test case though, the wide operator has a negative influence on the steady state convergence.

    Download full text (pdf)
    fulltext
  • 50.
    Erickson, Brittany. A.
    et al.
    Department of Geological Science, San Diego State University, 5500 Campanile Drive, San Diego, California, 92182-1020..
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Stable, High Order Accurate Adaptive Schemes for Long Time, Highly Intermittent Geophysics Problems2013Report (Other academic)
    Abstract [en]

    Many geophysical phenomena are characterized by properties that evolve over a wide range of scales which introduce difficulties when attempting to model these features in one computational method. We have developed a high-order finite difference method for the elastic wave equation that is able to efficiently handle varying temporal scales in a single, stand-alone framework. We apply this method to earthquake cycle models characterized by extremely long interseismic periods interspersed with abrupt, short periods of dynamic rupture. Through the use of summation-by-parts operators and weak enforcement of boundary conditions we derive a provably stable discretization. Time stepping is achieved through the implicit θ-method which allows us to take large time steps during the intermittent period between earthquakes and adapts appropriately to fully resolve rupture.

    Download full text (pdf)
    Stable, High Order Accurate Adaptive Schemes for Long Time, Highly Intermittent Geophysics Problem
12345 1 - 50 of 243
CiteExportLink to result list
Permanent link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • oxford
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf