...
首页> 外文期刊>Journal of Computational Physics >A fast numerical method for computing doubly-periodic regularized Stokes flow in 3D
【24h】

A fast numerical method for computing doubly-periodic regularized Stokes flow in 3D

机译:计算3D双周期正则Stokes流的快速数值方法

获取原文
获取原文并翻译 | 示例
           

摘要

A new numerical method for computing three-dimensional Stokes flow driven by a doubly-periodic array of regularized forces is presented. The method is based on deriving an analytical representation of a regularized Green's function in Fourier space. Then only an inverse fast Fourier transform (inverse FFT) has to be computed to determine the fluid velocity on a grid in the physical domain. The velocity at other points can be interpolated from this grid. Accuracy is verified by comparing numerical results to a solution that is independent of the method. Although the regularized forces lead to a smooth velocity field, the Green's function may contain rapid transitions that are not captured properly on a coarse grid. In that case, an Ewald splitting technique is used to compute the grid-resolved part of the flow using an inverse FFT and a sum in physical space for the localized part of the velocity. The splitting parameter can be chosen as small as a few grid cells, which makes the sum in physical space converge extremely fast. We present numerical examples that demonstrate that fact. In some cases, when the grid size is sufficiently small compared to the regularization parameter, the Ewald splitting is not needed.
机译:提出了一种新的数值方法,用于计算由正周期力的双周期阵列驱动的三维斯托克斯流。该方法基于推导傅立叶空间中正则化格林函数的解析表示。然后,仅需计算快速傅立叶逆变换(FFT)即可确定物理域中网格上的流体速度。可以从该网格内插其他点的速度。通过将数值结果与独立于该方法的解决方案进行比较,可以验证准确性。尽管规则力导致平滑的速度场,但格林函数可能包含快速跳变,而在粗糙的网格上无法正确捕获。在那种情况下,使用Ewald分裂技术使用逆FFT和速度局部化的物理空间总和来计算流的网格解析部分。分裂参数可以选择为只有几个网格单元,从而使物理空间中的总和收敛极快。我们提供了数值示例来证明这一事实。在某些情况下,当栅格大小与正则化参数相比足够小时,不需要进行Ewald分裂。

著录项

相似文献

  • 外文文献
  • 中文文献
  • 专利
获取原文

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号