Volume 15, Issue 1
Reduction of the Regularization Error of the Method of Regularized Stokeslets for a Rigid Object Immersed in a Three-Dimensional Stokes Flow

Commun. Comput. Phys., 15 (2014), pp. 126-152.

Published online: 2014-01

Cited by

Export citation
• Abstract

We focus on the problem of evaluating the velocity field outside a solid object moving in an incompressible Stokes flow using the boundary integral formulation. For points near the boundary, the integral is nearly singular, and accurate computation of the velocity is not routine. One way to overcome this problem is to regularize the integral kernel. The method of regularized Stokeslet (MRS) is a systematic way to regularize the kernel in this situation. For a specific blob function which is widely used, the error of the MRS is only of first order with respect to the blob parameter. We prove that this is the case for radial blob functions with decay property $ϕ(r)$=$\mathcal{O}$($r^{−3−α}$) when $r→∞$ for some constant $α>1$. We then find a class of blob functions for which the leading local error term can be removed to get second and third order errors with respect to blob parameter. Since the addition of these terms might give a flow field that is not divergence free, we introduce a modification of these terms to make the divergence of the corrected flow field close to zero while keeping the desired accuracy. Furthermore, these dominant terms are explicitly expressed in terms of blob function and so the computation time is negligible.

• Keywords

• AMS Subject Headings

• BibTex
• RIS
• TXT
@Article{CiCP-15-126, author = {}, title = {Reduction of the Regularization Error of the Method of Regularized Stokeslets for a Rigid Object Immersed in a Three-Dimensional Stokes Flow}, journal = {Communications in Computational Physics}, year = {2014}, volume = {15}, number = {1}, pages = {126--152}, abstract = {

We focus on the problem of evaluating the velocity field outside a solid object moving in an incompressible Stokes flow using the boundary integral formulation. For points near the boundary, the integral is nearly singular, and accurate computation of the velocity is not routine. One way to overcome this problem is to regularize the integral kernel. The method of regularized Stokeslet (MRS) is a systematic way to regularize the kernel in this situation. For a specific blob function which is widely used, the error of the MRS is only of first order with respect to the blob parameter. We prove that this is the case for radial blob functions with decay property $ϕ(r)$=$\mathcal{O}$($r^{−3−α}$) when $r→∞$ for some constant $α>1$. We then find a class of blob functions for which the leading local error term can be removed to get second and third order errors with respect to blob parameter. Since the addition of these terms might give a flow field that is not divergence free, we introduce a modification of these terms to make the divergence of the corrected flow field close to zero while keeping the desired accuracy. Furthermore, these dominant terms are explicitly expressed in terms of blob function and so the computation time is negligible.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.021112.290413a}, url = {http://global-sci.org/intro/article_detail/cicp/7090.html} }
TY - JOUR T1 - Reduction of the Regularization Error of the Method of Regularized Stokeslets for a Rigid Object Immersed in a Three-Dimensional Stokes Flow JO - Communications in Computational Physics VL - 1 SP - 126 EP - 152 PY - 2014 DA - 2014/01 SN - 15 DO - http://doi.org/10.4208/cicp.021112.290413a UR - https://global-sci.org/intro/article_detail/cicp/7090.html KW - AB -

We focus on the problem of evaluating the velocity field outside a solid object moving in an incompressible Stokes flow using the boundary integral formulation. For points near the boundary, the integral is nearly singular, and accurate computation of the velocity is not routine. One way to overcome this problem is to regularize the integral kernel. The method of regularized Stokeslet (MRS) is a systematic way to regularize the kernel in this situation. For a specific blob function which is widely used, the error of the MRS is only of first order with respect to the blob parameter. We prove that this is the case for radial blob functions with decay property $ϕ(r)$=$\mathcal{O}$($r^{−3−α}$) when $r→∞$ for some constant $α>1$. We then find a class of blob functions for which the leading local error term can be removed to get second and third order errors with respect to blob parameter. Since the addition of these terms might give a flow field that is not divergence free, we introduce a modification of these terms to make the divergence of the corrected flow field close to zero while keeping the desired accuracy. Furthermore, these dominant terms are explicitly expressed in terms of blob function and so the computation time is negligible.

Hoang-Ngan Nguyen & Ricardo Cortez. (2020). Reduction of the Regularization Error of the Method of Regularized Stokeslets for a Rigid Object Immersed in a Three-Dimensional Stokes Flow. Communications in Computational Physics. 15 (1). 126-152. doi:10.4208/cicp.021112.290413a
Copy to clipboard
The citation has been copied to your clipboard