TY - JOUR
T1 - Stencil Penalty approach based constraint immersed boundary method
AU - Bale, Rahul
AU - Patankar, Neelesh A.
AU - Jansson, Niclas
AU - Onishi, Keiji
AU - Tsubokura, Makoto
N1 - Funding Information:
R.B, K.O, N.J, and M.T acknowledge support of this research by MEXT as “Priority Issue on Post-K computer” (Development of innovative design and production processes) and the computational resources of the K computer provided by the RIKEN Center for Computational Science through the HPCI System Research project (Project ID: hp160033, hp160032, hp170276, hp190182 and hp180192). N.A.P acknowledges support from the National Science Foundation (NSF award SI2- SSI-1450374 ).
Publisher Copyright:
© 2020 Elsevier Ltd
PY - 2020/3/30
Y1 - 2020/3/30
N2 - The constraint-based immersed boundary (cIB) method has been shown to be accurate between low and moderate Reynolds number (Re) flows when the immersed body constraint is imposed as a volumetric constraint force. When the IB is modelled as a zero-thickness interface, where it is no longer possible to model a volumetric constraint force, we found that cIB is not able to produce accurate results. The main source of inaccuracies in the cIB method is the distribution of the pressure field around the IB surface. An IB surface results in a jump in the pressure field across the IB. Evaluation of the discrete gradient of pressure close to the IB leads to a pressure gradient that does not satisfy the Neumann boundary condition for pressure at the IB. Furthermore, a non-zero discrete pressure gradient on the IB results in spurious flow at grid points close to the IB. We present a novel numerical formulation which adapts the cIB formulation for ‘zero-thickness’ immersed bodies. To impose the Neumann boundary condition on pressure on the IB more accurately, we introduce an additional body force to the momentum equation. A WENO based stencil penalization technique is used to define the new force term. Due to the more accurate imposition on the Neumann pressure boundary condition on the IB, spurious flow is reduced and the accuracy of no penetration velocity boundary condition on the IB is improved.
AB - The constraint-based immersed boundary (cIB) method has been shown to be accurate between low and moderate Reynolds number (Re) flows when the immersed body constraint is imposed as a volumetric constraint force. When the IB is modelled as a zero-thickness interface, where it is no longer possible to model a volumetric constraint force, we found that cIB is not able to produce accurate results. The main source of inaccuracies in the cIB method is the distribution of the pressure field around the IB surface. An IB surface results in a jump in the pressure field across the IB. Evaluation of the discrete gradient of pressure close to the IB leads to a pressure gradient that does not satisfy the Neumann boundary condition for pressure at the IB. Furthermore, a non-zero discrete pressure gradient on the IB results in spurious flow at grid points close to the IB. We present a novel numerical formulation which adapts the cIB formulation for ‘zero-thickness’ immersed bodies. To impose the Neumann boundary condition on pressure on the IB more accurately, we introduce an additional body force to the momentum equation. A WENO based stencil penalization technique is used to define the new force term. Due to the more accurate imposition on the Neumann pressure boundary condition on the IB, spurious flow is reduced and the accuracy of no penetration velocity boundary condition on the IB is improved.
UR - http://www.scopus.com/inward/record.url?scp=85078984288&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85078984288&partnerID=8YFLogxK
U2 - 10.1016/j.compfluid.2020.104457
DO - 10.1016/j.compfluid.2020.104457
M3 - Article
AN - SCOPUS:85078984288
VL - 200
JO - Computers and Fluids
JF - Computers and Fluids
SN - 0045-7930
M1 - 104457
ER -