liu.seSearch for publications in DiVA
Change search
Refine search result
12 1 - 50 of 92
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.
    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
  • 2.
    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
  • 3.
    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
  • 4.
    Berg, Jens
    et al.
    Uppsala University, Department of Information Technology.
    Nordström, Jan
    Linköping University, The Institute of Technology. Linköping University, Department of Mathematics, Scientific Computing.
    Stable Robin solid wall boundary conditions for the Navier-Stokes equations2011In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 230, no 19, p. 7519-7532Article in journal (Refereed)
    Abstract [en]

    In this paper we prove stability of Robin solid wall boundary conditions for the compressible Navier–Stokes equations. Applications include the no-slip boundary conditions with prescribed temperature or temperature gradient and the first order slip-flow boundary conditions. The formulation is uniform and the transitions between different boundary conditions are done by a change of parameters. We give different sharp energy estimates depending on the choice of parameters.

    The discretization is done using finite differences on Summation-By-Parts form with weak boundary conditions using the Simultaneous Approximation Term. We verify convergence by the method of manufactured solutions and show computations of flows ranging from no-slip to almost full slip.

    Download full text (pdf)
    fulltext
  • 5.
    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
  • 6.
    Bohm, Marvin
    et al.
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Winters, Andrew Ross
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Gassner, Gregor J
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Derigs, Dominik
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    Hindenlang, Florian
    Max-Planck Institut für Plasmaphysik, Garching, Germany.
    Saur, Joachim
    Institut für Geophysik und Meteorologie, Universität zu Köln, Köln, Germany.
    An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and numerical verification2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716Article in journal (Refereed)
    Abstract [en]

    The first paper of this series presents a discretely entropy stable discontinuous Galerkin (DG) method for the resistive magnetohydrodynamics (MHD) equations on three-dimensional curvilinear unstructured hexahedral meshes. Compared to other fluid dynamics systems such as the shallow water equations or the compressible Navier-Stokes equations, the resistive MHD equations need special considerations because of the divergence-free constraint on the magnetic field. For instance, it is well known that for the symmetrization of the ideal MHD system as well as the continuous entropy analysis a non-conservative term proportional to the divergence of the magnetic field, typically referred to as the Powell term, must be included. As a consequence, the mimicry of the continuous entropy analysis in the discrete sense demands a suitable DG approximation of the non-conservative terms in addition to the ideal MHD terms.

    This paper focuses on the resistive MHD equations: Our first contribution is a proof that the resistive terms are symmetric and positive-definite when formulated in entropy space as gradients of the entropy variables, which enables us to show that the entropy inequality holds for the resistive MHD equations. This continuous analysis is the key for our DG discretization and guides the path for the construction of an approximation that discretely mimics the entropy inequality, typically termed entropy stability. Our second contribution is a detailed derivation and analysis of the discretization on three-dimensional curvilinear meshes. The discrete analysis relies on the summation-by-parts property, which is satisfied by the DG spectral element method (DGSEM) with Legendre-Gauss-Lobatto (LGL) nodes. Although the divergence-free constraint is included in the non-conservative terms, the resulting method has no particular treatment of the magnetic field divergence errors, which might pollute the solution quality. Our final contribution is the extension of the standard resistive MHD equations and our DG approximation with a divergence cleaning mechanism that is based on a generalized Lagrange multiplier (GLM).

    As a conclusion to the first part of this series, we provide detailed numerical validations of our DGSEM method that underline our theoretical derivations. In addition, we show a numerical example where the entropy stable DGSEM demonstrates increased robustness compared to the standard DGSEM.

  • 7.
    Burdakov, Oleg
    et al.
    Linköping University, Department of Mathematics, Optimization . Linköping University, The Institute of Technology.
    Kapyrin, Ivan
    Russian Academy of Science.
    Vassilevski, Yuri
    Russian Academy of Science.
    Monotonicity recovering and accuracy preserving optimization methods for postprocessing finite element solutions2012In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 231, no 8, p. 3126-3142Article in journal (Refereed)
    Abstract [en]

    We suggest here a least-change correction to available finite element (FE) solution. This postprocessing procedure is aimed at recovering the monotonicity and some other important properties that may not be exhibited by the FE solution. Although our approach is presented for FEs, it admits natural extension to other numerical schemes, such as finite differences and finite volumes. For the postprocessing, a priori information about the monotonicity is assumed to be available, either for the whole domain or for a subdomain where the lost monotonicity is to be recovered. The obvious requirement is that such information is to be obtained without involving the exact solution, e.g. from expected symmetries of this solution. less thanbrgreater than less thanbrgreater thanThe postprocessing is based on solving a monotonic regression problem with some extra constraints. One of them is a linear equality-type constraint that models the conservativity requirement. The other ones are box-type constraints, and they originate from the discrete maximum principle. The resulting postprocessing problem is a large scale quadratic optimization problem. It is proved that the postprocessed FE solution preserves the accuracy of the discrete FE approximation. less thanbrgreater than less thanbrgreater thanWe introduce an algorithm for solving the postprocessing problem. It can be viewed as a dual ascent method based on the Lagrangian relaxation of the equality constraint. We justify theoretically its correctness. Its efficiency is demonstrated by the presented results of numerical experiments.

    Download full text (pdf)
    fulltext
  • 8.
    Carpenter, Mark H.
    et al.
    Aerodynamic and Acoustic Methods Branch, NASA Langley Research Center, Hampton, Virginia, USA.
    Nordström, Jan
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Gottlieb, David
    Division of Applied Mathematics, Brown University, Providence, Rhode Island, USA.
    A stable and conservative interface treatment of arbitrary spatial accuracy1999In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 148, p. 341-365Article in journal (Refereed)
    Abstract [en]

    Stable and accurate interface conditions based on the SAT penalty method are derived for the linear advection–diffusion equation. The conditions are functionally independent of the spatial order of accuracy and rely only on the form of the discrete operator. We focus on high-order finite-difference operators that satisfy the summation-by-parts (SBP) property. We prove that stability is a natural consequence of the SBP operators used in conjunction with the new, penalty type, boundary conditions. In addition, we show that the interface treatments are conservative. The issue of the order of accuracy of the interface boundary conditions is clarified. New finite-difference operators of spatial accuracy up to sixth order are constructed which satisfy the SBP property. These finite-difference operators are shown to admit design accuracy (pth-order global accuracy) when (p−1)th-order stencil closures are used near the boundaries, if the physical boundary conditions and interface conditions are implemented to at leastpth-order accuracy. Stability and accuracy are demonstrated on the nonlinear Burgers' equation for a 12-subdomain problem with randomly distributed interfaces.

  • 9.
    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
  • 10.
    Derigs, Dominik
    et al.
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    Winters, Andrew Ross
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Gassner, Gregor J
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Walch, Stefanie
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    A novel averaging technique for discrete entropy-stable dissipation operators for ideal MHD2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 330, p. 624-632Article in journal (Refereed)
    Abstract [en]

    Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation of the discrete entropy. Present works in the field of entropy stable numerical schemes are concerned with thorough derivations of entropy conservative fluxes for ideal MHD. However, as we show in this work, if the dissipation operator is not constructed in a very specific way, it cannot lead to a generally stable numerical scheme. The two main findings presented in this paper are that the entropy conserving flux of Ismail & Roe can easily break down for certain initial conditions commonly found in astrophysical simulations, and that special care must be taken in the derivation of a discrete dissipation matrix for an entropy stable numerical scheme to be robust. We present a convenient novel averaging procedure to evaluate the entropy Jacobians of the ideal MHD and the compressible Euler equations that yields a discretization with favorable robustness properties.

    Download full text (pdf)
    A novel averaging technique for discrete entropy-stable dissipation operators for ideal MHD
  • 11.
    Derigs, Dominik
    et al.
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    Winters, Andrew Ross
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Gassner, Gregor J
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Walch, Stefanie
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure2016In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 317, p. 223-256Article in journal (Refereed)
    Abstract [en]

    We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code FLASH (http://flash.uchicago.edu). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in FLASH.

    Download full text (pdf)
    A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure
  • 12.
    Derigs, Dominik
    et al.
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    Winters, Andrew Ross
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Gassner, Gregor J
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Walch, Stefanie
    I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
    Bohm, Marvin
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 364, p. 420-467Article in journal (Refereed)
    Abstract [en]

    The paper presents two contributions in the context of the numerical simulation of magnetized fluid dynamics. First, we show how to extend the ideal magnetohydrodynamics (MHD) equations with an inbuilt magnetic field divergence cleaning mechanism in such a way that the resulting model is consistent with the second law of thermodynamics. As a byproduct of these derivations, we show that not all of the commonly used divergence cleaning extensions of the ideal MHD equations are thermodynamically consistent. Secondly, we present a numerical scheme obtained by constructing a specific finite volume discretization that is consistent with the discrete thermodynamic entropy. It includes a mechanism to control the discrete divergence error of the magnetic field by construction and is Galilean invariant. We implement the new high-order MHD solver in the adaptive mesh refinement code FLASH where we compare the divergence cleaning efficiency to the constrained transport solver available in FLASH (unsplit staggered mesh scheme).

    Download full text (pdf)
    Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations
  • 13.
    Eriksson, Sofia
    et al.
    Department of Information Technology, Scientific Computing, Uppsala University, SE-751 05 Uppsala, Sweden.
    Abbas, Qaisar
    Department of Information Technology, Scientific Computing, Uppsala University, SE-751 05 Uppsala, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Scientific Computing. Linköping University, The Institute of Technology.
    A stable and conservative method for locally adapting the design order of finite difference schemes2011In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 230, no 11, p. 4216-4231Article in journal (Refereed)
    Abstract [en]

    A procedure to locally change the order of accuracy of finite difference schemes is developed. The development is based on existing Summation-By-Parts operators and a weak interface treatment. The resulting scheme is proven to be accurate and stable.

     

    Numerical experiments verify the theoretical accuracy for smooth solutions. In addition, shock calculations are performed, using a scheme where the developed switching procedure is combined with the MUSCL technique.

    Download full text (pdf)
    fulltext
  • 14.
    Eriksson, Sofia
    et al.
    Linnaeus University, Växjö, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Finite difference schemes with transferable interfaces for parabolic problems2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 375, p. 935-949Article in journal (Refereed)
    Abstract [en]

    We derive a method to locally change the order of accuracy of finite difference schemes that approximate the second derivative. The derivation is based on summation-by-parts operators, which are connected at interfaces using penalty terms. At such interfaces, the numerical solution has a double representation, with one representation in each domain. We merge this double representation into a single one, yielding a new scheme with unique solution values in all grid points. The resulting scheme is proven to be stable, accurate and dual consistent.

    Download full text (pdf)
    fulltext
  • 15.
    Fisher, Travis C.
    et al.
    Computational Aerosciences Branch, NASA Langley Research Center, Hampton, VA 23681, USA.
    Carpenter, Mark H.
    School of Mechanical Engineering, Purdue University, West Lafayette, IN 47907, USA.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Yamaleev, Nail K.
    Department of Mathematics, North Carolina A&T State University, Greensboro, NC 27411, USA.
    Swanson, Charles
    Distinguished Research Associate, NASA Langley Research Center.
    Discretely Conservative Finite-Difference Formulations for Nonlinear Conservation Laws in Split Form: Theory and Boundary Conditions2013In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 234, p. 353-375Article in journal (Refereed)
    Abstract [en]

    The Lax-Wendroff theorem stipulates that a discretely conservative operator is necessary to accurately capture discontinuities. The discrete operator, however, need not be derived from the divergence form of the continuous equations. Indeed, conservation law equations that are split into linear combinations of the divergence and product rule form and then discretized using any diagonal-norm skew-symmetric summation-by-parts (SBP) spatial operator, yield discrete operators that are conservative. Furthermore, split-form, discretely conservation operators can be derived for periodic or finite-domain SBP spatial operators of any order. Examples are presented of a fourth-order, SBP finite-difference operator with second-order boundary closures. Sixth- and eighth-order constructions are derived, and are supplied in an accompanying text file.

    Download full text (pdf)
    fulltext
  • 16.
    Frenander, Hannes
    et al.
    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.
    Constructing non-reflecting boundary conditions using summation-by-parts in time2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 331, p. 38-48Article in journal (Refereed)
    Abstract [en]

    In this paper we provide a new approach for constructing non-reflecting boundary conditions. The boundary conditions are based on summation-by-parts operators and derived without Laplace transformation in time. We prove that the new non-reflecting boundary conditions yield a well-posed problem and that the corresponding numerical approximation is unconditionally stable. The analysis is demonstrated on a hyperbolic system in two space dimensions, and the theoretical results are confirmed by numerical experiments.

    Download full text (pdf)
    fulltext
  • 17.
    Gassner, Gregor J
    et al.
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Winters, Andrew Ross
    Mathematisches Institut, Universität zu Köln, Köln, Germany.
    Kopriva, David A
    Department of Mathematics, The Florida State University, Tallahassee, US.
    Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations2016In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 327, p. 39-66Article in journal (Refereed)
    Abstract [en]

    Fisher and Carpenter (High-order entropy stable finite difference schemes for non-linear conservation laws: Finite domains, Journal of Computational Physics, 252:518–557, 2013) found a remarkable equivalence of general diagonal norm high-order summation-by- parts operators to a subcell based high-order finite volume formulation. This equivalence enables the construction of provably entropy stable schemes by a specific choice of the sub-cell finite volume flux. We show that besides the construction of entropy stable high order schemes, a careful choice of subcell finite volume fluxes generates split formulations of quadratic or cubic terms. Thus, by changing the subcell finite volume flux to a specific choice, we are able to generate, in a systematic way, all common split forms of the compressible Euler advection terms, such as the Ducros splitting and the Kennedy and Gruber splitting. Although these split forms are not entropy stable, we present a systematic way to prove which of those split forms are at least kinetic energy preserving. With this, we show we construct a unified high-order split form DG framework. We investigate with three dimensional numerical simulations of the inviscid Taylor-Green vortex and show that the new split forms enhance the robustness of high order simulations in comparison to the standard scheme when solving turbulent vortex dominated flows. In fact, we show that for certain test cases, the novel split form discontinuous Galerkin schemes are more robust than the discontinuous Galerkin scheme with over-integration.

    Download full text (pdf)
    Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations
  • 18.
    Ghasemi, Fatemeh
    et al.
    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.
    An energy stable coupling procedure for the compressible and incompressible Navier-Stokes equations2019In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 396, p. 280-302Article in journal (Refereed)
    Abstract [en]

    The coupling of the compressible and incompressible Navier-Stokes equations is considered. Our ambition is to take a first step towards a provably well posed and stable coupling procedure. We study a simplified setting with a stationary planar interface and small disturbances from a steady background flow with zero velocity normal to the interface. The simplified setting motivates the use of the linearized equations, and we derive interface conditions such that the continuous problem satisfy an energy estimate. The interface conditions can be imposed both strongly and weakly. It is shown that the weak and strong interface imposition produce similar continuous energy estimates. We discretize the problem in time and space by employing finite difference operators that satisfy a summation-by-parts rule. The interface and initial conditions are imposed weakly using a penalty formulation. It is shown that the results obtained for the weak interface conditions in the continuous case, lead directly to stability of the fully discrete problem.

    Download full text (pdf)
    fulltext
  • 19.
    Ghasemi, Fatemeh
    et al.
    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.
    On conservation and dual consistency for summation-by-parts based approximations of parabolic problems2020In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 410, article id 109282Article in journal (Refereed)
    Abstract [en]

    n/a

    Download full text (pdf)
    fulltext
  • 20.
    Glaubitz, Jan
    et al.
    Massachusetts Institute of Technology, USA.
    Klein, Simon-Christian
    TU Braunschweig, Germany.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. University of Johannesburg, South Africa.
    Öffner, Philipp
    Johannes Gutenberg University Mainz, Germany.
    Multi-dimensional summation-by-parts operators for general function spaces: Theory and construction2023In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 491, article id 112370Article in journal (Refereed)
    Abstract [en]

    Summation-by-parts (SBP) operators allow us to systematically develop energy-stable and high-order accurate numerical methods for time-dependent differential equations. Until recently, the main idea behind existing SBP operators was that polynomials can accurately approximate the solution, and SBP operators should thus be exact for them. However, polynomials do not provide the best approximation for some problems, with other approximation spaces being more appropriate. We recently addressed this issue and developed a theory for one-dimensional SBP operators based on general function spaces, coined function-space SBP (FSBP) operators. In this paper, we extend the theory of FSBP operators to multiple dimensions. We focus on their existence, connection to quadratures, construction, and mimetic properties. A more exhaustive numerical demonstration of multi-dimensional FSBP (MFSBP) operators and their application will be provided in future works. Similar to the one-dimensional case, we demonstrate that most of the established results for polynomial-based multi-dimensional SBP (MSBP) operators carry over to the more general class of MFSBP operators. Our findings imply that the concept of SBP operators can be applied to a significantly larger class of methods than is currently done. This can increase the accuracy of the numerical solutions and/or provide stability to the methods. © 2023 The Author(s)

    Download full text (pdf)
    fulltext
  • 21.
    Gong, Jing
    et al.
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Nordström, Jan
    Uppsala universitet, Avdelningen för teknisk databehandling.
    A stable and efficient hybrid scheme for viscous problems in complex geometries2007In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 226, p. 1291-1309Article in journal (Refereed)
  • 22.
    Kopriva, David A.
    et al.
    Florida State University, USA; San Diego State University, USA.
    Gassner, Gregor J.
    University of Cologne, Germany.
    Nordström, Jan
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering. University of Johannesburg, South Africa.
    On the theoretical foundation of overset grid methods for hyperbolic problems II: Entropy bounded formulations for nonlinear conservation laws2022In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 471, article id 111620Article in journal (Refereed)
    Abstract [en]

    We derive entropy conserving and entropy dissipative overlapping domain formulations for systems of nonlinear hyperbolic equations in conservation form, such as would be approximated by overset mesh methods. The entropy conserving formulation imposes a two-way coupling at the artificial interface boundaries through nonlinear penalty functions that vanish when the solutions coincide. The penalty functions are expressed in terms of entropy conserving fluxes originally introduced for finite volume schemes. In addition to the interface coupling, which is required, entropy dissipation and coupling can optionally be added through the use of linear penalties within the overlap region.

    Download full text (pdf)
    fulltext
  • 23.
    Kopriva, David A.
    et al.
    The Florida State University, USA; San Diego State University, USA.
    Nordström, Jan
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering. University of Johannesburg, South Africa.
    Gassner, Gregor J.
    University of Cologne, Germany.
    On the theoretical foundation of overset grid methods for hyperbolic problems: Well-posedness and conservation2022In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 448, article id 110732Article in journal (Refereed)
    Abstract [en]

    We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling conditions are used. If a system cannot be diagonalized, as is usually the case in multiple space dimensions, then the energy method does not give proper bounds in terms of initial and boundary data. For those problems, we propose a novel penalty approach. We show, by using a global energy that accounts for the energy in the overlap region of the domains, that under well-defined conditions on the coupling matrices the penalized overset domain problems are energy bounded, conservative, well-posed and have solutions equivalent to the original single domain problem.

    Download full text (pdf)
    fulltext
  • 24.
    Laurén, Fredrik
    et al.
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. Department of Mathematics and Applied Mathematics, University of Johannesburg, Johannesburg, South Africa.
    Energy stable wall modeling for the Navier-Stokes equations2022In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 457, article id 111046Article in journal (Refereed)
    Abstract [en]

    Close to solid walls, at high Reynolds numbers, fluids may develop steep gradients which require a fine mesh for an accurate simulation of the turbulent boundary layer. An often used cure is to use a wall model instead of a fine mesh, with the drawback that modeling is introduced, leading to possibly unstable numerical schemes. In this paper, we leave the modeling aside, take it for granted, and propose a new set of provably energy stable boundary procedures for the incompressible Navier-Stokes equations. We show that these new boundary procedures lead to numerical results with high accuracy even for coarse meshes where data is partially obtained from a wall model.

    Download full text (pdf)
    fulltext
  • 25.
    Laurén, Fredrik
    et al.
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Nordström, Jan
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering. Department of Mathematics and AppliedMathematics, University of Johannesburg, South Africa.
    Spectral properties of the incompressible Navier-Stokes equations2021In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 429, article id 110019Article in journal (Refereed)
    Abstract [en]

    The influence of different boundary conditions on the spectral properties of the incompressible Navier-Stokes equations is investigated. By using the Fourier-Laplace transform technique, we determine the spectra, extract the decay rate in time, and investigate the dispersion relation. In contrast to an infinite domain, where only diffusion affects the convergence, we show that also the propagation speed influence the rate of convergence to steady state for a bounded domain. Once the continuous equations are analyzed, we discretize using high-order finite-difference operators on summation-by-parts form and demonstrate that the continuous analysis carries over to the discrete setting. The theoretical results are verified by numerical experiments, where we highlight the necessity of high accuracy for a correct description of time-dependent phenomena.

    Download full text (pdf)
    fulltext
  • 26.
    Linders, Viktor
    et al.
    Lund University, Sweden.
    Carpenter, Mark H.
    NASA Langley Research Center, USA.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. University of Johannesburg, South Africa.
    A superconvergent stencil-adaptive SBP-SAT finite difference scheme2024In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 501, article id 112794Article in journal (Refereed)
    Abstract [en]

    A stencil-adaptive SBP-SAT finite difference scheme is shown to display superconvergent behavior. As proof of concept, applied to the linear advection equation, it has a convergence rate Ox4) in contrast to a conventional scheme, which converges at a rate Ox3).

    Download full text (pdf)
    fulltext
  • 27.
    Linders, Viktor
    et al.
    Lund University, Lund, Sweden; Technion, Israel Institute of Technology, Haifa, Israel.
    Carpenter, Mark H.
    NASA Langley Research Center, Hampton, United States.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Accurate solution-adaptive finite difference schemes for coarse and fine grids2020In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 410, article id 109393Article in journal (Refereed)
    Abstract [en]

    We introduce solution dependent finite difference stencils whose coefficients adapt to the current numerical solution by minimizing the truncation error in the least squares sense. The resulting scheme has the resolution capacity of dispersion relation preserving difference stencils in under-resolved domains, together with the high order convergence rate of conventional central difference methods in well resolved regions. Numerical experiments reveal that the new stencils outperform their conventional counterparts on all grid resolutions from very coarse to very fine.

    Download full text (pdf)
    fulltext
  • 28.
    Linders, Viktor
    et al.
    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.
    Uniformly Best Wavenumber Approximations by Spatial Central Difference Operators2015In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 300, p. 695-709Article in journal (Refereed)
    Abstract [en]

    We construct accurate central difference stencils for problems involving high frequency waves or multi-frequency solutions over long time intervals with a relatively coarse spatial mesh, and with an easily obtained bound on the dispersion error. This is done by demonstrating that the problem of constructing central difference stencils that have minimal dispersion error in the infinity norm can be recast into a problem of approximating a continuous function from a finite dimensional subspace with a basis forming a Chebyshev set. In this new formulation, characterising and numerically obtaining optimised schemes can be done using established theory.

    Download full text (pdf)
    Uniformly Best Wavenumber Approximations by Spatial Central Difference Operators
  • 29.
    Linders, Viktor
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering. 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. Department of Mathematics and Applied Mathematics, University of Johannesburg, South Africa.
    Frankel, Steven H.
    Faculty of Mechanical Engineering, Technion - Israel Institute of Technology, Haifa, Israel.
    Properties of Runge-Kutta-Summation-By-Parts methods2020In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 419, article id 109684Article in journal (Refereed)
    Abstract [en]

    We review and extend the list of stability and convergence properties satisfied by Runge-Kutta (RK) methods that are associated with Summation-By-Parts (SBP) operators, herein called RK-SBP methods. The analysis covers classical, generalized as well as upwind SBP operators. Previous work on the topic has relied predominantly on energy estimates. In contrast, we derive all results using a purely algebraic approach that rests on the well-established theory of RK methods. The purpose of this paper is to provide a bottom-up overview of stability and convergence results for linear and non-linear problems that relate to general RK-SBP methods. To this end, we focus on the RK viewpoint, since this perspective so far is largely unexplored. This approach allows us to derive all results as simple consequences of the properties of SBP methods combined with well-known results from RK theory. In this way, new proofs of known results such as A-, L- and B-stability are given. Additionally, we establish previously unreported results such as strong S-stability, dissipative stability and stiff accuracy of certain RK-SBP methods. Further, it is shown that a subset of methods are B-convergent for strictly contractive non-linear problems and convergent for non-linear problems that are both contractive and dissipative.

    Download full text (pdf)
    fulltext
  • 30.
    Linders, Viktor
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Kupiainen, Marco
    Rossby Centre SMHI, Norrköping, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Summation-by-Parts Operators with Minimal Dispersion Error for Coarse Grid Flow Calculations2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 340, p. 34p. 160-176Article in journal (Refereed)
    Abstract [en]

    We present a procedure for constructing Summation-by-Parts operators with minimal dispersion error both near and far from numerical interfaces. Examples of such operators are constructed and compared with a higher order non-optimised Summation-by-Parts operator. Experiments show that the optimised operators are superior for wave propagation and turbulent flows involving large wavenumbers, long solution times and large ranges of resolution scales.

    Download full text (pdf)
    fulltext
  • 31.
    Lindström, Jens
    et al.
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Nordström, Jan
    Uppsala universitet, Avdelningen för teknisk databehandling.
    A stable and high-order accurate conjugate heat transfer problem2010In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 229, no 14, p. 5440-5456Article in journal (Refereed)
    Abstract [en]

    This paper analyzes well-posedness and stability of a conjugate heat transfer problem in one space dimension. We study a model problem for heat transfer between a fluid and a solid. The energy method is used to derive boundary and interface conditions that make the continuous problem well-posed and the semi-discrete problem stable. The numerical scheme is implemented using 2nd-, 3rd- and 4th-order finite difference operators on Summation-By-Parts (SBP) form. The boundary and interface conditions are implemented weakly. We investigate the spectrum of the spatial discretization to determine which type of coupling that gives attractive convergence properties. The rate of convergence is verified using the method of manufactured solutions.

    Download full text (pdf)
    fulltext
  • 32.
    Lundquist, Tomas
    et al.
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Laurén, Fredrik
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. University of Johannesburg, South Africa.
    A multi-domain summation-by-parts formulation for complex geometries2022In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 463, article id 111269Article in journal (Refereed)
    Abstract [en]

    We combine existing summation-by-parts discretization methods to obtain a simplified numerical framework for partial differential equations posed on complex multi-block/element domains. The interfaces (conforming or non-conforming) between blocks are treated with inner-product-preserving interpolation operators, and the result is a high-order multi-block operator on summation-by-parts form that encapsulates both the metric terms as well as the interface treatments. This enables for a compact description of the numerical scheme that mimics the essential features of its continuous counterpart. Furthermore, the stability analysis on a multi-block domain is simplified for both for linear and nonlinear equations, since no problem-specific interface conditions need to be derived and implemented. We exemplify the combined operator technique by considering a nonlinearly stable discrete formulation of the incompressible Navier-Stokes equations and perform calculations on an underlying multi-block domain. 

    Download full text (pdf)
    fulltext
  • 33.
    Lundquist, Tomas
    et al.
    Industrial CFD Research Group, Department of Mechanical Engineering, University of Cape Town, South Africa.
    Malan, Arnaud G.
    Industrial CFD Research Group, Department of Mechanical Engineering, University of Cape Town, South Africa.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Efficient and error minimized coupling procedures for unstructured and moving meshes2020In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 406, p. 1-21, article id 109158Article in journal (Refereed)
    Abstract [en]

    We present a methodology for automatic generation and optimization of interpolation operators for the coupling of general non-collocated and/or moving numerical interfaces. The discrete equations are solved in a method-of-lines fashion by assuming volume preserving mesh motions. Interface interpolation errors are minimized effectively in a global least-squares sense, while satisfying strict stability conditions. The proposed automatic interface procedure is both more versatile and more accurate compared to previous techniques. We apply the new method to interfaces between hybrid meshes undergoing relative rigid body motion, demonstrating the stability, conservation and superior accuracy.

    Download full text (pdf)
    Efficient and error minimized coupling procedures for unstructured and moving meshes
  • 34.
    Lundquist, Tomas
    et al.
    Industrial CFD Research Group, Department of Mechanical Engineering, University of Cape Town, South Africa.
    Malan, Arnaud
    Industrial CFD Research Group, Department of Mechanical Engineering, University of Cape Town, South Africa.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    A hybrid framework for coupling arbitrary summation-by-parts schemes on general meshes2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 362, p. 49-68Article in journal (Refereed)
    Abstract [en]

    We develop a general interface procedure to couple both structured and unstructured parts of a hybrid mesh in a non-collocated, multi-block fashion. The target is to gain optimal computational efficiency in fluid dynamics simulations involving complex geometries. While guaranteeing stability, the proposed procedure is optimized for accuracy and requires minimal algorithmic modifications to already existing schemes. Initial numerical investigations confirm considerable efficiency gains compared to non-hybrid calculations of up to an order of magnitude.

    Download full text (pdf)
    fulltext
  • 35.
    Lundquist, Tomas
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    The SBP-SAT Technique for Initial Value Problems2014In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 270, p. 86-104Article in journal (Refereed)
    Abstract [en]

    A detailed account of the stability and accuracy properties of the SBP-SAT technique for numerical time integration is presented. We show how the technique can be used to formulate both global and multi-stage methods with high order of accuracy for both stiff and non-stiff problems. Linear and non- linear stability results, including A-stability, L-stability and B-stability are proven using the energy method for general initial value problems. Numerical experiments corroborate the theoretical properties.

    Download full text (pdf)
    The SBP-SAT Technique for Initial Value Problems
  • 36.
    Lundquist, Tomas
    et al.
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering.
    Winters, Andrew Ross
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. University of Johannesburg, South Africa.
    Encapsulated generalized summation-by-parts formulations for curvilinear and non-conforming meshes2024In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 498, article id 112699Article in journal (Refereed)
    Abstract [en]

    We extend the construction of so-called encapsulated global summation-by-parts operators to the general case of a mesh which is not boundary conforming. Owing to this development, energy stable discretizations of nonlinear and variable coefficient initial boundary value problems can be formulated in simple and straightforward ways using high-order accurate operators of generalized summation-by-parts type. Encapsulated features on a single computational block or element may include polynomial bases, tensor products as well as curvilinear coordinate transformations. Moreover, through the use of inner product preserving interpolation or projection, the global summation-by-parts property is extended to arbitrary multi-block or multi-element meshes with non-conforming nodal interfaces.

    Download full text (pdf)
    fulltext
  • 37.
    Mattsson, Ken
    et al.
    Department of Information Technology, Uppsala University, PO Box 337, SE-75105 Uppsala, Sweden.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Corrigendum to “Summation by parts operators for finite difference approximations of second derivatives” [J.Comput.Phys.199 (2004) 503–540]2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 351, p. 535-Article in journal (Other academic)
    Download full text (pdf)
    fulltext
  • 38. Mattsson, Ken
    et al.
    Nordström, Jan
    Uppsala universitet, Avdelningen för teknisk databehandling.
    High order finite difference methods for wave propagation in discontinuous media2006In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 220, p. 249-269Article in journal (Refereed)
  • 39.
    Mattsson, Ken
    et al.
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Nordström, Jan
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Summation by parts operators for finite difference approximations of second derivatives2004In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 199, p. 503-540Article in journal (Refereed)
  • 40.
    Nikkar, Samira
    et al.
    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 dual consistent summation-by-parts formulation for the linearized incompressible Navier-Stokes equations posed on deforming domains2019In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 376, p. 26p. 322-338Article in journal (Refereed)
    Abstract [en]

    In this article, well-posedness and dual consistency of the linearized constant coefficient incompressible Navier–Stokes equations posed on time-dependent spatial domains are studied. To simplify the derivation of the dual problem and improve the accuracy of gradients, the second order formulation is transformed to first order form. Boundary conditions that simultaneously lead to boundedness of the primal and dual problems are derived.Fully discrete finite difference schemes on summation-by-parts form, in combination with the simultaneous approximation technique, are constructed. We prove energy stability and discrete dual consistency and show how to construct the penalty operators such that the scheme automatically adjusts to the variations of the spatial domain. As a result of the aforementioned formulations, stability and discrete dual consistency follow simultaneously.The method is illustrated by considering a deforming time-dependent spatial domain in two dimensions. The numerical calculations are performed using high order operators in space and time. The results corroborate the stability of the scheme and the accuracy of the solution. We also show that linear functionals are superconverging. Additionally, we investigate the convergence of non-linear functionals and the divergence of the solution.

    Download full text (pdf)
    fulltext
  • 41.
    Nikkar, Samira
    et al.
    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 fully discrete, stable and conservative summation-by-parts formulation for deforming interfaces2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 339, p. 25p. 500-524Article in journal (Refereed)
    Abstract [en]

    We introduce an interface/coupling procedure for hyperbolic problems posed on time-dependent curved multi-domains. First, we transform the problem from Cartesian to boundary-conforming curvilinear coordinates and apply the energy method to derive well-posed and conservative interface conditions. Next, we discretize the problem in space and time by employing finite difference operators that satisfy a summation-by-parts rule. The interface condition is imposed weakly using a penalty formulation. We show how to formulate the penalty operators such that the coupling procedure is automatically adjusted to the movements and deformations of the interface, while both stability and conservation conditions are respected. The developed techniques are illustrated by performing numerical experiments on the linearized Euler equations and the Maxwell equations. The results corroborate the stability and accuracy of the fully discrete approximations.

    Download full text (pdf)
    fulltext
  • 42.
    Nikkar, Samira
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Nordström, Jan
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
    Fully discrete energy stable high order finite difference methods for hyperbolic problems in deforming domains2015In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 291, p. 82-98Article in journal (Refereed)
    Abstract [en]

    A time-dependent coordinate transformation of a constant coefficient hyperbolic system of equations which results in a variable coefficient system of equations is considered. By applying the energy method, well-posed boundary conditions for the continuous problem are derived. Summation-by-Parts (SBP) operators for the space and time discretization, together with a weak imposition of boundary and initial conditions using Simultaneously Approximation Terms (SATs) lead to a provable fully-discrete energy-stable conservative finite difference scheme. We show how to construct a time-dependent SAT formulation that automatically imposes boundary conditions, when and where they are required. We also prove that a uniform flow field is preserved, i.e. the Numerical Geometric Conservation Law (NGCL) holds automatically by using SBP-SAT in time and space. The developed technique is illustrated by considering an application using the linearized Euler equations: the sound generated by moving boundaries. Numerical calculations corroborate the stability and accuracy of the new fully discrete approximations.

    Download full text (pdf)
    fulltext
  • 43.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. University of Johannesburg, South Africa.
    A skew-symmetric energy and entropy stable formulation of the compressible Euler equations2022In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 470, article id 111573Article in journal (Refereed)
    Abstract [en]

    We show that a specific skew-symmetric form of nonlinear hyperbolic problems leads to energy and entropy bounds. Next, we exemplify by considering the compressible Euler equations in primitive variables, transform them to skew-symmetric form and show how to obtain energy and entropy estimates. Finally we show that the skew-symmetric formulation lead to energy and entropy stable discrete approximations if the scheme is formulated on summation-by-parts form.

    Download full text (pdf)
    fulltext
  • 44.
    Nordström, Jan
    Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering. Department of Mathematics and Applied Mathematics, University of Johannesburg, South Africa.
    Nonlinear and linearised primal and dual initial boundary value problems: When are they bounded? How are they connected?2022In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 455, article id 111001Article in journal (Refereed)
    Abstract [en]

    Linearisation is often used as a first step in the analysis of nonlinear initial boundary value problems. The linearisation procedure frequently results in a confusing contradiction where the nonlinear problem conserves energy and has an energy bound but the linearised version does not (or vice versa). In this paper we attempt to resolve that contradiction and relate nonlinear energy conserving and bounded initial boundary value problems to their linearised versions and the related dual problems. We start by showing that a specific skew-symmetric form of the primal nonlinear problem leads to energy conservation and a bound. Next, we show that this specific form together with a non-standard linearisation procedure preserves these properties for the new slightly modified linearised problem. We proceed to show that the corresponding linear and nonlinear dual (or self-adjoint) problems also have bounds and conserve energy due to this specific formulation. Next, the implication of the new formulation on the choice of boundary conditions is discussed. A straightforward nonlinear and linear analysis may lead to a different number and type of boundary conditions required for an energy bound. We show that the new formulation sheds somelight on this contradiction. We conclude by illustrating that the new continuous formulation automatically leads toenergy stable and energy conserving numerical approximations for both linear and nonlinear primal and dual problems if the approximations are formulated on summation-by-parts form.

    Download full text (pdf)
    fulltext
  • 45.
    Nordström, Jan
    Linköping University, Faculty of Science & Engineering. Linköping University, Department of Mathematics, Applied Mathematics. University of Johannesburg, South Africa.
    Nonlinear Boundary Conditions for Initial Boundary Value Problems with Applications in Computational Fluid Dynamics2024In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 498, article id 112685Article in journal (Refereed)
    Abstract [en]

    We derive new boundary conditions and implementation procedures for nonlinear initial boundary value problems (IBVPs) with non-zero boundary data that lead to bounded solutions. The new boundary procedure is applied to nonlinear IBVPs in skew-symmetric form, including dissipative terms. The complete procedure has two main ingredients. Firstly, the energy rate in terms of a surface integral with boundary terms is derived. Secondly, we bound the surface integral by deriving new nonlinear boundary procedures for boundary conditions with non-zero data. The new nonlinear boundary procedure generalises the well known characteristic boundary procedure for linear problems to the nonlinear setting.

    To introduce the procedure, a skew-symmetric scalar IBVP encompassing the linear advection equation and Burgers equation is analysed. Once the continuous analysis is done, we show that energy stable nonlinear discrete approximations follow by using summation-by-parts operators combined with weak boundary conditions. The scalar analysis is subsequently repeated for general nonlinear systems of equations. Finally, the new boundary procedure is applied to four important IBVPs in computational fluid dynamics: the incompressible Euler and Navier-Stokes, the shallow water and the compressible Euler equations.

    Download full text (pdf)
    fulltext
  • 46.
    Nordström, Jan
    et al.
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Carpenter, Mark H.
    Boundary and interface conditions for high-order finite-difference methods applied to the Euler and Navier-Stokes equations1999In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 148, p. 621-645Article in journal (Refereed)
  • 47.
    Nordström, Jan
    et al.
    Uppsala universitet, Avdelningen för teknisk databehandling.
    Carpenter, Mark H.
    High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates2001In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 173, p. 149-174Article in journal (Refereed)
  • 48.
    Nordström, Jan
    et al.
    Linköping University, Department of Mathematics, Scientific Computing. Linköping University, The Institute of Technology.
    Eriksson, Sofia
    Uppsala University.
    Eliasson, Peter
    Swedish Def Research Agency.
    Weak and strong wall boundary procedures and convergence to steady-state of the Navier-Stokes equations2012In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 231, no 14, p. 4867-4884Article in journal (Refereed)
    Abstract [en]

    We study the influence of different implementations of no-slip solid wall boundary conditions on the convergence to steady-state of the Navier-Stokes equations. The various approaches are investigated using the energy method and an eigenvalue analysis. It is shown that the weak implementation is superior and enhances the convergence to steady-state for coarse meshes. It is also demonstrated that all the stable approaches produce the same convergence rate as the mesh size goes to zero. The numerical results obtained by using a fully nonlinear finite volume solver support the theoretical findings from the linear analysis. 

    Download full text (pdf)
    fulltext
  • 49.
    Nordström, Jan
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Ghasemi, Fatemeh
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Corrigendum to “On the relation between conservation and dual consistency for summation-by-parts schemes”[J. Comput. Phys. 344 (2017) 437–439]2018In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 360, p. 247-247Article in journal (Other academic)
    Download full text (pdf)
    fulltext
  • 50.
    Nordström, Jan
    et al.
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    Ghasemi, Fatemeh
    Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
    On the relation between conservation and dual consistency for summation-by-parts schemes2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 344, p. 3p. 437-439Article in journal (Refereed)
    Abstract [en]

    n/a

    Download full text (pdf)
    fulltext
12 1 - 50 of 92
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