Search

We derive boundary conditions and estimates based on the energy and entropy analysis of systems of the nonlinear shallow water equations in two spatial dimensions. It is shown that the energy method provides more details, but is fully consistent with the entropy analysis. The details brought forward by the nonlinear energy analysis allow us to pinpoint where the difference between the linear and nonlinear analysis originate. We find that the result from the linear analysis does not necessarily hold in the nonlinear case. The nonlinear analysis leads in general to a different minimal number of boundary conditions compared with the linear analysis. In particular, and contrary to the linear case, the magnitude of the flow does not influence the number of required boundary conditions. Funding: Vetenskapsradet, Sweden [2018-05084 VR, 2020-03642 VR]

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. Funding: Vetenskapsradet, Sweden [2020-03642, 2018-05084, 2021-05484]; Swedish e-Science Research Center (SeRC)

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. Funding agencies: Vetenskapsradet, Sweden [2018-05084 VR]; Swedish e-Science Research Center (SeRC)

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. Funding: Simons Foundation [426393]; Vetenskapsradet, SwedenSwedish Research Council [2018-05084 VR]; Swedish e-Science Research Center (SeRC); Klaus-Tschira Stiftung; European Research CouncilEuropean Research Council (ERC)European Commission [71448]

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. Funding: Vetenskapsradet, SwedenSwedish Research Council [2018-05084 VR]; Swedish e-Science Research Center (SeRC)

We extend the convergence results in Svärd and Nordström (2019) [7] for single-domain energy-stable high-order finite difference schemes, to include domains split into several grid blocks. The analysis also demonstrates that reflective boundary conditions enjoy the same convergence properties. Finally, we briefly indicate that these results (and the previous ones in [7]) also hold in multiple dimensions. publishedVersion

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. Funding agency: The Swedish e-Science Research Centre (SeRC)

We develop a novel numerical scheme for the simulation of dissipative quantum dynamics following from two-body Lindblad master equations. All defining continuum properties of the Lindblad dynamics, hermiticity, positivity and in particular trace conservation of the evolved density matrix are preserved. The central ingredient is a new spatial difference operator, which not only fulfils the summation by parts (SBP) property but also implements a continuum reparametrization property. Using the time evolution of a heavy-quark anti-quark bound state in a hot thermal medium as an explicit example, we show how the reparametrization neutral summation-by-parts (RN-SBP) operator preserves the continuum properties of the theory., 34 pages, 7 figures, open-access code available via https://doi.org/10.5281/zenodo.3744460

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. Funding: Vetenskapsradet, SwedenSwedish Research Council [2018-05084 VR]; Swedish e-Science Research Center (SeRC)

We study an incompletely parabolic system in three space dimensions with stochastic boundary and initial data. We show how the variance of the solution can be manipulated by the boundary conditions, while keeping the mean value of the solution unaffected. Estimates of the variance of the solution is presented both analytically and numerically. We exemplify the technique by applying it to an incompletely parabolic model problem, as well as the one-dimensional compressible Navier–Stokes equations. Funding agencies: European Commission [ACP3-GA-2013-605036]

We introduce the notion of a transmission problem to describe a general class of problems where different dynamics are coupled in time. Well-posedness and stability are analysed for continuous and discrete problems using both strong and weak formulations, and a general transmission condition is obtained. The theory is applied to the coupling of fluid-acoustic models, multi-grid implementations, adaptive mesh refinements, multi-block formulations and numerical filtering. Funding agencies: Swedish Meteorological and Hydrological Institute (SMHI)

Multigrid schemes for high order finite difference methods on summation-by-parts form are studied by comparing the effect of different interpolation operators. By using the standard linear prolongation and restriction operators, the Galerkin condition leads to inaccurate coarse grid discretizations. In this paper, an alternative class of interpolation operators that bypass this issue and preserve the summation-by-parts property on each grid level is considered. Clear improvements of the convergence rate for relevant model problems are achieved. Funding agencies: VINNOVA, the Swedish Governmental Agency for Innovation Systems [2013-01209]

We investigate the spatial operator in the incompressible Navier–Stokes, Oseen and Stokes equations and show how to avoid singularities associated with null spaces by choosing specific boundary conditions. The theoretical results are derived for a general form of energy stable boundary conditions, and applied to a few commonly used ones. The analysis is done on a system that simultaneously covers the nonlinear incompressible Navier–Stokes, the Oseen and the Stokes equations. When the spectrum of the spatial operator is investigated, we restrict the analysis to the Oseen and Stokes equations. The continuous analysis carries over to the discrete setting by using the summation-by-parts framework.

We consider the coupling of parabolic problems discretized using difference operators on summation-by-parts (SBP) form with interface conditions imposed weakly. In [1, 2], it was shown that conserv ...

In this paper we study boundary conditions for linear hyperbolic systems of equations and the corresponding dual problem. In particular, we show that the primal and dual boundary conditions are related by a simple scaling relation. It is also shown that the weak dual problem can be derived directly from the weak primal problem. Based on the continuous analysis, we discretize and perform computations with a high-order finite difference scheme on summation- by-parts form with weak boundary conditions. It is shown that the results obtained in the continuous analysis lead directly to stability results for the primal and dual discrete problems. Numerical experiments corroborate the theoretical results.

The effect of an inaccurate geometry description on the solution accuracy of a hyperbolic problem is discussed. The inaccurate geometry can for example come from an imperfect CAD system, a faulty mesh generator, bad measurements or simply a misconception. We show that inaccurate geometry descriptions might lead to the wrong wave speeds, a misplacement of the boundary conditions, to the wrong boundary operator and a mismatch of boundary data. The errors caused by an inaccurate geometry description may affect the solution more than the accuracy of the specific discretization techniques used. In extreme cases, the order of accuracy goes to zero. Numerical experiments corroborate the theoretical results.

Unresolved gradients produce numerical oscillations and inaccurate results. The most straightforward solution to such a problem is to increase the resolution of the computational grid. However, this is often prohibitively expensive and may lead to ecessive execution times. By training a neural network to predict the shape of the solution, we show that it is possible to reduce numerical oscillations and increase both accuracy and efficiency. Data from the neural network prediction is imposed using multiple penalty terms inside the domain.

Artificial neural networks together with associated computational libraries provide a powerful framework for constructing both classification and regression algorithms. In this paper we use neural networks to design linear and non-linear discrete differential operators. We show that neural network based operators can be used to construct stable discretizations of initial boundary-value problems by ensuring that the operators satisfy a discrete analogue of integration-by-parts known as summation-by-parts. Our neural network approach with linear activation functions is compared and contrasted with a more traditional linear algebra approach. An application to overlapping grids is explored. The strategy developed in this work opens the door for constructing stable differential operators on general meshes.

n/a Classified in the journal as "Short note"

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.

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.

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.

Stochastic problems governed by nonlinear conservation laws are challenging due to solution discontinuities in stochastic and physical space. In this paper, we present a level set method to track discontinuities in stochastic space by solving a Hamilton-Jacobi equation. By introducing a speed function that vanishes at discontinuities, the iso-zeros of the level set problem coincide with the discontinuities of the conservation law. The level set problem is solved on a sequence of successively finer grids in stochastic space. The method is adaptive in the sense that costly evaluations of the conservation law of interest are only performed in the vicinity of the discontinuities during the refinement stage. In regions of stochastic space where the solution is smooth, a surrogate method replaces expensive evaluations of the conservation law. The proposed method is tested in conjunction with different sets of localized orthogonal basis functions on simplex elements, as well as frames based on piecewise polynomials conforming to the level set function. The performance of the proposed method is compared to existing adaptive multi-element generalized polynomial chaos methods. Funding agencies: Research Council of Norway through CLIMIT program [244035/E20 CONQUER]; U.S. Department of Energy Office of Science, Office of Advanced Scientific Computing Research [DE-SC0006402]; NSF [CMMI-1454601]

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.

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.

This work is devoted to the application of the super compact finite difference method (SCFDM) and the combined compact finite difference method (CCFDM) for spatial differencing of the spherical sha ...

We consider a hyperbolic system with uncertainty in the boundary and initial data. Our aim is to show that different boundary conditions gives different convergence rates of the variance of the solution. This means that we can with the same knowledge of data get a more or less accurate description of the uncertainty in the solution. A variety of boundary conditions are compared and both analytical and numerical estimates of the variance of the solution is presented. As applications, we study the effect of this technique on Maxwell's equations as well as on a subsonic outflow boundary for the Euler equations.

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.

Corrigendum to “On the relation between conservation and dual consistency for summation-by-parts schemes”[J. Comput. Phys. 344 (2017) 437–439]

Fully-implicit discrete formulations in summation-by-parts form for initial-boundary value problems must be invertible in order to provide well functioning procedures. We prove that, under mild assumptions, pseudo-spectral collocation methods for the time derivative lead to invertible discrete systems when energy-stable spatial discretizations are used.

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.

In this paper we derive new far-field 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-posedness 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.

The Euler equations subject to uncertainty in the initial and boundary conditions are investigated via the stochastic Galerkin approach. We present a new fully intrusive method based on a variable transformation of the continuous equations. Roe variables are employed to get quadratic dependence in the flux function and a well-defined Roe average matrix that can be determined without matrix inversion. In previous formulations based on generalized polynomial chaos expansion of the physical variables, the need to introduce stochastic expansions of inverse quantities, or square-roots of stochastic quantities of interest, adds to the number of possible different ways to approximate the original stochastic problem. We present a method where the square roots occur in the choice of variables, resulting in an unambiguous problem formulation. The Roe formulation saves computational cost compared to the formulation based on expansion of conservative variables. Moreover, the Roe formulation is more robust and can handle cases of supersonic flow, for which the conservative variable formulation fails to produce a bounded solution. For certain stochastic basis functions, the proposed method can be made more effective and well-conditioned. This leads to increased robustness for both choices of variables. We use a multi-wavelet basis that can be chosen to include a large number of resolution levels to handle more extreme cases (e.g. strong discontinuities) in a robust way. For smooth cases, the order of the polynomial representation can be increased for increased accuracy.

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.

The stochastic Galerkin and collocation methods are used to solve an advection–diffusion equation with uncertain and spatially varying viscosity. We investigate well-posedness, monotonicity and stability for the extended system resulting from the Galerkin projection of the advection–diffusion equation onto the stochastic basis functions. High-order summation-by-parts operators and weak imposition of boundary conditions are used to prove stability of the semi-discrete system. It is essential that the eigenvalues of the resulting viscosity matrix of the stochastic Galerkin system are positive and we investigate conditions for this to hold. When the viscosity matrix is diagonalizable, stochastic Galerkin and stochastic collocation are similar in terms of computational cost, and for some cases the accuracy is higher for stochastic Galerkin provided that monotonicity requirements are met. We also investigate the total spatial operator of the semi-discretized system and its impact on the convergence to steady-state.

Corrigendum to “Summation by parts operators for finite difference approximations of second derivatives” [J.Comput.Phys.199 (2004) 503–540]

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.

For wave propagation over distances of many wavelengths, high-order finite difference methods on staggered grids are widely used due to their excellent dispersion properties. However, the enforcement of boundary conditions in a stable manner and treatment of interface problems with discontinuous coefficients usually pose many challenges. In this work, we construct a provably stable and high-order-accurate finite difference method on staggered grids that can be applied to a broad class of boundary and interface problems. The staggered grid difference operators are in summation-by-parts form and when combined with a weak enforcement of the boundary conditions, lead to an energy stable method on multiblock grids. The general applicability of the method is demonstrated by simulating an explosive acoustic source, generating waves reflecting against a free surface and material discontinuity. Funding agencies: Department of Geophysics at Stanford University

We discuss conservative and stable numerical approximations in summation-by-parts form for linear hyperbolic problems with variable coefficients. An extended setting, where the boundary or interface may or may not be included in the grid, is considered. We prove that conservative and stable formulations for variable coefficient problems require a boundary and interface conforming grid and exact numerical mimicking of integration-by-parts. Finally, we comment on how the conclusions from the linear analysis carry over to the nonlinear setting. Funding agencies: VINNOVA [2013-01209]

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. Original Publication:Jens Berg and Jan Nordström, Stable Robin solid wall boundary conditions for the Navier-Stokes equations, 2011, Journal of Computational Physics, (230), 19, 7519-7532.http://dx.doi.org/10.1016/j.jcp.2011.06.027Copyright: Elsevierhttp://www.elsevier.com/

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. Original Publication:Jens Lindström and Jan Nordström, A stable and high-order accurate conjugate heat transfer problem, 2010, Journal of Computational Physics, (229), 5440-5456.http://dx.doi.org/10.1016/j.jcp.2010.04.010Copyright: Elsevier Science B.V., Amsterdamhttp://www.elsevier.com/

A stable wall boundary procedure is derived for the discretized compressible Navier-Stokes equations. The procedure leads to an energy estimate for the linearized equations. We discretize the equations using high-order accurate finite difference summation-by-parts (SBP) operators. The boundary conditions are imposed weakly with penalty terms. We prove linear stability for the scheme including the wall boundary conditions. The penalty imposition of the boundary conditions is tested for the flow around a circular cylinder at Ma=0.1 and Re=100. We demonstrate the robustness of the SBP-SAT technique by imposing incompatible initial data and show the behavior of the boundary condition implementation. Using the errors at the wall we show that higher convergence rates are obtained for the high-order schemes. We compute the vortex shedding from a circular cylinder and obtain good agreement with previously published (computational and experimental) results for lift, drag and the Strouhal number. We use our results to compare the computational time for a given for a accuracy and show the superior efficiency of the 5th-order scheme.

We present a well-posed stochastic Galerkin formulation of the incompressible Navier–Stokes equations with uncertainty in model parameters or the initial and boundary conditions. The stochastic Galerkin method involves representation of the solution through generalized polynomial chaos expansion and projection of the governing equations onto stochastic basis functions, resulting in an extended system of equations. A relatively low-order generalized polynomial chaos expansion is sufficient to capture the stochastic solution for the problem considered. We derive boundary conditions for the continuous form of the stochastic Galerkin formulation of the velocity and pressure equations. The resulting problem formulation leads to an energy estimate for the divergence. With suitable boundary data on the pressure and velocity, the energy estimate implies zero divergence of the velocity field. Based on the analysis of the continuous equations, we present a semi-discretized system where the spatial derivatives are approximated using finite difference operators with a summation-by-parts property. With a suitable choice of dissipative boundary conditions imposed weakly through penalty terms, the semi-discrete scheme is shown to be stable. Numerical experiments in the laminar flow regime corroborate the theoretical results and we obtain high-order accurate results for the solution variables and the velocity divergence converges to zero as the mesh is refined. Funding agencies: SUPRI-B at Stanford University; U.S. Department of Energy Office of Science, Office of Advanced Scientific Computing Research [DE-SC0006402]; Uni Research, Norway

In this paper, we present a stable hybrid scheme for viscous problems. The hybrid method combines the unstructured finite volume method with high-order finite difference methods on complex geometries. The coupling procedure between the two numerical methods is based on energy estimates and stable interface conditions are constructed. Numerical calculations show that the hybrid method is efficient and accurate.

We construct a stable high-order finite difference scheme for the compressible Navier-Stokes equations, that satisfy an energy estimate. The equations are discretized with high-order accurate finite difference methods that satisfy a Summation-By-Parts rule. The boundary conditions are imposed with penalty terms known as the Simultaneous Approximation Term technique. The main result is a stability proof for the full three-dimensional Navier-Stokes equations, including the boundary conditions. We show the theoretical third-, fourth-, and fifth-order convergence rate, for a viscous shock, where the analytic solution is known. We demonstrate the stability and discuss the non-reflecting properties of the outflow conditions for a vortex in free space. Furthermore, we compute the three-dimensional vortex shedding behind a circular cylinder in an oblique free stream for Mach number 0.5 and Reynolds number 500.

New sets of boundary conditions for the velocity-pressure formulation of the incompressible Navier-Stokes equations are derived. The boundary conditions have the same form on both inflow and outflow boundaries and lead to a divergence free solution. Moreover, the specific form of the boundary condition makes it possible derive a symmetric positive definite equation system for the internal pressure. Numerical experiments support the theoretical conclusions.

High order finite difference approximations are derived for the second order wave equation with discontinuous coefficients, on rectangular geometries. The discontinuity is treated by splitting the domain at the discontinuities in a multi block fashion. Each sub-domain is discretized with compact second derivative summation by parts operators and the blocks are patched together to a global domain using the projection method. This guarantees a conservative, strictly stable and high order accurate scheme. The analysis is verified by numerical simulations in one and two spatial dimensions.

Finite difference approximations of the second derivative in space appearing in, parabolic, incompletely parabolic systems of, and 2nd-order hyperbolic, partial differential equations are considered. If the solution is pointwise bounded, we prove that finite difference approximations of those classes of equations can be closed with two orders less accuracy at the boundary without reducing the global order of accuracy.This result is generalised to initial-boundary value problems with an mth-order principal part. Then, the boundary accuracy can be lowered m orders.Further, it is shown that schemes using summation-by-parts operators that approximate second derivatives are pointwise bounded. Linear and nonlinear computations, including the two-dimensional Navier-Stokes equations, corroborate the theoretical results.

A stable hybrid method for hyperbolic problems that combines the unstructured finite volume method with high-order finite difference methods has been developed. The coupling procedure is based on energy estimates and stability can be guaranteed. Numerical calculations verify that the hybrid method is efficient and accurate.