1
|
Syed F, Khan S, Toma M. Modeling Dynamics of the Cardiovascular System Using Fluid-Structure Interaction Methods. BIOLOGY 2023; 12:1026. [PMID: 37508455 PMCID: PMC10376821 DOI: 10.3390/biology12071026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2023] [Revised: 07/07/2023] [Accepted: 07/18/2023] [Indexed: 07/30/2023]
Abstract
Using fluid-structure interaction algorithms to simulate the human circulatory system is an innovative approach that can provide valuable insights into cardiovascular dynamics. Fluid-structure interaction algorithms enable us to couple simulations of blood flow and mechanical responses of the blood vessels while taking into account interactions between fluid dynamics and structural behaviors of vessel walls, heart walls, or valves. In the context of the human circulatory system, these algorithms offer a more comprehensive representation by considering the complex interplay between blood flow and the elasticity of blood vessels. Algorithms that simulate fluid flow dynamics and the resulting forces exerted on vessel walls can capture phenomena such as wall deformation, arterial compliance, and the propagation of pressure waves throughout the cardiovascular system. These models enhance the understanding of vasculature properties in human anatomy. The utilization of fluid-structure interaction methods in combination with medical imaging can generate patient-specific models for individual patients to facilitate the process of devising treatment plans. This review evaluates current applications and implications of fluid-structure interaction algorithms with respect to the vasculature, while considering their potential role as a guidance tool for intervention procedures.
Collapse
Affiliation(s)
- Faiz Syed
- College of Osteopathic Medicine, New York Institute of Technology, Northern Boulevard, Old Westbury, NY 11568, USA
| | - Sahar Khan
- College of Osteopathic Medicine, New York Institute of Technology, Northern Boulevard, Old Westbury, NY 11568, USA
| | - Milan Toma
- College of Osteopathic Medicine, New York Institute of Technology, Northern Boulevard, Old Westbury, NY 11568, USA
| |
Collapse
|
2
|
Fluid-Structure Interaction Analyses of Biological Systems Using Smoothed-Particle Hydrodynamics. BIOLOGY 2021; 10:biology10030185. [PMID: 33801566 PMCID: PMC8001855 DOI: 10.3390/biology10030185] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Revised: 02/23/2021] [Accepted: 02/26/2021] [Indexed: 12/21/2022]
Abstract
Due to the inherent complexity of biological applications that more often than not include fluids and structures interacting together, the development of computational fluid-structure interaction models is necessary to achieve a quantitative understanding of their structure and function in both health and disease. The functions of biological structures usually include their interactions with the surrounding fluids. Hence, we contend that the use of fluid-structure interaction models in computational studies of biological systems is practical, if not necessary. The ultimate goal is to develop computational models to predict human biological processes. These models are meant to guide us through the multitude of possible diseases affecting our organs and lead to more effective methods for disease diagnosis, risk stratification, and therapy. This review paper summarizes computational models that use smoothed-particle hydrodynamics to simulate the fluid-structure interactions in complex biological systems.
Collapse
|
3
|
Wang X, Das R, Fernandez J. A surrogate mechanostatistical microstructural model to inform whole hip cortical bone remodelling. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3183. [PMID: 30701692 DOI: 10.1002/cnm.3183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2018] [Revised: 01/05/2019] [Accepted: 01/20/2019] [Indexed: 06/09/2023]
Abstract
Multiscale models of bone remodelling which incorporate clinically important features at a wide range of spatial scales are computationally intractable. In addition, there is a large gap in existing bone remodelling literature which deals with Haversian-level cortical bone remodelling behaviour and its effect on whole organ bone strength. In this study, we present an integration between continuum mechanics and surrogate modelling as an efficient circumvention to the problem of intractability, linked with an intricate 3D cortical bone remodelling algorithm which captures realistic microanatomical features at the Haversian level. The surrogate model predicted clear differentiation in long-term bone strength changes in clinically significant cases with extremely high computational efficiency.
Collapse
Affiliation(s)
- Xiaoming Wang
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Raj Das
- Department of Mechanical Engineering, The University of Auckland, Auckland, New Zealand
| | - Justin Fernandez
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
- Department of Engineering Science, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
4
|
Idkaidek A, Jasiuk I. Cortical bone fracture analysis using XFEM - case study. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2017; 33:e2809. [PMID: 27287280 DOI: 10.1002/cnm.2809] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Revised: 05/31/2016] [Accepted: 06/06/2016] [Indexed: 06/06/2023]
Abstract
We aim to achieve an accurate simulation of human cortical bone fracture using the extended finite element method within a commercial finite element software abaqus. A two-dimensional unit cell model of cortical bone is built based on a microscopy image of the mid-diaphysis of tibia of a 70-year-old human male donor. Each phase of this model, an interstitial bone, a cement line, and an osteon, are considered linear elastic and isotropic with material properties obtained by nanoindentation, taken from literature. The effect of using fracture analysis methods (cohesive segment approach versus linear elastic fracture mechanics approach), finite element type, and boundary conditions (traction, displacement, and mixed) on cortical bone crack initiation and propagation are studied. In this study cohesive segment damage evolution for a traction separation law based on energy and displacement is used. In addition, effects of the increment size and mesh density on analysis results are investigated. We find that both cohesive segment and linear elastic fracture mechanics approaches within the extended finite element method can effectively simulate cortical bone fracture. Mesh density and simulation increment size can influence analysis results when employing either approach, and using finer mesh and/or smaller increment size does not always provide more accurate results. Both approaches provide close but not identical results, and crack propagation speed is found to be slower when using the cohesive segment approach. Also, using reduced integration elements along with the cohesive segment approach decreases crack propagation speed compared with using full integration elements. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Ashraf Idkaidek
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, 1206 West Green Street, Urbana, IL, 61801, USA
| | - Iwona Jasiuk
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, 1206 West Green Street, Urbana, IL, 61801, USA
| |
Collapse
|
5
|
Role of cortical bone in hip fracture. BONEKEY REPORTS 2017; 6:867. [PMID: 28277562 DOI: 10.1038/bonekey.2016.82] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/11/2016] [Accepted: 09/03/2016] [Indexed: 12/23/2022]
Abstract
In this review, I consider the varied mechanisms in cortical bone that help preserve its integrity and how they deteriorate with aging. Aging affects cortical bone in two ways: extrinsically through its effects on the individual that modify its mechanical loading experience and 'milieu interieur'; and intrinsically through the prolonged cycle of remodelling and renewal extending to an estimated 20 years in the proximal femur. Healthy femoral cortex incorporates multiple mechanisms that help prevent fracture. These have been described at multiple length scales from the individual bone mineral crystal to the scale of the femur itself and appear to operate hierarchically. Each cortical bone fracture begins as a sub-microscopic crack that enlarges under mechanical load, for example, that imposed by a fall. In these conditions, a crack will enlarge explosively unless the cortical bone is intrinsically tough (the opposite of brittle). Toughness leads to microscopic crack deflection and bridging and may be increased by adequate regulation of both mineral crystal size and the heterogeneity of mineral and matrix phases. The role of osteocytes in optimising toughness is beginning to be worked out; but many osteocytes die in situ without triggering bone renewal over a 20-year cycle, with potential for increasing brittleness. Furthermore, the superolateral cortex of the proximal femur thins progressively during life, so increasing the risk of buckling during a fall. Besides preserving or increasing hip BMD, pharmaceutical treatments have class-specific effects on the toughness of cortical bone, although dietary and exercise-based interventions show early promise.
Collapse
|
6
|
Fernandez J, Zhang J, Heidlauf T, Sartori M, Besier T, Röhrle O, Lloyd D. Multiscale musculoskeletal modelling, data-model fusion and electromyography-informed modelling. Interface Focus 2016; 6:20150084. [PMID: 27051510 DOI: 10.1098/rsfs.2015.0084] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
This paper proposes methods and technologies that advance the state of the art for modelling the musculoskeletal system across the spatial and temporal scales; and storing these using efficient ontologies and tools. We present population-based modelling as an efficient method to rapidly generate individual morphology from only a few measurements and to learn from the ever-increasing supply of imaging data available. We present multiscale methods for continuum muscle and bone models; and efficient mechanostatistical methods, both continuum and particle-based, to bridge the scales. Finally, we examine both the importance that muscles play in bone remodelling stimuli and the latest muscle force prediction methods that use electromyography-assisted modelling techniques to compute musculoskeletal forces that best reflect the underlying neuromuscular activity. Our proposal is that, in order to have a clinically relevant virtual physiological human, (i) bone and muscle mechanics must be considered together; (ii) models should be trained on population data to permit rapid generation and use underlying principal modes that describe both muscle patterns and morphology; and (iii) these tools need to be available in an open-source repository so that the scientific community may use, personalize and contribute to the database of models.
Collapse
Affiliation(s)
- J Fernandez
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand; Department of Engineering Science, University of Auckland, Auckland, New Zealand
| | - J Zhang
- Auckland Bioengineering Institute , University of Auckland , Auckland , New Zealand
| | - T Heidlauf
- Institut für Mechanik (Bau) , University of Stuttgart , Stuttgart , Germany
| | - M Sartori
- Department of Neurorehabilitation Engineering , University Medical Center Göttingen , Göttingen , Germany
| | - T Besier
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand; Department of Engineering Science, University of Auckland, Auckland, New Zealand
| | - O Röhrle
- Institut für Mechanik (Bau) , University of Stuttgart , Stuttgart , Germany
| | - D Lloyd
- Centre for Musculoskeletal Research, Menzies Health Institute Queensland, Griffith University, Queensland, Australia; School of Rehabilitation Sciences, Griffith University, Queensland, Australia
| |
Collapse
|
7
|
Dalbeth N, Deacon M, Gamble GD, Mithraratne K, Fernandez J. Relationship between tissue stress during gait in healthy volunteers and patterns of urate deposition and bone erosion in gout: a biomechanical computational modelling study. RMD Open 2015; 1:e000101. [PMID: 26535140 PMCID: PMC4623368 DOI: 10.1136/rmdopen-2015-000101] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2015] [Revised: 06/29/2015] [Accepted: 07/02/2015] [Indexed: 01/19/2023] Open
Abstract
Objectives To determine whether patterns of high internal tissue stress during gait are associated with patterns of monosodium urate crystal deposition and bone erosion in gout. Methods We compared patterns of foot von Mises stress predicted computationally during gait in volunteers of normal and high body mass index (BMI) with patterns of urate deposition in gout and asymptomatic hyperuricaemia, and bone erosion in gout using dual-energy and conventional CT data. Results The highest average and peak von Mises stress during gait was observed at the third metatarsal (MT) head. Similar stress patterns were observed for high and low BMI groups. In contrast, for both urate deposition and bone erosion, the first MT head was most frequently affected, with very infrequent involvement of the third MT head. There was no clear relationship between average or peak von Mises stress patterns with patterns of urate deposition or bone erosion (−0.29>r<0.16). Addition of BMI into linear regression models did not alter the findings. Conclusions These data do not support the concept that elevated internal tissue stress during biomechanical loading plays an important role in patterns of monosodium urate crystal deposition or structural damage in gout.
Collapse
Affiliation(s)
- Nicola Dalbeth
- Department of Medicine , University of Auckland , Auckland , New Zealand
| | - Michelle Deacon
- Auckland Bioengineering Institute, University of Auckland , Auckland , New Zealand
| | - Gregory D Gamble
- Department of Medicine , University of Auckland , Auckland , New Zealand
| | - Kumar Mithraratne
- Auckland Bioengineering Institute, University of Auckland , Auckland , New Zealand
| | - Justin Fernandez
- Auckland Bioengineering Institute, University of Auckland , Auckland , New Zealand ; Department of Engineering Science , University of Auckland , Auckland , New Zealand
| |
Collapse
|
8
|
Wang X, Thomas CDL, Clement JG, Das R, Davies H, Fernandez JW. A mechanostatistical approach to cortical bone remodelling: an equine model. Biomech Model Mechanobiol 2015; 15:29-42. [PMID: 25862068 DOI: 10.1007/s10237-015-0669-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2014] [Accepted: 03/20/2015] [Indexed: 11/26/2022]
Abstract
In this study, the development of a mechanostatistical model of three-dimensional cortical bone remodelling informed with in vivo equine data is presented. The equine model was chosen as it is highly translational to the human condition due to similar Haversian systems, availability of in vivo bone strain and biomarker data, and furthermore, equine models are recommended by the US Federal Drugs Administration for comparative joint research. The model was derived from micro-computed tomography imaged specimens taken from the equine third metacarpal bone, and the Frost-based 'mechanostat' was informed from both in vivo strain gauges and biomarkers to estimate bone growth rates. The model also described the well-known 'cutting cone' phenomena where Haversian canals tunnel and replace bone. In order to make this model useful in practice, a partial least squares regression (PLSR) surrogate model was derived based on training data from finite element simulations with different loads. The PLSR model was able to predict microstructure and homogenised Young's modulus with errors less than 2.2% and 0.6%, respectively.
Collapse
Affiliation(s)
- X Wang
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.
| | - C D L Thomas
- Melbourne Dental School, University of Melbourne, Parkville, VIC, Australia
| | - J G Clement
- Melbourne Dental School, University of Melbourne, Parkville, VIC, Australia
| | - R Das
- Department of Mechanical Engineering, The University of Auckland, Auckland, New Zealand
| | - H Davies
- Faculty of Veterinary and Agricultural Science, University of Melbourne, Parkville, VIC, Australia
| | - J W Fernandez
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
- Department of Engineering Science, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
9
|
Das R, Collins A, Verma A, Fernandez J, Taylor M. Evaluating simulant materials for understanding cranial backspatter from a ballistic projectile. J Forensic Sci 2015; 60:627-37. [PMID: 25739515 DOI: 10.1111/1556-4029.12701] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2013] [Revised: 04/03/2014] [Accepted: 04/17/2014] [Indexed: 11/29/2022]
Abstract
In cranial wounds resulting from a gunshot, the study of backspatter patterns can provide information about the actual incidents by linking material to surrounding objects. This study investigates the physics of backspatter from a high-speed projectile impact and evaluates a range of simulant materials using impact tests. Next, we evaluate a mesh-free method called smoothed particle hydrodynamics (SPH) to model the splashing mechanism during backspatter. The study has shown that a projectile impact causes fragmentation at the impact site, while transferring momentum to fragmented particles. The particles travel along the path of least resistance, leading to partial material movement in the reverse direction of the projectile motion causing backspatter. Medium-density fiberboard is a better simulant for a human skull than polycarbonate, and lorica leather is a better simulant for a human skin than natural rubber. SPH is an effective numerical method for modeling the high-speed impact fracture and fragmentations.
Collapse
Affiliation(s)
- Raj Das
- Department of Mechanical Engineering, Centre for Advanced Composite Materials, University of Auckland, Auckland, 1010, New Zealand
| | | | | | | | | |
Collapse
|