skip to main content
research-article

A PML-based nonreflective boundary for free surface fluid animation

Published: 05 November 2010 Publication History

Abstract

This article presents a novel nonreflective boundary condition for the free surface incompressible Euler and Navier-Stokes equations. Boundaries of this type are very useful when, for example, simulating water flow around a ship moving over a wide ocean. Normally waves generated by the ship will reflect off of the boundaries of the simulation domain and as these reflected waves return towards the ship they will cause undesired interference patterns. By employing a Perfectly Matched Layer (PML) approach we have derived a boundary condition that absorbs incoming waves and thus efficiently prevents these undesired wave reflections. To solve the resulting boundary equations we present a fast and stable algorithm based on the stable fluids approach. Through numerical experiments we then show that our boundaries are significantly more effective than simpler reflection preventing techniques. We also provide a thorough analysis of the parameters involved in our boundary formulation and show how they effect wave absorption efficiency.

Supplementary Material

MP4 File (tp085_11.mp4)

References

[1]
Abarbanel, S. and Gottlieb, D. 1997. A mathematical analysis of the pml method. J. Comput. Phys. 134, 2, 357--363.
[2]
Adalsteinsson, D. and Sethian, J. A. 1995. A fast level set method for propagating interfaces. J. Comput. Phys. 118, 2, 269--277.
[3]
Baldauf, M. 2008. Stability analysis for linear discretisations of the advection equation with runge-kutta time integration. J. Comput. Phys. 227, 13, 6638--6659.
[4]
Bcache, E., Fauqueux, S., and Joly, P. 2003. Stability of perfectly matched layers, group velocities and anisotropic waves. J. Comput. Phys. 188, 2, 399--433.
[5]
Bcache, E., Dhia, A.-S. B.-B., and Legendre, G. 2004. Perfectly matched layers for the convected helmholtz equation. SIAM J. Numer. Anal. 42, 1, 409--433.
[6]
Berenger, J.-P. 1994. A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. 114, 2, 185--200.
[7]
Bcache, E., Fauqueux, S., and Joly, P. 2003. Stability of perfectly matched layers, group velocities and anisotropic waves. J. Comput. Phys. 188, 2, 399--433.
[8]
Chew, W. C. and Weedon, W. H. 1994. A 3-d perfectly matched medium from modified maxwell's equations with stretched coordinates. Microwave Opt. Tech. Lett. 7, 599--604.
[9]
Chorin, A. 1968. Numerical solution of navier-stokes equations. Math. Comp. 22, 745--762.
[10]
Enright, D., Fedkiw, R., Ferziger, J., and Mitchell, I. 2002. A hybrid particle level set method for improved interface capturing. J. Comput. Phys. 183, 1, 83--116.
[11]
Fedkiw, R., Stam, J., and Jensen, H. W. 2001. Visual simulation of smoke. In Proceedings of ACM SIGGRAPH Conference. 15--22.
[12]
Foster, N. and Fedkiw, R. 2001. Practical animation of liquids. In Proceedings of ACM SIGGRAPH Conference. ACM Press, 23--30.
[13]
Gottlieb, S. and Shu, C.-W. 1998. Total variation diminishing runge-kutta schemes. Math. Comput. 67, 221, 73--85.
[14]
Hagstrom, T., Goodrich, J., Nazarov, I., and Dodson, C. 2005. High-order methods and boundary conditions for simulating subsonic flows. In Proceedings of the 11th AIAA/CEAS Aeroacoustics Conference. AIAA--2005--2869.
[15]
Hagstrom, T. and Nazarov, I. 2002. Absorbing layers and radiation boundary conditions for jet flow simulations. In Proceedings of the 8th AIAA/CEAS Aeroacoustics Conference and Exhibit. 2002--2606.
[16]
Hagstrom, T. and Nazarov, I. 2003. Perfectly matched layers and radiation boundary conditions for shear flow calculations. In Proceedings of the 9th AIAA/CEAS Aeroacoustics Conference and Exhibit. 2003--3298.
[17]
He, X. and Luo, L.-S. 1997. Lattice boltzmann model for the incompressible navierstokes equation. J. Statis. Phys. 88, 3, 927--944.
[18]
Hein, S., Hohage, T., Koch, W., and Schberl, J. 2007. Acoustic resonances in a high-lift configuration. J. Fluid Mech. 582, -1, 179--202.
[19]
Hu, F. 2006. On the construction of pml absorbing boundary condition for the non-linear euler equations. In Proceedings of the 44th AIAA Aerospace Sciences Meeting and Exhibit. AIAA--2006--798.
[20]
Hu, F. Q. 1995. On absorbing boundary conditions for linearized euler equations by a perfectly matched layer. J. Comput. Phys. 129, 201--219.
[21]
Hu, F. Q. 1996. On absorbing boundary conditions for linearized euler equations by a perfectly matched layer. J. Comput. Phys. 129, 1, 201--219.
[22]
Hu, F. Q. 2001a. A stable, perfectly matched layer for linearized euler equations in unslit physical variables. J. Comput. Phys. 173, 2, 455--480.
[23]
Hu, F. Q. 2001b. A stable, perfectly matched layer for linearized euler equations in unsplit physical variables. J. Comput. Phys. 173, 2, 455--480.
[24]
Hu, F. Q. 2008. Development of pml absorbing boundary conditions for computational aeroacoustics: A progress review. Comput. Fluids 37, 4, 336--348.
[25]
Hu, F. Q., Hussaini, M. Y., and Manthey, J. L. 1996. Low-dissipation and low-dispersion runge-kutta schemes for computational acoustics. J. Comput. Phys. 124, 1, 177--191.
[26]
Hu, F. Q., Li, X. D., and Lin, D. K. 2008. Absorbing boundary conditions for nonlinear euler and navier-stokes equations based on the perfectly matched layer technique. J. Comput. Phys. 227, 9, 4398--4424.
[27]
Johnson, S. G. 2007. Notes on perfectly matched layers. http://www-math.mit.edu/~stevenj/18.369/pml.pdf
[28]
Li, X. D. and Gao, J. H. 2005. Numerical simulation of the generation mechanism of axisymmetric supersonic jet screech tones. Phys. Fluids 17, 8, 085105.
[29]
Li, X. D. and Gao, J. H. 2008. Numerical simulation of the three-dimensional screech phenomenon from a circular jet. Phys. Fluids 20, 3, 035101.
[30]
Liu, X., Osher, S., and Chan, T. 1994. Weighted essentially nonoscillatory schemes. J. Comput. Phys. 115, 200--212.
[31]
Monaghan, J. 1988. An introduction to sph. Comput. Phys. Comm. 48, 89--96.
[32]
Museth, K. and Leforestier, C. 1996. On the direct complex scaling of matrix elements expressed in a discrete variable representation: Application to molecular resonances. J. Chem. Phys. 104, 18, 7008--7014.
[33]
Neuhasuer, D. and Baer, M. 1989. The time?dependent schrždinger equation: Application of absorbing boundary conditions. J. Chem. Phys. 90, 8, 4351--4355.
[34]
Nielsen, M. B. and Museth, K. 2006. Dynamic Tubular Grid: An efficient data structure and algorithms for high resolution level sets. J. Sci. Comput. 26, 3, 261--299.
[35]
Nielsen, M. B., Nilsson, O., Söderström, A., and Museth, K. 2007. Out-of-core and compressed level set methods. ACM Trans. Graph. 26, 4, 16.
[36]
Osher, S. and Fedkiw, R. 2002. Level Set Methods and Dynamic Implicit Surfaces. Springer, Berlin.
[37]
Özyörük, Y. 2009. Numerical prediction of aft radiation of turbofan tones through exhaust jets. J. Sound Vibrat. 325, 1-2, 122--144.
[38]
Reinhardt, W. P. 1982. Complex coordinates in the theory of atomic and molecular structure and dynamics. Ann. Rev. Phys. Chem. 33, 223--255.
[39]
Richards, S. K., Zhang, X., Chen, X. X., and Nelson, P. A. 2004. The evaluation of non-reflecting boundary conditions for duct acoustic computation. J. Sound Vibrat. 270, 3, 539--557.
[40]
Shu, C. and Osher, S. 1988. Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys. 77, 439--471.
[41]
Söderström, A. and Museth, K. 2009. Nonreflective boundary conditions for incompressible free surface fluids. In Proceedings of the SIGGRAPH'09 Talks. ACM, New York 1--1.
[42]
Stam, J. 1999. Stable fluids. In Proceedings of SIGGRAPH Conference. 121--128.
[43]
Tam, C. K. and Webb, J. C. 1993. Dispersion-relation-preserving finite difference schemes for computational acoustics. J. Comput. Phys. 107, 2, 262--281.
[44]
Turkel, E. and Yefet, A. 1998. Absorbing pml boundary layers for wave-like equations. Appl. Numer. Math. 27, 4, 533--557.
[45]
Zhao, L. and Cangellaris, A. 1996. Gt-pml: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids. IEEE Trans. Microw. Theory Tech. 44, 12, 2555--2563.
[46]
Zhu, Y. and Bridson, R. 2005. Animating sand as a fluid. In Proceedings of the ACM SIGGRAPH Papers. ACM, New York. 965--972.

Cited By

View all

Recommendations

Comments

Information & Contributors

Information

Published In

cover image ACM Transactions on Graphics
ACM Transactions on Graphics  Volume 29, Issue 5
October 2010
58 pages
ISSN:0730-0301
EISSN:1557-7368
DOI:10.1145/1857907
Issue’s Table of Contents
Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected]

Publisher

Association for Computing Machinery

New York, NY, United States

Publication History

Published: 05 November 2010
Accepted: 01 August 2010
Revised: 01 April 2010
Received: 01 August 2009
Published in TOG Volume 29, Issue 5

Permissions

Request permissions for this article.

Check for updates

Author Tags

  1. Computational fluid dynamics
  2. Euler equations
  3. Navier-Stokes equations
  4. free surface
  5. nonreflecting boundary condition
  6. perfectly matched layer
  7. stable fluids

Qualifiers

  • Research-article
  • Research
  • Refereed

Contributors

Other Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

  • Downloads (Last 12 months)8
  • Downloads (Last 6 weeks)5
Reflects downloads up to 19 Feb 2025

Other Metrics

Citations

Cited By

View all
  • (2024)Non-Hermitian Absorbing Layers for Schrödinger's SmokeACM SIGGRAPH 2024 Posters10.1145/3641234.3671033(1-2)Online publication date: 25-Jul-2024
  • (2024)Trapped free surface waves for a Lamb–Oseen vortex flowJournal of Fluid Mechanics10.1017/jfm.2024.645997Online publication date: 17-Oct-2024
  • (2021)Ships, splashes, and waves on a vast oceanACM Transactions on Graphics10.1145/3478513.348049540:6(1-15)Online publication date: 10-Dec-2021
  • (2021)Kelvin transformations for simulations on infinite domainsACM Transactions on Graphics10.1145/3450626.345980940:4(1-15)Online publication date: 19-Jul-2021
  • (2021)A flux-interpolated advection scheme for fluid simulationThe Visual Computer10.1007/s00371-021-02187-2Online publication date: 10-Jun-2021
  • (2020)A practical octree liquid simulator with adaptive surface resolutionACM Transactions on Graphics10.1145/3386569.339246039:4(32:1-32:17)Online publication date: 12-Aug-2020
  • (2020)Asynchronous Eulerian Liquid SimulationComputer Graphics Forum10.1111/cgf.1390739:2(1-8)Online publication date: 13-Jul-2020
  • (2019)Robust high‐order semi‐implicit backward difference formulas for Navier‐Stokes equationsInternational Journal for Numerical Methods in Fluids10.1002/fld.472190:4(195-216)Online publication date: 6-Mar-2019
  • (2017)Fluxed animated boundary methodACM Transactions on Graphics10.1145/3072959.307359736:4(1-8)Online publication date: 20-Jul-2017
  • (2016)Generalized non-reflecting boundaries for fluid re-simulationACM Transactions on Graphics10.1145/2897824.292596335:4(1-7)Online publication date: 11-Jul-2016
  • Show More Cited By

View Options

Login options

Full Access

View options

PDF

View or Download as a PDF file.

PDF

eReader

View online with eReader.

eReader

Figures

Tables

Media

Share

Share

Share this Publication link

Share on social media