Endre søk
RefereraExporteraLink to record
Permanent link

Direct link
Referera
Referensformat
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annet format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annet språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf
Corrected trapezoidal rules for a class of singular functions
KTH, Skolan för teknikvetenskap (SCI), Matematik (Inst.), Numerisk analys, NA.
KTH, Skolan för teknikvetenskap (SCI), Matematik (Inst.), Numerisk analys, NA. KTH, Centra, SeRC - Swedish e-Science Research Centre.ORCID-id: 0000-0002-6321-8619
KTH, Skolan för teknikvetenskap (SCI), Matematik (Inst.), Numerisk analys, NA. KTH, Centra, SeRC - Swedish e-Science Research Centre.
2014 (engelsk)Inngår i: IMA Journal of Numerical Analysis, ISSN 0272-4979, E-ISSN 1464-3642, Vol. 34, nr 4, s. 1509-1540Artikkel i tidsskrift (Fagfellevurdert) Published
Abstract [en]

A set of accurate quadrature rules applicable to a class of integrable functions with isolated singularities is designed and analysed theoretically in one and two dimensions. These quadrature rules are based on the trapezoidal rule with corrected quadrature weights for points in the vicinity of the singularity. To compute the correction weights, small-size ill-conditioned systems have to be solved. The convergence of the correction weights is accelerated by the use of compactly supported functions that annihilate boundary errors. Convergence proofs with error estimates for the resulting quadrature rules are given in both one and two dimensions. The tabulated weights are specific for the singularities under consideration, but the methodology extends to a large class of functions with integrable isolated singularities. Furthermore, in one dimension we have obtained a closed form expression based on which the modified weights can be computed directly.

sted, utgiver, år, opplag, sider
2014. Vol. 34, nr 4, s. 1509-1540
Emneord [en]
singular functions, quadrature methods, high order
HSV kategori
Identifikatorer
URN: urn:nbn:se:kth:diva-29760DOI: 10.1093/imanum/drt046ISI: 000343320900008Scopus ID: 2-s2.0-84904197227OAI: oai:DiVA.org:kth-29760DiVA, id: diva2:397500
Merknad

QC 20141121. Updated from submitted to published.

Tilgjengelig fra: 2011-02-15 Laget: 2011-02-15 Sist oppdatert: 2018-01-12bibliografisk kontrollert
Inngår i avhandling
1. Quadrature rules for boundary integral methods applied to Stokes flow
Åpne denne publikasjonen i ny fane eller vindu >>Quadrature rules for boundary integral methods applied to Stokes flow
2011 (engelsk)Licentiatavhandling, med artikler (Annet vitenskapelig)
Abstract [en]

 

Fluid phenomena dominated by viscous effects can, in many cases, be modeled by the Stokes equations. The boundary integral form of the Stokes equations reduces the number of degrees of freedom in a numerical discretization by reformulating the threedimensional problem to two-dimensional integral equations to be discretized over the boundaries of the domain. Hence for the study of objects immersed in a fluid, such as drops or elastic or solid particles, integral equations are to be discretized over the surfaces of these objects only. As outer boundaries or confinements are added these must also be included in the formulation. This work is focused on the development and validation of such a wall treatment. An inherent difficulty in the numerical treatment of boundary integrals for Stokes flow is the integration of the singular fundamental solution of the Stokes equations – the so called Stokeslet. To alleviate this problem we developed a set of high-order quadrature rules for the numerical integration of the Stokeslet over a flat surface. Such a quadrature rule was first designed for singularities of the type 1/|x|. To assess the convergence properties of this quadrature rule a theoretical analysis has been performed. The slightly more complicated singularity of the Stokeslet required certain modifications of the integration rule developed for 1/|x|. To validate the quadrature rule developed for the Stokeslet against a physical model we use it in a classical problem in fluid dynamics, the sedimentation of a sphere onto a flat plate. This involves a direct discretization of the plane wall and at the same time of the immersed sphere. Without any special treatment the algebraic system given by the discrete problem is quite memory consuming since matrix blocks are full. By exploring the structure of the block matrices that build up the system we have found that the wall discretization leads to a matrix which is generated by only three of its columns. This information together with certain preconditioning considerations allowed us to use the Schur complement method thus leading to a less memory expensive solution to the algebraic system. As a final step it is shown that the numerical simulations match the analytical solution, within the limitations of the model. This wall treatment can be easily extended to the problem of two parallel walls, and it is also shown that the simulation is in good agreement with some known results for the two parallel walls problem.

 

sted, utgiver, år, opplag, sider
Stockholm: KTH Royal Institute of Technology, 2011. s. vii, 30
Serie
Trita-CSC-A, ISSN 1653-5723 ; 2011:01
HSV kategori
Identifikatorer
urn:nbn:se:kth:diva-29763 (URN)978-91-7415-862-5 (ISBN)
Presentation
2011-02-15, Sal E32, KTH, Lindstedtsvägen 3, Stockholm, 13:00 (engelsk)
Opponent
Veileder
Forskningsfinansiär
Swedish e‐Science Research Center
Merknad
QC 20110215Tilgjengelig fra: 2011-02-15 Laget: 2011-02-15 Sist oppdatert: 2018-01-12bibliografisk kontrollert
2. Boundary integral methods for Stokes flow: Quadrature techniques and fast Ewald methods
Åpne denne publikasjonen i ny fane eller vindu >>Boundary integral methods for Stokes flow: Quadrature techniques and fast Ewald methods
2012 (engelsk)Doktoravhandling, med artikler (Annet vitenskapelig)
Abstract [en]

Fluid phenomena dominated by viscous effects can, in many cases, be modeled by the Stokes equations. The boundary integral form of the Stokes equations reduces the number of degrees of freedom in a numerical discretization by reformulating the three-dimensional problem to two-dimensional integral equations to be discretized over the boundaries of the domain.

Hence for the study of objects immersed in a fluid, such as drops or elastic/solid particles, integral equations are to be discretized over the surfaces of these objects only. As outer boundaries or confinements are added these must also be included in the formulation.

An inherent difficulty in the numerical treatment of boundary integrals for Stokes flow is the integration of the singular fundamental solution of the Stokes equations, e.g. the so called Stokeslet. To alleviate this problem we developed a set of high-order quadrature rules for the numerical integration of the Stokeslet over a flat surface. Such a quadrature rule was first designed for singularities of the type . To assess the convergence properties of this quadrature rule a theoretical analysis has been performed. The slightly more complicated singularity of the Stokeslet required certain modifications of the integration rule developed for . An extension of this type of quadrature rule to a cylindrical surface is also developed. These quadrature rules are tested also on physical problems that have an analytic solution in the literature.

Another difficulty associated with boundary integral problems is introduced by periodic boundary conditions. For a set of particles in a periodic domain periodicity is imposed by requiring that the motion of each particle has an added contribution from all periodic images of all particles all the way up to infinity. This leads to an infinite sum which is not absolutely convergent, and an additional physical constraint which removes the divergence needs to be imposed. The sum is decomposed into two fast converging sums, one that handles the short range interactions in real space and the other that sums up the long range interactions in Fourier space. Such decompositions are already available in the literature for kernels that are commonly used in boundary integral formulations. Here a decomposition in faster decaying sums than the ones present in the literature is derived for the periodic kernel of the stress tensor.

However the computational complexity of the sums, regardless of the decomposition they stem from, is . This complexity can be lowered using a fast summation method as we introduced here for simulating a sedimenting fiber suspension. The fast summation method was initially designed for point particles, which could be used for fibers discretized numerically almost without any changes. However, when two fibers are very close to each other, analytical integration is used to eliminate numerical inaccuracies due to the nearly singular behavior of the kernel and the real space part in the fast summation method was modified to allow for this analytical treatment. The method we have developed for sedimenting fiber suspensions allows for simulations in large periodic domains and we have performed a set of such simulations at a larger scale (larger domain/more fibers) than previously feasible.

sted, utgiver, år, opplag, sider
Stockholm: KTH Royal Institute of Technology, 2012. s. ix, 59
Serie
Trita-NA, ISSN 0348-2952 ; 2012:14
Emneord
boundary integral, Stokes flow, quadrature rule, Ewald decomposition
HSV kategori
Identifikatorer
urn:nbn:se:kth:diva-105540 (URN)978-91-7501-578-1 (ISBN)
Disputas
2012-12-14, D2, Lindstedtsvägen, 5, Stockholm, 10:15 (engelsk)
Opponent
Veileder
Forskningsfinansiär
Swedish e‐Science Research Center
Merknad

QC 20121122

Tilgjengelig fra: 2012-11-22 Laget: 2012-11-22 Sist oppdatert: 2013-04-09bibliografisk kontrollert

Open Access i DiVA

Fulltekst mangler i DiVA

Andre lenker

Forlagets fulltekstScopus

Personposter BETA

Runborg, Olof

Søk i DiVA

Av forfatter/redaktør
Marin, OanaRunborg, OlofTornberg, Anna-Karin
Av organisasjonen
I samme tidsskrift
IMA Journal of Numerical Analysis

Søk utenfor DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric

doi
urn-nbn
Totalt: 396 treff
RefereraExporteraLink to record
Permanent link

Direct link
Referera
Referensformat
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annet format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annet språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf