Change search
ReferencesLink to record
Permanent link

Direct link
A fully discrete, stable and conservative summation-by-parts formulation for deforming interfaces
Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, Faculty of Science & Engineering.
2016 (English)Report (Other academic)
Abstract [en]

We introduce an interface/coupling procedure for hyperbolic problems posedon 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.

Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2016. , 38 p.
Series
LiTH-MAT-R, ISSN 0348-2960 ; 2016:9
Keyword [en]
Finite difference, High order accuracy, Deforming domains, Time-dependent interface, Well-posedness, Conservation, Summation-by-parts, Stability, Hyperbolic problems
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:liu:diva-130583ISRN: LiTH-MAT-R--2016/09--SEOAI: oai:DiVA.org:liu-130583DiVA: diva2:953284
Available from: 2016-08-17 Created: 2016-08-17 Last updated: 2016-09-27Bibliographically approved
In thesis
1. Stable High Order Finite Difference Methods for Wave Propagation and Flow Problems on Deforming Domains
Open this publication in new window or tab >>Stable High Order Finite Difference Methods for Wave Propagation and Flow Problems on Deforming Domains
2016 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

We construct stable, accurate and efficient numerical schemes for wave propagation and flow problems posed on spatial geometries that are moving, deforming, erroneously described or non-simply connected. The schemes are on Summation-by-Parts (SBP) form, combined with the Simultaneous Approximation Term (SAT) technique for imposing initial and boundary conditions. The main analytical tool is the energy method, by which well-posedness, stability and conservation are investigated. To handle the deforming domains, time-dependent coordinate transformations are used to map the problem to fixed geometries.

The discretization is performed in such a way that the Numerical Geometric Conservation Law (NGCL) is satisfied. Additionally, even though the schemes are constructed on fixed domains, time-dependent penalty formulations are necessary, due to the originally moving boundaries. We show how to satisfy the NGCL and present an automatic formulation for the penalty operators, such that the correct number of boundary conditions are imposed, when and where required.

For problems posed on erroneously described geometries, we investigate how the accuracy of the solution is affected. It is shown that the inaccurate geometry descriptions may lead to wrong wave speeds, a misplacement of the boundary condition, the wrong boundary operator or a mismatch of data. Next, the SBP-SAT technique is extended to time-dependent coupling procedures for deforming interfaces in hyperbolic problems. We prove conservation and stability and show how to formulate the penalty operators such that the coupling procedure is automatically adjusted to the variations of the interface location while the NGCL is preserved.

Moreover, dual consistent SBP-SAT schemes for the linearized incompressible Navier-Stokes equations posed on deforming domains are investigated. To simplify the derivations of the dual problem and incorporate the motions of the boundaries, the second order formulation is reduced to first order and the problem is transformed to a fixed domain. We prove energy stability and dual consistency. It is shown that the solution as well as the divergence of the solution converge with the design order of accuracy, and that functionals of the solution are superconverging.

Finally, initial boundary value problems posed on non-simply connected spatial domains are investigated. The new formulation increases the accuracy of the scheme by minimizing the use of multi-block couplings. In order to show stability, the spectrum of the semi-discrete SBP-SAT formulation is studied. We show that the eigenvalues have the correct sign, which implies stability, in combination with the SBP-SAT technique in time.

Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2016. 23 p.
Series
Linköping Studies in Science and Technology. Dissertations, ISSN 0345-7524 ; 1774
National Category
Computational Mathematics Mathematical Analysis Fluid Mechanics and Acoustics
Identifiers
urn:nbn:se:liu:diva-130928 (URN)10.3384/diss.diva-130928 (DOI)9789176857373 (Print) (ISBN)
Public defence
2016-09-30, Visionen, Hus B, Campus Valla, Linköping, 13:15 (English)
Opponent
Supervisors
Available from: 2016-08-31 Created: 2016-08-31 Last updated: 2016-09-01Bibliographically approved

Open Access in DiVA

Version 2. Fully discrete, stable and conservative summation-by-parts formulation for deforming interfaces(2446 kB)0 downloads
File information
File name FULLTEXT03.pdfFile size 2446 kBChecksum SHA-512
eb4032cfb8df2bde01b310781aaeb674d1be3738321434e44cac2614ff12c3b5003ad8ad98c9b5226106facda775e48125714b60b5912acf7735e1429c588750
Type fulltextMimetype application/pdf
Version 1. Fully discrete, stable and conservative summation-by-parts formulation for deforming interfaces(1985 kB)10 downloads
File information
File name FULLTEXT02.pdfFile size 1985 kBChecksum SHA-512
07ab734935cbd5e9fa297f630d170fd3df54dab40b5ecf3b1319a3f8a9cbd9a36b2cfe2930b7fa9b34e180a022f8b29ed00df183558c3ee2efc0b16883c886a1
Type fulltextMimetype application/pdf

Search in DiVA

By author/editor
Nikkar, SamiraNordström, Jan
By organisation
Computational MathematicsFaculty of Science & Engineering
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
Total: 27 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Total: 96 hits
ReferencesLink to record
Permanent link

Direct link