A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc

Background and Objective: The finite element method is widely used for studying the intervertebral disc at the organ level due to its ability to model complex geometries. An indispensable requirement for proper modelling of the intervertebral disc is a reliable porohyperelastic framework that captur...

ver descrição completa

Detalhes bibliográficos
Autores: Lialios, Dimitrios, Eguzkitza, Ane Betriz, Houzeaux, Guillaume, Casoni, Eva, Baumgartner, Laura, Noailly, Jérôme, Muñoz-Moya, Estefano, Gantenbein, Benjamin, Vázquez, Mariano
Formato: artículo
Estado:Versión publicada
Fecha de publicación:2025
País:España
Recursos:Universitat Pompeu Fabra
Repositorio:Repositorio Digital de la UPF
OAI Identifier:oai:repositori.upf.edu:10230/72094
Acesso em linha:http://hdl.handle.net/10230/72094
http://dx.doi.org/10.1016/j.cmpb.2024.108493
Access Level:acceso abierto
Palavra-chave:Intervertebral disc modelling
Porohyperelasticity
HPC systems
Finite element methods
Multiphysics coupling
id ES_1f844ebd884b37754a2fe342c166dbd5
oai_identifier_str oai:repositori.upf.edu:10230/72094
network_acronym_str ES
network_name_str España
repository_id_str
spelling A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral discLialios, DimitriosEguzkitza, Ane BetrizHouzeaux, GuillaumeCasoni, EvaBaumgartner, LauraNoailly, JérômeMuñoz-Moya, EstefanoGantenbein, BenjaminVázquez, MarianoIntervertebral disc modellingPorohyperelasticityHPC systemsFinite element methodsMultiphysics couplingBackground and Objective: The finite element method is widely used for studying the intervertebral disc at the organ level due to its ability to model complex geometries. An indispensable requirement for proper modelling of the intervertebral disc is a reliable porohyperelastic framework that captures the elaborate underlying mechanics. The increased complexity of such models requires significant computational power that is available within high-performance computing systems. The objective of this study is to present such a framework, validated both against literature and experiments, aiming to enable intervertebral disc research to benefit from state-of-the-art computational resources. Methods: In the context of this work, we implement a biphasic model that captures the mechanical response of the intricate, tissue-dependent models of the solid phase along with the hydrostatic pressure effects of the fluid phase. The tissue-dependent models involve the hyperelastic ground substance, fibrillar reinforcement, and osmotic swelling. The derived porohyperelastic, staggered scheme is implemented in Alya, a finite element code targeted at high-performance computing applications. The formulation is subsequently verified and validated by comparing the results of consolidation simulations with literature data for simulations and experiments using either generic or patient-specific geometries. Additionally, in-house experiments are replicated, evaluating the model’s ability to simulate alternating loading. Finally, the implementation’s circadian response is compared to previous implementation of similar material models in commercial software. Results: Results align well with experimental and literature findings in terms of disc height reduction (4% error), intradiscal pressure (14% error) and disc bulging. Validating the patient-specific geometry results in 4% and 7% deviation in measuring height loss. Simulations show excellent agreement with in-house experimental results, with less than 1% error regarding height reduction. Finally, the comparison to similar, published, earlier implementation in commercial software unveils excellent agreement of less than 1% error for the water content during circadian simulations. Simulation times are reported at 4 min per circadian cycle in the supercomputer Marenostrum V. Conclusions: This work presents a clear and validated formulation for simulating porohyperelastic materials based on assumptions that comply with the non-linear elasticity theory. The implementation in Alya enables intervertebral disc research to benefit from high-performance computing systems.Financial support was received from the Marie Skłodowska Curie International Training Network (ITN) “disc4all”, grant agreement #955735.Elsevier202520252025info:eu-repo/semantics/articleinfo:eu-repo/semantics/publishedVersionapplication/pdfapplication/pdfhttp://hdl.handle.net/10230/72094http://dx.doi.org/10.1016/j.cmpb.2024.108493reponame:Repositorio Digital de la UPFinstname:Universitat Pompeu FabraInglésComputer Methods and Programs in Biomedicine. 2025 Feb;259:108493info:eu-repo/grantAgreement/EC/FP7/955735© 2024 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).http://creativecommons.org/licenses/by/4.0/info:eu-repo/semantics/openAccessoai:repositori.upf.edu:10230/720942026-06-12T07:21:37Z
dc.title.none.fl_str_mv A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
title A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
spellingShingle A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
Lialios, Dimitrios
Intervertebral disc modelling
Porohyperelasticity
HPC systems
Finite element methods
Multiphysics coupling
title_short A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
title_full A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
title_fullStr A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
title_full_unstemmed A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
title_sort A porohyperelastic scheme targeted at high-performance computing frameworks for the simulation of the intervertebral disc
dc.creator.none.fl_str_mv Lialios, Dimitrios
Eguzkitza, Ane Betriz
Houzeaux, Guillaume
Casoni, Eva
Baumgartner, Laura
Noailly, Jérôme
Muñoz-Moya, Estefano
Gantenbein, Benjamin
Vázquez, Mariano
author Lialios, Dimitrios
author_facet Lialios, Dimitrios
Eguzkitza, Ane Betriz
Houzeaux, Guillaume
Casoni, Eva
Baumgartner, Laura
Noailly, Jérôme
Muñoz-Moya, Estefano
Gantenbein, Benjamin
Vázquez, Mariano
author_role author
author2 Eguzkitza, Ane Betriz
Houzeaux, Guillaume
Casoni, Eva
Baumgartner, Laura
Noailly, Jérôme
Muñoz-Moya, Estefano
Gantenbein, Benjamin
Vázquez, Mariano
author2_role author
author
author
author
author
author
author
author
dc.subject.none.fl_str_mv Intervertebral disc modelling
Porohyperelasticity
HPC systems
Finite element methods
Multiphysics coupling
topic Intervertebral disc modelling
Porohyperelasticity
HPC systems
Finite element methods
Multiphysics coupling
description Background and Objective: The finite element method is widely used for studying the intervertebral disc at the organ level due to its ability to model complex geometries. An indispensable requirement for proper modelling of the intervertebral disc is a reliable porohyperelastic framework that captures the elaborate underlying mechanics. The increased complexity of such models requires significant computational power that is available within high-performance computing systems. The objective of this study is to present such a framework, validated both against literature and experiments, aiming to enable intervertebral disc research to benefit from state-of-the-art computational resources. Methods: In the context of this work, we implement a biphasic model that captures the mechanical response of the intricate, tissue-dependent models of the solid phase along with the hydrostatic pressure effects of the fluid phase. The tissue-dependent models involve the hyperelastic ground substance, fibrillar reinforcement, and osmotic swelling. The derived porohyperelastic, staggered scheme is implemented in Alya, a finite element code targeted at high-performance computing applications. The formulation is subsequently verified and validated by comparing the results of consolidation simulations with literature data for simulations and experiments using either generic or patient-specific geometries. Additionally, in-house experiments are replicated, evaluating the model’s ability to simulate alternating loading. Finally, the implementation’s circadian response is compared to previous implementation of similar material models in commercial software. Results: Results align well with experimental and literature findings in terms of disc height reduction (4% error), intradiscal pressure (14% error) and disc bulging. Validating the patient-specific geometry results in 4% and 7% deviation in measuring height loss. Simulations show excellent agreement with in-house experimental results, with less than 1% error regarding height reduction. Finally, the comparison to similar, published, earlier implementation in commercial software unveils excellent agreement of less than 1% error for the water content during circadian simulations. Simulation times are reported at 4 min per circadian cycle in the supercomputer Marenostrum V. Conclusions: This work presents a clear and validated formulation for simulating porohyperelastic materials based on assumptions that comply with the non-linear elasticity theory. The implementation in Alya enables intervertebral disc research to benefit from high-performance computing systems.
publishDate 2025
dc.date.none.fl_str_mv 2025
2025
2025
dc.type.none.fl_str_mv info:eu-repo/semantics/article
info:eu-repo/semantics/publishedVersion
format article
status_str publishedVersion
dc.identifier.none.fl_str_mv http://hdl.handle.net/10230/72094
http://dx.doi.org/10.1016/j.cmpb.2024.108493
url http://hdl.handle.net/10230/72094
http://dx.doi.org/10.1016/j.cmpb.2024.108493
dc.language.none.fl_str_mv Inglés
language_invalid_str_mv Inglés
dc.relation.none.fl_str_mv Computer Methods and Programs in Biomedicine. 2025 Feb;259:108493
info:eu-repo/grantAgreement/EC/FP7/955735
dc.rights.none.fl_str_mv http://creativecommons.org/licenses/by/4.0/
info:eu-repo/semantics/openAccess
rights_invalid_str_mv http://creativecommons.org/licenses/by/4.0/
eu_rights_str_mv openAccess
dc.format.none.fl_str_mv application/pdf
application/pdf
dc.publisher.none.fl_str_mv Elsevier
publisher.none.fl_str_mv Elsevier
dc.source.none.fl_str_mv reponame:Repositorio Digital de la UPF
instname:Universitat Pompeu Fabra
instname_str Universitat Pompeu Fabra
reponame_str Repositorio Digital de la UPF
collection Repositorio Digital de la UPF
repository.name.fl_str_mv
repository.mail.fl_str_mv
_version_ 1869404404957839360
score 15,198674

Registos relacionados