A consistent sharp interface fictitious domain method for moving boundary problems with arbitrarily polyhedral mesh

A consistent, sharp interface fully Eulerian fictitious domain method is proposed in this article for moving boundary problems. In this method, a collocated finite volume method is used for the continuous phase; a geometry intersection method is employed for numerical integrals over the solid domain...

Full description

Saved in:
Bibliographic Details
Published in:International journal for numerical methods in fluids Vol. 93; no. 7; pp. 2065 - 2088
Main Authors: Chai, Guoliang, Wang, Le, Gu, Zhaolin, Yu, Chunlei, Zhang, Yigen, Shu, Qinglin, Su, Junwei
Format: Journal Article
Language:English
Published: Hoboken, USA John Wiley & Sons, Inc 01.07.2021
Wiley Subscription Services, Inc
Subjects:
Algorithms
Circular cylinders
Collocation methods
Continuity equation
Cylinders
Domains
Drag coefficient
Drag coefficients
fictitious domain method
Finite volume method
Fluctuations
Flux
Magnetism
Mathematical analysis
Momentum
Momentum equation
moving boundary problems
Numerical prediction
OpenFOAM
Oscillations
polyhedral mesh
sharp interface
spurious force oscillations
Tests
ISSN:0271-2091, 1097-0363
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:A consistent, sharp interface fully Eulerian fictitious domain method is proposed in this article for moving boundary problems. In this method, a collocated finite volume method is used for the continuous phase; a geometry intersection method is employed for numerical integrals over the solid domain and transport of the body force; the pseudo body force defined at “solid centers” ensures the algorithm consists of the body force between the continuous form and its discretization counterpart; an explicit flux correction on cell faces and resulting mass source is introduced into the continuity equation to lower noncontinuity errors in the velocity correction step. This method is valid for stationary and moving boundary problems with arbitrarily polyhedral mesh. Several numerical tests are carried out to validate the proposed method. A second‐order spatial accuracy is found in the flow around a cylinder case, and the spurious force oscillation is well suppressed for the in‐line oscillation of a circular cylinder case. The performances on different meshes are tested, and structured mesh yields the best result, polyhedral next, and tetrahedral worst. A serial of tests is further performed on structured mesh to verify the effect of three different features (i.e., storing the body force at the solid centers, flux correction, and whether including the body force in the momentum equation) on the numerical predictions. Numerical results show that, in the in‐line oscillation of a circular cylinder, “flux correction” can eliminate the large spikes in the drag coefficient, and “including the body force in the momentum equation” helps suppress the small oscillations. For other tests, “storing the body force at the solid centers” has enormous impacts on the final results of moving boundary problems, “flux correction” has little effects and the necessity of “including the body force in the momentum equation” is case dependent. Arbitrarily polyhedral mesh can be used for the proposed fictitious domain method. The sharpness of the solid–fluid boundary is acquired directly on the Eulerian background mesh without any additional mesh manipulation. Subgrid information is used for the discretization of the body force term to maintain the consistence, and for face flux correction to suppress the spurious force oscillations.
Bibliography:Funding information
National Major Science and Technology Projects of China, 2016ZX05011001; National Natural Science Foundation of China, 21306145
ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0271-2091
1097-0363
DOI:10.1002/fld.4965