Finite Volume Method
Dotaz
Zobrazit nápovědu
BACKGROUND: Total knee arthroplasty (TKA) with all-polyethylene tibial (APT) components has shown comparable survivorship and clinical outcomes to that with metal-backed tibial (MBT). Although MBT is more frequently implanted, APT equivalents are considered a low-cost variant for elderly patients. A biomechanical analysis was assumed to be suitable to compare the response of the periprosthetic tibia after implantation of TKA NexGen APT and MBT equivalent. METHODS: A standardised load model was used representing the highest load achieved during level walking. The geometry and material models were created using computed tomography data. In the analysis, a material model was created that represents a patient with osteopenia. RESULTS: The equivalent strain distribution in the models of cancellous bone with an APT component showed values above 1000 με in the area below the medial tibial section, with MBT component were primarily localised in the stem tip area. For APT variants, the microstrain values in more than 80% of the volume were in the range from 300 to 1500 με, MBT only in less than 64% of the volume. CONCLUSION: The effect of APT implantation on the periprosthetic tibia was shown as equal or even superior to that of MBT despite maximum strain values occurring in different locations. On the basis of the strain distribution, the state of the bone tissue was analysed to determine whether bone tissue remodelling or remodelling would occur. Following clinical validation, outcomes could eventually modify the implant selection criteria and lead to more frequent implantation of APT components.
- Klíčová slova
- All-polyethylene tibial component, Computational modeling, FEA, Finite element method, Knee replacement, Metal-backed tibial component, TKR, Total knee arthroplasty,
- MeSH
- analýza metodou konečných prvků MeSH
- kovy MeSH
- lidé MeSH
- polyethylen MeSH
- protézy - design MeSH
- protézy kolene * MeSH
- senioři MeSH
- tibie diagnostické zobrazování chirurgie MeSH
- totální endoprotéza kolene * MeSH
- Check Tag
- lidé MeSH
- senioři MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- kovy MeSH
- polyethylen MeSH
BACKGROUND AND OBJECTIVE: Total knee arthroplasty (TKA) with modern all-polyethylene tibial (APT) components has shown high long-term survival rates and comparable results to those with metal-backed tibial components. Nevertheless, APT components are primarily recommended for older and low-demand patients. There are no evidence-based biomechanical guidelines for orthopaedic surgeons to determine the appropriate lower age limit for implantation of APT components. A biomechanical analysis was assumed to be suitable to evaluate the clinical results in patients under 70 years. The scope of this study was to determine biomechanically the appropriate lower age limit for implantation of APT components. METHODS: To generate data of the highest possible quality, the geometry of the computational models was created based on computed tomography (CT) images of a representative patient. The cortical bone tissue model distinguishes the change in mechanical properties described in three parts from the tibial cut. The cancellous bone material model has a heterogeneous distribution of mechanical properties. The values used to determine the material properties of the tissues were obtained from measurements of a CT dataset comprising 45 patients. RESULTS: Computational modeling showed that in the majority of the periprosthetic volume, the von Mises strain equivalent ranges from 200 to 2700 με; these strain values induce bone modeling and remodeling. The highest measured deformation value was 2910 με. There was no significant difference in the induced mechanical response between bone models of the 60-year and 70-year age groups, and there was <3% difference from the 65-year age group. CONCLUSIONS: Considering in silico limitations, we suggest that APT components could be conveniently used on a bone with mechanical properties of the examined age categories. Under defined loading conditions, implantation of TKA with APT components is expected to induce modeling and remodeling of the periprosthetic tibia. Following clinical validation, the results of our study could modify the indication criteria of the procedure, and lead to more frequent implantation of all-polyethylene TKA in younger patients.
- Klíčová slova
- All-polyethylene tibial component, Bone mechanics, Finite element model, Orthopaedic biomechanics, Total knee arthroplasty, von Mises strain,
- MeSH
- analýza metodou konečných prvků MeSH
- biomechanika MeSH
- kovy MeSH
- lidé MeSH
- mechanický stres MeSH
- polyethylen MeSH
- protézy - design MeSH
- protézy kolene * MeSH
- tibie diagnostické zobrazování chirurgie MeSH
- totální endoprotéza kolene * metody MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- kovy MeSH
- polyethylen MeSH
We have investigated different approaches to handling parallel-tempering Monte Carlo (PTMC) simulations in the isothermal-isobaric ensemble of molecular cluster/nanoparticle systems for predicting structural phase diagram transitions. We have implemented various methodologies that consist of treating pressure implicitly through its effect on the volume. Thus, the main problem in the simulations under nonzero pressure becomes the volume definition of the finite nonperiodic system, and we considered approaches based on the particles' coordinates. Various volume models, namely container-volume, particle-volume, average-volume, ellipsoids-volume, and convex hull-volume, were employed, and the required corrections for each of them in the Monte Carlo computations were introduced. Finally, we explored the effects of volume/pressure changes for all models on structural phase transitions of a test system, such as the small "icelike" (H2O)12 water cluster. The temperature and pressure dependence of the cluster's heat capacity and energy-volume Pearson correlation coefficient were studied, phase diagrams were constructed using a multiple-histogram method, and attempts were made to identify phase transitions to particular cluster structures. Our results show significant differences between the employed volume models, and we discuss all pressure-induced, such as solid-solid-, solid-liquid-, and liquid-gas-like, phase transformations in the present study.
- Publikační typ
- časopisecké články MeSH
We present comprehensive biomechanical analyses of abdominal aortic aneurysms (AAA) for 43 patients. We compare stress magnitudes and stress distributions within arterial walls of abdominal aortic aneurysms (AAA) obtained using two simulation and modelling methods: (a) Fully automated and computationally very efficient linear method embedded in the software platform Biomechanics based Prediction of Aneurysm Rupture Risk (BioPARR), freely available from https://bioparr.mech.uwa.edu.au/; (b) More complex and much more computationally demanding Non-Linear Iterative Stress Analysis (Non-LISA) that uses a non-linear inverse iterative approach and strongly non-linear material model. Both methods predicted localised high stress zones with over 90% of AAA model volume fraction subjected to stress below 20% of the 99th percentile maximum principal stress. However, for the non-linear iterative method, the peak maximum principal stress (and 99th percentile maximum principal stress) was higher and the stress magnitude in the low stress area lower than for the automated linear method embedded in BioPARR. Differences between the stress distributions obtained using the two methods tended to be particularly pronounced in the areas where the AAA curvature was large. Performance of the selected characteristic features of the stress fields (we used 99th percentile maximum principal stress) obtained using BioPARR and Non-LISA in distinguishing between the AAAs that would rupture and remain intact was for practical purposes the same for both methods.
- Klíčová slova
- abdominal aortic aneurysm, finite element method, patient-specific modelling, rupture, stress,
- MeSH
- analýza metodou konečných prvků MeSH
- aneurysma břišní aorty * MeSH
- aorta abdominalis MeSH
- biomechanika MeSH
- lidé MeSH
- mechanický stres MeSH
- modely kardiovaskulární MeSH
- ruptura aorty * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
PURPOSE: The cement augmentation of a conventional anterior screw fixation in type II odontoid process fractures for elderly patients significantly increased stiffness and load to failure under anterior-posterior load in comparison with non-augmented fixation. The amount and quality of bone cement are usually taken ad hoc in clinical practise. In this study, we wanted to clarify the role of bone cement amount and its quality to the stiffness of odontoid and vertebrae body junction. METHODS: Finite-element method was used to achieve different scenarios of cement augmentation. For all models, an initial stiffness was calculated. Model (1) the intact vertebrae were virtually potted into a polymethylmethacrylate base via the posterior vertebral arches. A V-shaped punch was used for loading the odontoid in an anterior-posterior direction. (2) The odontoid fracture type IIa (Anderson-D'Alonzo classification) was achieved by virtual transverse osteotomy. Anterior screw fixation was virtually performed by putting self-drilling titanium alloy 3.5 mm diameter anterior cannulated lag screw with a 12 mm thread into the inspected vertebrae. A V-shaped punch was used for loading the odontoid in an anterior-posterior direction. The vertebrae body was assumed to be non-cemented and cemented with different volume. RESULTS: The mean cement volume was lowest for body base filling with 0.47 ± 0.03 ml. The standard body filling corresponds to 0.95 ± 0.15 ml. The largest volume corresponds to 1.62 ± 0.12 ml in the presence of cement leakage. The initial stiffness of the intact C2 vertebrae was taken as the reference value. The mean initial stiffness for non-porous cement (E = 3000 MPa) increased linearly (R2 = 0.98). The lowest stiffness (123.3 ± 5.8 N/mm) was measured in the intact C2 vertebrae. However, the highest stiffness (165.2 ± 5.2 N/mm) was measured when cement leakage out of the odontoid peg occurred. The mean initial stiffness of the base-only cemented group was 147.2 ± 8.4 N/mm compared with 157.9 ± 6.6 N/mm for the base and body cemented group. This difference was statistically significant (p < 0.0061). The mean initial stiffness for porous cement (E = 500 MPa) remains constant. Therefore, there is no difference between cemented and non-cemented junction. This difference was not statistically significant (p < 0.18). CONCLUSION: The present study showed that the low porous cement was able to significantly influence the stiffness of the augmented odontoid screw fixation in vitro, although further in vivo clinical studies should be undertaken. Our results suggest that only a small amount of non-porous cement is needed to restore stiffness at least to its pre-fracture level and this can be achieved with the injection of 0.7-1.2 ml of cement. These slides can be retrieved under Electronic Supplementary Material.
- Klíčová slova
- Bone cement, Finite element, Fracture, Odontoid, Spine,
- MeSH
- dens axis * diagnostické zobrazování zranění chirurgie MeSH
- fraktury páteře * diagnostické zobrazování chirurgie MeSH
- kostní cementy terapeutické užití MeSH
- kostní šrouby MeSH
- lidé MeSH
- senioři MeSH
- vnitřní fixace fraktury MeSH
- Check Tag
- lidé MeSH
- senioři MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- kostní cementy MeSH
We determined the total system elastic Helmholtz free energy, under the constraints of constant temperature and volume, for systems comprised of one or more perfectly bonded hard spherical inclusions (i.e. "hard spheres") embedded in a finite spherical elastic solid. Dirichlet boundary conditions were applied both at the surface(s) of the hard spheres, and at the outer surface of the elastic solid. The boundary conditions at the surface of the spheres were used to describe the rigid displacements of the spheres, relative to their initial location(s) in the unstressed initial state. These displacements, together with the initial positions, provided the final shape of the strained elastic solid. The boundary conditions at the outer surface of the elastic medium were used to ensure constancy of the system volume. We determined the strain and stress tensors numerically, using a method that combines the Neuber-Papkovich spherical harmonic decomposition, the Schwartz alternating method, and Least-squares for determining the spherical harmonic expansion coefficients. The total system elastic Helmholtz free energy was determined by numerically integrating the elastic Helmholtz free energy density over the volume of the elastic solid, either by a quadrature, or a Monte Carlo method, or both. Depending on the initial position of the hard sphere(s) (or equivalently, the shape of the un-deformed stress-free elastic solid), and the displacements, either stationary or non-stationary Helmholtz free energy minima were found. The non-stationary minima, which involved the hard spheres nearly in contact with one another, corresponded to lower Helmholtz free energies, than did the stationary minima, for which the hard spheres were further away from one another.
- Klíčová slova
- Flowing Matter: Liquids and Complex Fluids,
- Publikační typ
- časopisecké články MeSH
The numerical simulation of inhaled aerosols in medical research starts to play a crucial role in understanding local deposition within the respiratory tract, a feat often unattainable experimentally. Research on children is particularly challenging due to the limited availability of in vivo data and the inherent morphological intricacies. CFD solvers based on Finite Volume Methods (FVM) have been widely employed to solve the flow field in such studies. Recently, Lattice Boltzmann Methods (LBM), a mesoscopic approach, have gained prominence, especially for their scalability on High-Performance Computers. This study endeavours to compare the effectiveness of LBM and FVM in simulating particulate flows within a child's respiratory tract, supporting research related to particle deposition and medication delivery using LBM. Considering a 5-year-old child's airway model at a steady inspiratory flow, the results are compared with in vitro experiments. Notably, both LBM and FVM exhibit favourable agreement with experimental data for the mean velocity field and the turbulence intensity. For particle deposition, both numerical methods yield comparable results, aligning well with in vitro experiments across a particle size range of 0.1-20 µm. Discrepancies are identified in the upper airways and trachea, indicating a lower deposition fraction than in the experiment. Nonetheless, both LBM and FVM offer invaluable insights into particle behaviour for different sizes, which are not easily achievable experimentally. In terms of practical implications, the findings of this study hold significance for respiratory medicine and drug delivery systems - potential health impacts, targeted drug delivery strategies or optimisation of respiratory therapies.
- Klíčová slova
- Child airways, Finite Volume Method, In vitro measurement, Lattice Boltzmann Method, Particle deposition,
- MeSH
- aerosoly MeSH
- hydrodynamika * MeSH
- lidé MeSH
- počítačová simulace MeSH
- předškolní dítě MeSH
- trachea * anatomie a histologie MeSH
- velikost částic MeSH
- Check Tag
- lidé MeSH
- předškolní dítě MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- aerosoly MeSH
Relatively few polymer gels are known to show volume phase transition where the gels undergo an abrupt change in the degree of swelling by passing through a three-phase equilibrium. Characteristic for such transition is the existence of van der Waals (vdW) loop on the dependence of solvent chemical potential versus polymer concentration. For the χ-induced transition, the existence of vdW loop is determined by the concentration dependence of the interaction function. It is shown that expansive mechanical strains can assist in development of the vdW loop. Systems characterized by continuous change of the degree of swelling transform upon such strain into ones where the degree of swelling changes much and abruptly. Also, expansive modes of strain can make the transition wider and more robust in gels where transition is already observed under free swelling condition. The possibility to induce the volume phase transition by external stresses can be utilized for finding other stimuli sensitive gels, strengthening of gel response, and in modeling of properties of gel constructs by Finite Element Method.
- Publikační typ
- časopisecké články MeSH
A fast and accurate fusion of intra-operative images with a pre-operative data is a key component of computer-aided interventions which aim at improving the outcomes of the intervention while reducing the patient's discomfort. In this paper, we focus on the problematic of the intra-operative navigation during abdominal surgery, which requires an accurate registration of tissues undergoing large deformations. Such a scenario occurs in the case of partial hepatectomy: to facilitate the access to the pathology, e.g. a tumor located in the posterior part of the right lobe, the surgery is performed on a patient in lateral position. Due to the change in patient's position, the resection plan based on the pre-operative CT scan acquired in the supine position must be updated to account for the deformations. We suppose that an imaging modality, such as the cone-beam CT, provides the information about the intra-operative shape of an organ, however, due to the reduced radiation dose and contrast, the actual locations of the internal structures necessary to update the planning are not available. To this end, we propose a method allowing for fast registration of the pre-operative data represented by a detailed 3D model of the liver and its internal structure and the actual configuration given by the organ surface extracted from the intra-operative image. The algorithm behind the method combines the iterative closest point technique with a biomechanical model based on a co-rotational formulation of linear elasticity which accounts for large deformations of the tissue. The performance, robustness and accuracy of the method is quantitatively assessed on a control semi-synthetic dataset with known ground truth and a real dataset composed of nine pairs of abdominal CT scans acquired in supine and flank positions. It is shown that the proposed surface-matching method is capable of reducing the target registration error evaluated of the internal structures of the organ from more than 40 mm to less then 10 mm. Moreover, the control data is used to demonstrate the compatibility of the method with intra-operative clinical scenario, while the real datasets are utilized to study the impact of parametrization on the accuracy of the method. The method is also compared to a state-of-the art intensity-based registration technique in terms of accuracy and performance.
- Klíčová slova
- Computer-aided interventions, Deformable image registration, Finite element method, Human liver,
- MeSH
- algoritmy MeSH
- analýza metodou konečných prvků MeSH
- biomechanika MeSH
- břicho diagnostické zobrazování chirurgie MeSH
- chirurgie s pomocí počítače metody MeSH
- elastografie * MeSH
- lidé MeSH
- nemoci jater diagnostické zobrazování chirurgie MeSH
- peroperační doba MeSH
- počítačová rentgenová tomografie * MeSH
- počítačová simulace MeSH
- počítačová tomografie s kuželovým svazkem * MeSH
- polohování pacienta MeSH
- rentgenový obraz - interpretace počítačová metody MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- validační studie MeSH
Considering the fact that hemodynamics plays an important role in the patency and overall performance of implanted bypass grafts, this work presents a numerical investigation of pulsatile non-Newtonian blood flow in three different patient-specific aorto-coronary bypasses. The three bypass models are distinguished from each other by the number of distal side-to-side and end-to-side anastomoses and denoted as single, double and triple bypasses. The mathematical model in the form of time-dependent nonlinear system of incompressible Navier-Stokes equations is coupled with the Carreau-Yasuda model describing the shear-thinning property of human blood and numerically solved using the principle of the SIMPLE algorithm and cell-centred finite volume method formulated for hybrid unstructured tetrahedral grids. The numerical results computed for non-Newtonian and Newtonian blood flow in the three aorto-coronary bypasses are compared and analysed with emphasis placed on the distribution of cycle-averaged wall shear stress and oscillatory shear index. As shown in this study, the non-Newtonian blood flow in all of the considered bypass models does not significantly differ from the Newtonian one. Our observations further suggest that, especially in the case of sequential grafts, the resulting flow field and shear stimulation are strongly influenced by the diameter of the vessels involved in the bypassing.
- Klíčová slova
- Carreau-Yasuda model, SIMPLE algorithm, finite volume method, hybrid grid system, patient-specific aorto-coronary bypass, pulsatile blood flow,
- MeSH
- algoritmy MeSH
- hemodynamika fyziologie MeSH
- koronární bypass metody MeSH
- koronární cévy chirurgie MeSH
- lidé MeSH
- mechanický stres MeSH
- rychlost toku krve fyziologie MeSH
- teoretické modely MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH