Display options
Share it on

J Biomech Eng. 2018 Feb 01;140(2). doi: 10.1115/1.4038716.

Finite Element Framework for Computational Fluid Dynamics in FEBio.

Journal of biomechanical engineering

Gerard A Ateshian, Jay J Shim, Steve A Maas, Jeffrey A Weiss

Affiliations

  1. Department of Mechanical Engineering, Columbia University, New York, NY 10027.
  2. Department of Bioengineering, University of Utah, Salt Lake City, UT 84112.

PMID: 29238817 PMCID: PMC5816258 DOI: 10.1115/1.4038716

Abstract

The mechanics of biological fluids is an important topic in biomechanics, often requiring the use of computational tools to analyze problems with realistic geometries and material properties. This study describes the formulation and implementation of a finite element framework for computational fluid dynamics (CFD) in FEBio, a free software designed to meet the computational needs of the biomechanics and biophysics communities. This formulation models nearly incompressible flow with a compressible isothermal formulation that uses a physically realistic value for the fluid bulk modulus. It employs fluid velocity and dilatation as essential variables: The virtual work integral enforces the balance of linear momentum and the kinematic constraint between fluid velocity and dilatation, while fluid density varies with dilatation as prescribed by the axiom of mass balance. Using this approach, equal-order interpolations may be used for both essential variables over each element, contrary to traditional mixed formulations that must explicitly satisfy the inf-sup condition. The formulation accommodates Newtonian and non-Newtonian viscous responses as well as inviscid fluids. The efficiency of numerical solutions is enhanced using Broyden's quasi-Newton method. The results of finite element simulations were verified using well-documented benchmark problems as well as comparisons with other free and commercial codes. These analyses demonstrated that the novel formulation introduced in FEBio could successfully reproduce the results of other codes. The analogy between this CFD formulation and standard finite element formulations for solid mechanics makes it suitable for future extension to fluid-structure interactions (FSIs).

References

  1. J Biomech Eng. 1993 Feb;115(1):1-12 - PubMed
  2. Comput Struct. 2007 Jun;85(11-14):1114-1120 - PubMed
  3. J Biomech. 2012 Apr 5;45(6):1023-7 - PubMed
  4. Biorheology. 1991;28(3-4):241-62 - PubMed
  5. J Comput Phys. 2008 Nov 20;227(22):9303-9332 - PubMed
  6. Comput Biol Med. 2011 Feb;41(2):67-75 - PubMed
  7. J Biomech Eng. 2012 Jan;134(1):011005 - PubMed
  8. J Biomech Eng. 2004 Dec;126(6):760-9 - PubMed
  9. J Biomech. 2015 Apr 13;48(6):941-7 - PubMed
  10. J Biomech Eng. 2011 Jul;133(7):071003 - PubMed
  11. ASAIO J. 2002 May-Jun;48(3):260-7 - PubMed
  12. Cell Mol Bioeng. 2010 Sep 1;3(3):296-306 - PubMed
  13. J Biomech Eng. 2011 Feb;133(2):021001 - PubMed
  14. Annu Rev Fluid Mech. 2009 Jan 1;41:347-374 - PubMed
  15. Interface Focus. 2016 Feb 6;6(1):20150063 - PubMed
  16. Ann Biomed Eng. 2017 Mar;45(3):525-541 - PubMed
  17. J Theor Biol. 2001 May 21;210(2):221-36 - PubMed
  18. Med Eng Phys. 2005 Dec;27(10):871-83 - PubMed
  19. Ann Biomed Eng. 2011 Jan;39(1):484-96 - PubMed
  20. J Biomech Eng. 2011 Aug;133(8):081005 - PubMed
  21. J Biomech Eng. 2013 Nov;135(11):111001 - PubMed
  22. J Biomech Eng. 2010 Jun;132(6):061006 - PubMed
  23. J Biomech. 1976;9(7):439-48 - PubMed
  24. Biomech Model Mechanobiol. 2016 Feb;15(1):9-28 - PubMed
  25. Ann Biomed Eng. 2010 Apr;38(4):1426-41 - PubMed
  26. Proc Inst Mech Eng H. 2014 Oct;228(10):985-95 - PubMed
  27. Biomech Model Mechanobiol. 2014 Oct;13(5):1105-20 - PubMed
  28. Basic Clin Neurosci. 2013 Winter;4(1):64-75 - PubMed

MeSH terms

Publication Types

Grant support