Akshay Chandran , Fabien Evrard , Berend van Wachem
{"title":"A semi-analytical transient undisturbed velocity correction scheme for wall-bounded two-way coupled Euler-Lagrange simulations","authors":"Akshay Chandran , Fabien Evrard , Berend van Wachem","doi":"10.1016/j.jcp.2024.113496","DOIUrl":null,"url":null,"abstract":"<div><div>Force closure models employed in Euler-Lagrange (EL) point-particle simulations rely on the accurate estimation of the undisturbed fluid velocity at the particle center to evaluate the fluid forces on each particle. Due to the self-induced velocity disturbance of the particle in the fluid, two-way coupled EL simulations only have access to the disturbed velocity. The undisturbed velocity can be recovered if the particle-generated disturbance is estimated. In the present paper, we model the velocity disturbance generated by a regularized forcing near a planar wall, which, along with the temporal nature of the forcing, provides an estimate of the unsteady velocity disturbance of the particle near a planar wall. We use the analytical solution for a singular in-time transient Stokeslet near a planar wall (Felderhof <span><span>[1]</span></span>) and derive the corresponding time-persistent Stokeslets. The velocity disturbance due to a regularized forcing is then obtained numerically via a discrete convolution with the regularization kernel. The resulting Green's functions for parallel and perpendicular regularized forcing to the wall are stored as pre-computed temporal correction maps. By storing the time-dependent particle force on the fluid as fictitious particles, we estimate the unsteady velocity disturbance generated by the particle as a scalar product between the stored forces and the pre-computed Green's functions. Since the model depends on the analytical Green's function solution of the singular Stokeslet near a planar wall, the obtained velocity disturbance exactly satisfies the no-slip condition and does not require any fitted parameters to account for the rapid decay of the disturbance near the wall. The numerical evaluation of the convolution integral makes the present method suitable for arbitrary regularization kernels. Additionally, the generation of parallel and perpendicular correction maps enables to estimate the velocity disturbance due to particle motion in arbitrary directions relative to the local flow. The convergence of the method is studied on a fixed particle near a planar wall, and verification tests are performed in the Stokes regime on a settling particle parallel to a wall and a free-falling particle perpendicular to the wall.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113496"},"PeriodicalIF":3.8000,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999124007447","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Force closure models employed in Euler-Lagrange (EL) point-particle simulations rely on the accurate estimation of the undisturbed fluid velocity at the particle center to evaluate the fluid forces on each particle. Due to the self-induced velocity disturbance of the particle in the fluid, two-way coupled EL simulations only have access to the disturbed velocity. The undisturbed velocity can be recovered if the particle-generated disturbance is estimated. In the present paper, we model the velocity disturbance generated by a regularized forcing near a planar wall, which, along with the temporal nature of the forcing, provides an estimate of the unsteady velocity disturbance of the particle near a planar wall. We use the analytical solution for a singular in-time transient Stokeslet near a planar wall (Felderhof [1]) and derive the corresponding time-persistent Stokeslets. The velocity disturbance due to a regularized forcing is then obtained numerically via a discrete convolution with the regularization kernel. The resulting Green's functions for parallel and perpendicular regularized forcing to the wall are stored as pre-computed temporal correction maps. By storing the time-dependent particle force on the fluid as fictitious particles, we estimate the unsteady velocity disturbance generated by the particle as a scalar product between the stored forces and the pre-computed Green's functions. Since the model depends on the analytical Green's function solution of the singular Stokeslet near a planar wall, the obtained velocity disturbance exactly satisfies the no-slip condition and does not require any fitted parameters to account for the rapid decay of the disturbance near the wall. The numerical evaluation of the convolution integral makes the present method suitable for arbitrary regularization kernels. Additionally, the generation of parallel and perpendicular correction maps enables to estimate the velocity disturbance due to particle motion in arbitrary directions relative to the local flow. The convergence of the method is studied on a fixed particle near a planar wall, and verification tests are performed in the Stokes regime on a settling particle parallel to a wall and a free-falling particle perpendicular to the wall.
欧拉-拉格朗日(EL)点粒子模拟中采用的力闭合模型依赖于对粒子中心未受扰动的流体速度的精确估计,以评估每个粒子上的流体力。由于粒子在流体中的自致速度扰动,双向耦合 EL 模拟只能获得扰动速度。如果对粒子产生的扰动进行估计,则可以恢复未扰动的速度。在本文中,我们模拟了由平面壁附近的正则化作用力产生的速度扰动,它与作用力的时间性质一起,提供了对平面壁附近粒子的非稳态速度扰动的估计。我们使用平面壁附近的奇异时瞬态斯托克斯小波的解析解(Felderhof [1]),并推导出相应的时驻斯托克斯小波。然后,通过与正则化核的离散卷积,数值求得正则化强迫引起的速度扰动。由此得到的平行和垂直于壁面的正则化作用力的格林函数被存储为预先计算的时间校正图。通过将随时间变化的粒子对流体的作用力存储为虚构粒子,我们将粒子产生的非稳态速度扰动估算为存储的作用力与预先计算的格林函数之间的标量乘积。由于模型取决于平面壁附近奇异斯托克斯小波的格林函数解析解,因此得到的速度扰动完全满足无滑动条件,并且不需要任何拟合参数来解释壁附近扰动的快速衰减。卷积积分的数值评估使本方法适用于任意正则化核。此外,通过生成平行和垂直修正图,可以估算粒子在相对于局部流动的任意方向运动所产生的速度扰动。研究了该方法在平面壁附近的固定粒子上的收敛性,并在斯托克斯体系中对平行于壁的沉降粒子和垂直于壁的自由落体粒子进行了验证测试。
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.