Multiphysics simulation of the effect of leaflet thickness inhomogeneity and material anisotropy on the stress-strain distribution on the aortic valve

The impact of thickness heterogeneity on soft tissue biomechanics: a novel measurement technique and a demonstration on heart valve tissue

The mechanical properties of soft tissues are driven by their complex, heterogeneous composition and structure. Interestingly, studies of soft tissue biomechanics often ignore spatial heterogeneity. In our work, we are therefore interested in exploring the impact of tissue heterogeneity on the mechanical properties of soft tissues. Therein, we specifically focus on soft tissue heterogeneity arising from spatially varying thickness. To this end, our first goal is to develop a non-destructive measurement technique that has a high spatial resolution, provides continuous thickness maps, and is fast. Our secondary goal is to demonstrate that including spatial variation in thickness is important to the accuracy of biomechanical analyses. To this end, we use mitral valve leaflet tissue as our model system. To attain our first goal, we identify a soft tissue-specific contrast protocol that enables thickness measurements using a Keyence profilometer. We also show that this protocol does not affect our tissues' mechanical properties. To attain our second goal, we conduct virtual biaxial, bending, and buckling tests on our model tissue both ignoring and considering spatial variation in thickness. Thereby, we show that the assumption of average, homogeneous thickness distributions significantly alters the results of biomechanical analyses when compared to including true, spatially varying thickness distributions. In conclusion, our work provides a novel measurement technique that can capture continuous thickness maps non-invasively, at high resolution, and in a short time. Our work also demonstrates the importance of including heterogeneous thickness in biomechanical analyses of soft tissues.

Strategies for Development of Synthetic Heart Valve Tissue Engineering Scaffolds

The current clinical solutions, including mechanical and bioprosthetic valves for valvular heart diseases, are plagued by coagulation, calcification, nondurability, and the inability to grow with patients. The tissue engineering approach attempts to resolve these shortcomings by producing heart valve scaffolds that may deliver patients a life-long solution. Heart valve scaffolds serve as a three-dimensional support structure made of biocompatible materials that provide adequate porosity for cell infiltration, and nutrient and waste transport, sponsor cell adhesion, proliferation, and differentiation, and allow for extracellular matrix production that together contributes to the generation of functional neotissue. The foundation of successful heart valve tissue engineering is replicating native heart valve architecture, mechanics, and cellular attributes through appropriate biomaterials and scaffold designs. This article reviews biomaterials, the fabrication of heart valve scaffolds, and their in-vitro and in-vivo evaluations applied for heart valve tissue engineering.

A Doppler-exclusive non-invasive computational diagnostic framework for personalized transcatheter aortic valve replacement

Given the associated risks with transcatheter aortic valve replacement (TAVR), it is crucial to determine how the implant will affect the valve dynamics and cardiac function, and if TAVR will improve or worsen the outcome of the patient. Effective treatment strategies, indeed, rely heavily on the complete understanding of the valve dynamics. We developed an innovative Doppler-exclusive non-invasive computational framework that can function as a diagnostic tool to assess valve dynamics in patients with aortic stenosis in both pre- and post-TAVR status. Clinical Doppler pressure was reduced by TAVR (52.2 ± 20.4 vs. 17.3 ± 13.8 [mmHg], p < 0.001), but it was not always accompanied by improvements in valve dynamics and left ventricle (LV) hemodynamics metrics. TAVR had no effect on LV workload in 4 patients, and LV workload post-TAVR significantly rose in 4 other patients. Despite the group level improvements in maximum LV pressure (166.4 ± 32.2 vs 131.4 ± 16.9 [mmHg], p < 0.05), only 5 of the 12 patients (41%) had a decrease in LV pressure. Moreover, TAVR did not always improve valve dynamics. TAVR did not necessarily result in a decrease (in 9 out of 12 patients investigated in this study) in major principal stress on the aortic valve leaflets which is one of the main contributors in valve degeneration and, consequently, failure of heart valves. Diastolic stresses increased significantly post-TAVR (34%, 109% and 81%, p < 0.001) for each left, right and non-coronary leaflets respectively. Moreover, we quantified the stiffness and material properties of aortic valve leaflets which correspond with the reduced calcified region average stiffness among leaflets (66%, 74% and 62%; p < 0.001; N = 12). Valve dynamics post-intervention should be quantified and monitored to ensure the improvement of patient conditions and prevent any further complications. Improper evaluation of biomechanical valve features pre-intervention as well as post-intervention may result in harmful effects post-TAVR in patients including paravalvular leaks, valve degeneration, failure of TAVR and heart failure.

An ultrasound-exclusive non-invasive computational diagnostic framework for personalized cardiology of aortic valve stenosis

Aortic stenosis (AS) is an acute and chronic cardiovascular disease and If left untreated, 50% of these patients will die within two years of developing symptoms. AS is characterized as the stiffening of the aortic valve leaflets which restricts their motion and prevents the proper opening under transvalvular pressure. Assessments of the valve dynamics, if available, would provide valuable information about the patient's state of cardiac deterioration as well as heart recovery and can have incredible impacts on patient care, planning interventions and making critical clinical decisions with life-threatening risks. Despite remarkable advancements in medical imaging, there are no clinical tools available to quantify valve dynamics invasively or noninvasively. In this study, we developed a highly innovative ultrasound-based non-invasive computational framework that can function as a diagnostic tool to assess valve dynamics (e.g. transient 3-D distribution of stress and displacement, 3-D deformed shape of leaflets, geometric orifice area and angular positions of leaflets) for patients with AS at no risk to the patients. Such a diagnostic tool considers the local valve dynamics and the global circulatory system to provide a platform for testing the intervention scenarios and evaluating their effects. We used clinical data of 12 patients with AS not only to validate the proposed framework but also to demonstrate its diagnostic abilities by providing novel analyses and interpretations of clinical data in both pre and post intervention states. We used transthoracic echocardiogram (TTE) data for the developments and transesophageal echocardiography (TEE) data for validation.

Design of an aortic polymeric valve with asymmetric leaflets and evaluation of its performance by finite element method

BACKGROUND

The geometry of leaflets plays a significant role in prosthetic valves' (PVs) performance. Typically, natural aortic valves have three unequal leaflets, which differ in size. The present study aims to design an asymmetric tri-leaflet polymeric valve with one large and two small leaflets based on commissure lengths and leaflet eccentricities.

METHODS

Eccentricity was related to commissure lengths based on the deformation of the free margins for the fully-opened state of leaflets. The polystyrene-block-polyethylene-polypropylene-block-polystyrene polymer characterized the material properties of the leaflets. The Finite Element Method (FEM) was used to evaluate performance parameters, including maximum geometric orifice area (GOA), average GOA, maximum von Mises stress, and leaflet's coaptation surface area (CSA).

RESULTS

Asymmetric valves with no eccentricity provided a low level of GOA because the asymmetric form of small leaflets caused them to close faster than the large leaflet, leading to a sudden drop in the GOA during systole. As the radial curve tends towards a straight line, an undesirable coaptation occurs, and peak stress increases despite higher GOAs. A new radial curve consisting of two straight lines connected by an arc that provided 25.64 mm² coaptation surface area (CAS) and 117.54 mm² average GOA, was proposed to improve coaptation and GOA.

CONCLUSION

The radial curve of leaflets affects the valve's performance more than other geometric parameters. The combination of straight lines and arcs for radial curves was selected as the reference model for asymmetric valves with one large and two small leaflets.

An study on the influence of collagen fiber directions in TAVs performance using FEM

Transcatheter Aortic Valve Implantation (TAVI) or Replacement (TAVR) is a promising treatment for aortic valve stenosis, consisting of a procedure to replace a damaged native aortic valve by a bioprosthetic one. This replacement valve control the flow of blood using leaflets that are similar to the ones of a native aortic valve. Commonly manufactured using bovine or porcine pericardium, it is a tissue histologically composed of collagen fibers embedded into a nearly-isotropic matrix, where their distribution makes the pericardium behave as an anisotropic hyperelastic material. Because of such complicated behavior, bioprosthetic pericardium valves are, as expected, sensitive to the distribution and orientation of these fibers in such device. Therefore, the objective of this work is a thorough systematic study on the influence of these fibers' distribution. First, a Finite Element model of a bioprosthetic valve is generated; then, a material routine to accurately describe the behavior of pericardium is implemented in a commercial software package; in addition, a dedicated algorithm to specify the direction of fibers is developed. Finally, a systematic study on the influence that fiber orientations have on the overall behavior of the TAV is performed. As a result of this study, two extreme behaviors are highlighted depending on the preferential orientation of collagen fibers; namely, one with fibers in circumferential direction and the opposite with fibers in an axial orientation. Then, it is concluded that the behavior of fibers in circumferential direction is very sensitive to small variations of the orientation angle, whereas such orientation is not as determining when the aim is to achieve a behavior near to the one corresponding with axial orientation.

Local and global growth and remodeling in calcific aortic valve disease and aging

Aging and calcific aortic valve disease (CAVD) are the main factors leading to aortic stenosis. Both processes are accompanied by growth and remodeling pathways that play a crucial role in aortic valve pathophysiology. Herein, a computational growth and remodeling (G&R) framework was developed to investigate the effects of aging and calcification on aortic valve dynamics. Particularly, an algorithm was developed to couple the global growth and stiffening of the aortic valve due to aging and the local growth and stiffening due to calcification with the aortic valve transient dynamics. The aortic valve dynamics during baseline were validated with available data in the literature. Subsequently, the changes in aortic valve dynamic patterns during aging and CAVD progression were studied. The results revealed the patterns in geometric orifice area reduction and an increase in the valve stress during local and global growth and remodeling of the aortic valve. The proposed algorithm provides a framework to couple mechanobiology models of disease growth with tissue-scale transient structural mechanics models to study the biomechanical changes during cardiovascular disease growth and aging.

Porcine and bovine aortic valve comparison for surgical optimization: A fluid-structure interaction modeling study

BACKGROUND

Porcine aortic valve (PAV) and bovine aortic valve (BAV) are commonly used in aortic valve replacement (AVR) surgeries. A detailed comparison for their hemodynamic and structural stress/strain performances would help to better understand valve cardiac function and select valve type and size for AVR outcome optimizations.

METHODS

Eight fluid-structure interaction models were constructed to compare hemodynamic and stress/strain behaviors of PAV and BAV with 4 sizes (19, 21, 23, and 25 mm). Blood flow velocity, systolic cross-valve pressure gradient (SCVPG), geometric orifice area (GOA), flow shear stresses (FSS), and stress/strain were obtained for comparison.

RESULTS

Compared with PAV, BAV has better hemodynamic performance, with lower maximum flow velocity (7.17%) and pressure (9.82%), smaller pressure gradient (mean and peak SCVPG: 8.92% and 9.28%), larger GOA (9.56%) and lower FSS (6.61%). The averages of the mean and peak net pressure gradient values from 4 BAV models were 8.10% and 8.35% lower than that from PAV models. Larger valve sizes for both PAV and BAV had improved hemodynamic performance. Maximum flow velocity, pressure, mean SCVPG and maximum FSS from 25 mm BAV were 36.80%, 15.81%, 39.05% and 38.83% lower than those from 19 mm BAV. The GOA of PAV and BAV 25 mm Valve were 43.75% and 33.07% larger than 19 mm valves, respectively. BAV has lower stress on the leaflets than PAV.

CONCLUSIONS

BAV had better hemodynamic performance and lower leaflets stress than PAV. More patient studies are needed to validate our findings.

Subject-specific multiscale modeling of aortic valve biomechanics

A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.36-0.87$$\end{document}0.36-0.87) with a strong correlation with the organ length-scale radial strain (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^{2}= 0.95$$\end{document}R2=0.95); conversely, circumferential cell strains spanned a very narrow range (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.75-0.88$$\end{document}0.75-0.88) showing no correlation with the circumferential strain at the organ level (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^{2}= 0.02$$\end{document}R2=0.02). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.

Hemodynamic effects of support modes of LVADs on the aortic valve

As the alternative treatment for heart failure, left ventricular assist devices (LVADs) have been widely applied to clinical practice. However, the effects of the support modes of LVADs on the biomechanical states of the aortic valve are still poorly understood. Hence, the present study investigates such effects and proposes a novel fluid-structure interaction (FSI) approach that combines the lattice Boltzmann method (LBM) and finite element (FE) method. Two support modes of LVADs, namely constant speed mode and constant flow mode, which have been widely applied to clinical practice, are also designed. Results demonstrate that the support modes of LVADs could significantly affect the biomechanical states of the aortic valve and the blood flow pattern of the ascending aorta. Compared with those in the constant flow mode, the leaflets in the constant speed mode could achieve better dynamic performance and lower stress during the systolic phase. The max radial displacement of the leaflets in the constant speed mode is at 8 mm, whereas that in the constant flow mode is at 0.8 mm. Furthermore, the outflow of LVADs directly impacts the aortic surfaces of the leaflets during the diastolic phase by increasing the level of wall shear stress of the leaflets. The leaflets in the constant speed mode receive less impact than those in the constant flow mode. The condition with such minimal impact is conducive to maintaining the normal structure of leaflets and benefits the reduction of the risk of valvular diseases. In sum, the support modes of LVADs exert a crucial effect on the biomechanical environment of the aortic valve. The constant speed mode is better than the constant flow mode in terms of providing a good hemodynamic environment for the aortic valve.

Comparison of a fixed-grid and arbitrary Lagrangian–Eulerian methods on modelling fluid–structure interaction of the aortic valve

This study was aimed at assessing the robustness of a fixed-grid fluid–structure interaction method (Multi-Material Arbitrary Lagrangian–Eulerian) to modelling the two-dimensional native aortic valve dynamics and comparing it to the Arbitrary Lagrangian–Eulerian method. For the fixed-grid method, the explicit finite element solver LS-DYNA was utilized, where two independent meshes for the fluid and structure were generated and the penalty method was used to handle the coupling between the fluid and structure domains. For the Arbitrary Lagrangian–Eulerian method, the implicit finite element solver ADINA was used where two separate conforming meshes were used for the valve structure and the fluid domains. The comparison demonstrated that both fluid–structure interaction methods predicted accurately the valve dynamics, fluid flow, and stress distribution, implying that fixed-grid methods can be used in situations where the Arbitrary Lagrangian–Eulerian method fails.

Validation and Extension of a Fluid-Structure Interaction Model of the Healthy Aortic Valve

Purpose

The understanding of the optimum function of the healthy aortic valve is essential in interpreting the effect of pathologies in the region, and in devising effective treatments to restore the physiological functions. Still, there is no consensus on the operating mechanism that regulates the valve opening and closing dynamics. The aim of this study is to develop a numerical model that can support a better comprehension of the valve function and serve as a reference to identify the changes produced by specific pathologies and treatments.

Methods

A numerical model was developed and adapted to accurately replicate the conditions of a previous in vitro investigation into aortic valve dynamics, performed by means of particle image velocimetry (PIV). The resulting velocity fields of the two analyses were qualitatively and quantitatively compared to validate the numerical model. In order to simulate more physiological operating conditions, this was then modified to overcome the main limitations of the experimental setup, such as the presence of a supporting stent and the non-physiological properties of the fluid and vessels.

Results

The velocity fields of the initial model resulted in good agreement with those obtained from the PIV, with similar flow structures and about 90% of the computed velocities after valve opening within the standard deviation of the equivalent velocity measurements of the in vitro model. Once the experimental limitations were removed from the model, the valve opening dynamics changed substantially, with the leaflets opening into the sinuses to a much greater extent, enlarging the effective orifice area by 11%, and reducing greatly the vortical structures previously observed in proximity of the Valsalva sinuses wall.

Conclusions

The study suggests a new operating mechanism for the healthy aortic valve leaflets considerably different from what reported in the literature to date and largely more efficient in terms of hydrodynamic performance. This work also confirms the crucial role that numerical approaches, complemented with experimental findings, can play in overcoming some of the limitations inherent in experimental techniques, supporting the full understanding of complex physiological phenomena.

Electronic supplementary material

The online version of this article (doi:10.1007/s13239-018-00391-1) contains supplementary material, which is available to authorized users.

Fluid-structure interaction simulation of artificial textile reinforced aortic heart valve: Validation with an in-vitro test

Prosthetic heart valves deployed in the left heart (aortic and mitral) are subjected to harsh hemodynamical conditions. Most of the tissue engineered heart valves have been developed for the low pressure pulmonary position because of the difficulties in fabricating a mechanically strong valve, able to withstand the systemic circulation. This necessitates the use of reinforcing scaffolds, resulting in a tissue-engineered textile reinforced tubular aortic heart valve. Therefore, to better design these implants, material behaviour of the composite, valve kinematics and its hemodynamical response need to be evaluated. Experimental assessment can be immensely time consuming and expensive, paving way for numerical studies. In this work, the material properties obtained using the previously proposed multi-scale numerical method for textile composites was evaluated for its accuracy. An in silico immersed boundary (IB) fluid structure interaction (FSI) simulation emulating the in vitro experiment was set-up to evaluate and compare the geometric orifice area and flow rate for one beat cycle. Results from the in silico FSI simulation were found to be in good coherence with the in vitro test during the systolic phase, while mean deviation of approximately 9% was observed during the diastolic phase of a beat cycle. Merits and demerits of the in silico IB-FSI method for the presented case study has been discussed with the advantages outweighing the drawbacks, indicating the potential towards an effective use of this framework in the development and analysis of heart valves.

Comparative Fluid-Structure Interaction Analysis of Polymeric Transcatheter and Surgical Aortic Valves' Hemodynamics and Structural Mechanics

Transcatheter aortic valve replacement (TAVR) has emerged as an effective alternative to conventional surgical aortic valve replacement (SAVR) in high-risk elderly patients with calcified aortic valve disease. All currently FDA-approved TAVR devices use tissue valves that were adapted to but not specifically designed for TAVR use. Emerging clinical evidence indicates that these valves may get damaged during crimping and deployment- leading to valvular calcification, thrombotic complications, and limited durability. This impedes the expected expansion of TAVR to lower-risk and younger patients. Viable polymeric valves have the potential to overcome such limitations. We have developed a polymeric SAVR valve, which was optimized to reduce leaflet stresses and offer a thromboresistance profile similar to that of a tissue valve. This study compares the polymeric SAVR valve's hemodynamic performance and mechanical stresses to a new version of the valve- specifically designed for TAVR. Fluid-structure interaction (FSI) models were utilized and the valves' hemodynamics, flexural stresses, strains, orifice area, and wall shear stresses were compared. The TAVR valve had 42% larger opening area and 27% higher flow rate versus the SAVR valve, while wall shear stress distribution and mechanical stress magnitudes were of the same order, demonstrating the enhanced performance of the TAVR valve prototype. The TAVR valve FSI simulation and Vivitro pulse duplicator experiments were compared in terms of the leaflets' kinematics and the effective orifice area. The numerical methodology presented can be further used as a predictive tool for valve design optimization for enhanced hemodynamics and durability.

The Effects of Left Ventricular Assist Device Support Level on the Biomechanical States of Aortic Valve

Background

Although aortic valve disease caused by left ventricular assist device (LVAD) support has attracted more and more attention, the precise biomechanical effects of LVAD support level on the aortic valve are still unclear.

Material/Methods

A structural finite element models study was conducted using an ideal aortic valve geometric model. Four different study conditions were designed, according to the reduction of the open duration of the aortic valve. The isotropic hyperelastic constitutive equation was chosen to reflect the mechanical property of the leaflets. The distribution of the stress, strain, and transient dynamics of the leaflet were calculated.

Results

Along with the increase of LVAD support level, the open duration of the aortic valve was also reduced by the increase of LVAD support (low support level case 0.23 seconds versus middle support level case 0.2 seconds versus high support level case 0.14 seconds). Moreover, along with the increase of support mode of LVAD, the von Mises stress in most leaflet areas was increased from the low stress level (0–0.4 MPa) to the middle region (0.4–0.8 MPa). Once the leaflets were continuously closed, the high stress level (larger than 0.8 MPa) was observed. In contrast, the support level of LVAD only had slight effects on the distribution of von Mises strain. According to the aforementioned results, maintaining the open duration of aortic valve longer than 0.2 seconds could achieve better performance of biomechanical states of leaflets.

Conclusions

This study could provide useful information on the determination of optimal LVAD support strategy.

Computational methods for the aortic heart valve and its replacements

INTRODUCTION

Replacement with a prosthetic device remains a major treatment option for the patients suffering from heart valve disease, with prevalence growing resulting from an ageing population. While the most popular replacement heart valve continues to be the bioprosthetic heart valve (BHV), its durability remains limited. There is thus a continued need to develop a general understanding of the underlying mechanisms limiting BHV durability to facilitate development of a more durable prosthesis. In this regard, computational models can play a pivotal role as they can evaluate our understanding of the underlying mechanisms and be used to optimize designs that may not always be intuitive. Areas covered: This review covers recent progress in computational models for the simulation of BHV, with a focus on aortic valve (AV) replacement. Recent contributions in valve geometry, leaflet material models, novel methods for numerical simulation, and applications to BHV optimization are discussed. This information should serve not only to infer reliable and dependable BHV function, but also to establish guidelines and insight for the design of future prosthetic valves by analyzing the influence of design, hemodynamics and tissue mechanics. Expert commentary: The paradigm of predictive modeling of heart valve prosthesis are becoming a reality which can simultaneously improve clinical outcomes and reduce costs. It can also lead to patient-specific valve design.

3D physiological model of the aortic valve incorporating small coronary arteries

The diseases of the coronary arteries and the aortic root are still the leading causes of mortality and morbidity worldwide. In this study, a 3D global fluid-structure interaction of the aortic root with inclusion of anatomically inspired small coronary arteries using the finite element method is presented. This innovative model allows to study the impact and interaction of root biomechanics on coronary hemodynamics and brings a new understanding to small coronary vessels hemodynamics. For the first time, the velocity profiles and shear stresses are reported in distal coronary arteries as a result of the aortic flow conditions in a global fluid-structure interaction model.

Effects of Biowastes Released by Mechanically Damaged Muscle Cells on the Propagation of Deep Tissue Injury: A Multiphysics Study

Deep tissue injuries occur in muscle tissues around bony prominences under mechanical loading leading to severe pressure ulcers. Tissue compression can potentially compromise lymphatic transport and cause accumulation of metabolic biowastes, which may cause further cell damage under continuous mechanical loading. In this study, we hypothesized that biowastes released by mechanically damaged muscle cells could be toxic to the surrounding muscle cells and could compromise the capability of the surrounding muscle cells to withstand further mechanical loadings. In vitro, we applied prolonged low compressive stress (PLCS) and short-term high compressive stress to myoblasts to cause cell damage and collected the biowastes released by the damaged cells under the respective loading scenarios. In silico, we used COMSOL to simulate the compressive stress distribution and the diffusion of biowastes in a semi-3D buttock finite element model. In vitro results showed that biowastes collected from cells damaged under PLCS were more toxic and could compromise the capability of normal myoblasts to resist compressive damage. In silico results showed that higher biowastes diffusion coefficient, higher biowastes release rate, lower biowastes tolerance threshold and earlier timeline of releasing biowastes would cause faster propagation of tissue damage. This study highlighted the importance of biowastes in the development of deep tissue injury to clinical pressure ulcers under prolonged skeletal compression.