An investigation of the effect of brain atrophy on brain injury in multiple sclerosis.

Multiple sclerosis (MS) is a disease of the central nervous system (CNS) that affects the brain and spinal cord. It is estimated that the average prevalence of MS is 35.9 cases per 100,000 and a total of 2.8 million people worldwide have MS. Brain atrophy is usually seen in the early stages of MS, and its progress is faster than healthy people. The present study was a numerical study that uses the Fluid-structure interaction (FSI) model to investigate the effect of brain atrophy on brain injury in MS. Firstly, a healthy model was constructed from MRI images and validated by experimental data. Then three models with different degrees of brain atrophy, which showed the rate of brain atrophy in different years in MS patients, were developed to model the brain atrophy in MS. The models were subjected to two different types of impact conditions. Type I, which only produced a translational motion and the HIC value of 744, was applied to each model. Type II produced both translational and rotational motion. In this type of impact, the experimental kinematics, with peaks of 450 g for the translational acceleration and 26.2 krad/s2 for the rotational acceleration, were applied to the nodes that located in the center of gravity of the head models and the results were extracted from each one. According to the results of impact type I, the pressure of the frontal lobe of the brain is 149,647 Pa in the health model and 137,690 Pa in the model with severe atrophy.

A mathematical model for biomechanical behavior of the aortic arch.

The aortic arch plays a significant role in homeostatic mechanisms to retain blood pressure at stable balance in the cardiovascular system. Therefore, the objective is to estimate and identify cardiovascular illness imposed by the abnormal blood hemodynamic domain. In this regard, hemodynamic forces are monitored by the baroreflex of the artery wall. Therefore, these receptors quickly detect the abnormal stress magnitudes in the aortic arterial wall. The present study presents a 3D aortic arch model extracted by a Computerized tomography scan. Also, the numerical solution was carried out by ANSYS 2020 R1 in view of Fluid-Structure Interaction After that, we found wall shear stress (WSS), pressure, and velocity in the fluid domain. Also, the normal stress was analyzed to determine the aortic arch baroreflex location in the solid range. In this regard, higher WSS values are measured at the supra-aortic branches going out the aortic arch that reached 42.5 Pa. Also, higher normal stress happened at the aortic root and the supra-aortic branches and reached approximately 200 kPa at peak systole.

Modeling of aortic valve stenosis using fluid-structure interaction method.

BACKGROUND AND OBJECTIVES: Fluid structure interaction (FSI) is defined as interaction of the structures with contacting fluids. The aortic valve experiences the interaction with blood flow in systolic phase. In this study, we have tried to predict the hemodynamics of blood flow through a normal and stenotic aortic valve in two relaxation and exercise conditions using a three-dimensional FSI method. METHODS: The aorta valve was modeled as a three-dimensional geometry including a normal model and two others with 25% and 50% stenosis. The geometry of the aortic valve was extracted from CT images and the models were generated by MMIMCS software and then they were implemented in ANSYS software. The pulsatile flow rate was used for all cases and the numerical simulations were conducted based on a time-dependent domain. RESULTS: The obtained results including the velocity, pressure, and shear stress contours in different systolic time sequences were explained and discussed. The maximum blood flow velocity in relaxation phase was obtained 1.62 m/s (normal valve), 3.78 m/s (25% stenosed valve), and 4.73 m/s (50% stenosed valve). In exercise condition, the maximum velocities are 2.86, 4.32, and 5.42 m/s respectively. The maximum blood pressure in relaxation phase was calculated 111.45 mmHg (normal), 148.66 mmHg (25% stenosed), and 164.21 mmHg (50% stenosed). However, the calculated values in exercise situation were 129.57, 163.58, and 191.26 mmHg.The validation of the predicted results was also conducted using existing literature. CONCLUSIONS: We believe that such model are useful tools for biomechanical experts. The further studies should be done using experimental data and the data are implemented on the boundary conditions for better comparison of the results.

A New Template and Teleoperation System for Human-Guided Spine Surgery.

Fluoroscopy-guided placement of pedicle screws is usually used to position pedicle screws, although it is highly risky due to lack of accuracy during operation and surgeon's intensive labor requirement during surgery. In this study, a new method was introduced to consider the issues and we tried to reduce the risk of pedicle screw placement and drilling process. To put pedicle screws in the correct position and orientation, a specific drill guide is designed and fabricated by additive manufacturing technology. In addition, since the drilling process is remarkably important and it is usually dependent on surgeon skill, therefore a teleoperation system is proposed to perform this task. In order to let the surgeon to have better control on the patient, a control scheme including position and velocity signals along with surgeon force and reaction force of vertebra was proposed. This helped the surgeon for proper control on the patient during surgery. A force estimation algorithm was presented to eliminate measuring external force signal. Consequently, 10 vertebras were used to evaluate specific drill guide and teleoperation system simultaneously. Then, the computed tomography evaluation demonstrated that Kirschner wire trajectories followed the planned axes with an error between 1 and 2 mm in 20% of cases and less than 1 mm in 80% of cases. The results obtained by the teleoperation system showed that the surgeon, in the master side, has proper control over the patient body in the slave side because the slave robot followed the master position. Furthermore, the surgeon in the master side sensed the reaction force of vertebra in the slave side appropriately.

Development of a fluid-structure interaction model to simulate mitral valve malcoaptation.

OBJECT: Mitral regurgitation (MR) is a condition in which the mitral valve does not prevent the reversal of blood flow from the left ventricle into the left atrium. This study aimed at numerically developing a model to mimic MR and poor leaflet coaptation and also comparing the performance of a normal mitral valve to that of the MR conditions at different gap junctions of 1, 3 and 5 mm between the anterior and posterior leaflets. RESULTS: The results revealed no blood flow to the left ventricle when a gap between the leaflets was 0 mm. However, MR increased this blood flow, with increases in the velocity and pressure within the atrium. However, the pressure within the aorta did not vary meaningfully (ranging from 22 kPa for a 'healthy' model to 25 kPa for severe MR). CONCLUSIONS: The findings from this study have implications not only for understanding the changes in pressure and velocity as a result of MR in the ventricle, atrium or aorta, but also for the development of a computational model suitable for clinical translation when diagnosing and determining treatment for MR.

Color spectrographic respiratory monitoring from the external ear canal.

The need for simple and reliable means of respiratory monitoring has existed since the beginnings of medicine. In the present study, we describe the use of color spectrographic analysis of breathing sounds recorded from the external ear canal as a candidate technology to meet this need. A miniature electret microphone was modified with the addition of an adapter to allow it to be placed comfortably in the external ear canal. The amplified signal was then connected to a real-time color spectrogram program running on a laptop personal computer utilizing the Windows operating system. Based on the results obtained, we hypothesize that the real-time display of color spectrogram breathing patterns locally or at a central monitoring station may turn out to be a useful means of respiratory monitoring in patients at increased risk of respiratory depression or other respiratory problems. Finally, we conducted a statistical analysis that suggests that significant spectrogram differences may exist among some groups investigated in the study.

A medium invasiveness multi-level patient's specific template for pedicle screw placement in the scoliosis surgery.

BACKGROUND: Several methods including free-hand technique, fluoroscopic guidance, image-guided navigation, computer-assisted surgery system, robotic platform and patient's specific templates are being used for pedicle screw placement. These methods have screw misplacements and are not always easy to be applied. Furthermore, it is necessary to expose completely a large portions of the spine in order to access fit entirely around the vertebrae. METHODS: In this study, a multi-level patient's specific template with medium invasiveness was proposed for pedicle screw placement in the scoliosis surgery. It helps to solve the problems related to the soft tissues removal. After a computer tomography (CT) scan of the spine, the templates were designed based on surgical considerations. Each template was manufactured using three-dimensional printing technology under a semi-flexible post processing. The templates were placed on vertebras at four points-at the base of the superior-inferior articular processes on both left-right sides. This helps to obtain less invasive and more accurate procedure as well as true-stable and easy placement in a unique position. The accuracy of screw positions was confirmed by CT scan after screw placement. RESULTS: The result showed the correct alignment in pedicle screw placement. In addition, the template has been initially tested on a metal wire series Moulage (height 70 cm and material is PVC). The results demonstrated that it could be possible to implement it on a real patient. CONCLUSIONS: The proposed template significantly reduced screw misplacements, increased stability, and decreased the sliding & the intervention invasiveness.

Analysis of central venous pressure (CVP) signals using mathematical methods.

Central venous pressure (CVP) is an important clinical parameter for physicians but only the absolute CVP value is typically monitored in the intensive care unit (ICU). In this study, we propose a novel mathematical method to present and analyze CVP signals. A total of 44 suitable samples were chosen from a total of 65 collected in an ICU. Pre-processing of the samples included rate reduction and digital filtering. The statistical features of time and frequency domain, wavelet, and empirical mode decomposition of these signals were extracted. We found no significant difference among the CVP signals regarding sex, smoking, coronary disease, and respiration mode of the samples.

A multi-channel acoustics monitor for perioperative respiratory monitoring: preliminary data.

This study pertains to a six-channel acoustic monitoring system for use in patient monitoring during or after surgery. The base hardware consists of a USB data acquisition system, a custom-built six-channel amplification system, and a series of microphones of various designs. The software is based on the MATLAB platform with data acquisition drivers installed. The displayed information includes: time domain signals, frequency domain signals, and tools to aid in the detection of endobronchial intubation. We hypothesize that the above mentioned arrangement may be helpful to the anesthesiologist in recognizing clinical conditions like wheezing, bronchospasm, endobronchial intubation, and apnea. The study also evaluated various types of microphone designs used to transduce breath sounds. The system also features selectable band-pass filtering using MATLAB algorithms as well as a collection of recordings obtained with the system to establish what respiratory acoustic signals look like under various conditions.

Estimation of maximum intraventricular pressure: a three-dimensional fluid-structure interaction model.

BACKGROUND: The aim of this study was to propose a method to estimate the maximum pressure in the left ventricle (MPLV) for a healthy subject, based on cardiac outputs measured by echo-Doppler (non-invasive) and catheterization (invasive) techniques at rest and during exercise. METHODS: Blood flow through aortic valve was measured by Doppler flow echocardiography. Aortic valve geometry was calculated by echocardiographic imaging. A Fluid-structure Interaction (FSI) simulation was performed, using an Arbitrary Lagrangian-Eulerian (ALE) mesh. Boundary conditions were defined as pressure loads on ventricular and aortic sides during ejection phase. The FSI simulation was used to determine a numerical relationship between the cardiac output to aortic diastolic and left ventricular pressures. This relationship enabled the prediction of pressure loads from cardiac outputs measured by invasive and non-invasive clinical methods. RESULTS: Ventricular systolic pressure peak was calculated from cardiac output of Doppler, Fick oximetric and Thermodilution methods leading to a 22%, 18% and 24% increment throughout exercise, respectively. The mean gradients obtained from curves of ventricular systolic pressure based on Doppler, Fick oximetric and Thermodilution methods were 0.48, 0.41 and 0.56 mmHg/heart rate, respectively. Predicted Fick-MPLV differed by 4.7%, Thermodilution-MPLV by 30% and Doppler-MPLV by 12%, when compared to clinical reports. CONCLUSIONS: Preliminary results from one subject show results that are in the range of literature values. The method needs to be validated by further testing, including independent measurements of intraventricular pressure. Since flow depends on the pressure loads, measuring more accurate intraventricular pressures helps to understand the cardiac flow dynamics for better clinical diagnosis. Furthermore, the method is non-invasive, safe, cheap and more practical. As clinical Fick-measured values have been known to be more accurate, our Fick-based prediction could be the most applicable.

A finite element investigation on plaque vulnerability in realistic healthy and atherosclerotic human coronary arteries.

Atherosclerosis is the most common arterial disease. It has been shown that stresses that are induced during blood circulation can cause plaque rupture and, in turn, lead to thrombosis and stroke. In this study, finite element method is used to predict plaque vulnerability based on peak plaque stress using human samples. A total of 23 healthy and atherosclerotic human coronary arteries of 14 healthy and 9 atherosclerotic patients are excised within 5 h postmortem. The samples are mounted on an uniaxial tensile test machine, and the obtained mechanical properties are used in two-dimensional and three-dimensional finite element models. The results including the Neo-Hookean hyperelastic coefficients of the samples as well as peak plaque stresses are analyzed. The results indicate that the atherosclerotic human coronary arteries have significantly (p < 0.05) higher stiffness compared with the healthy ones. The hypocellular plaque also has the highest stress values and, as a result, is most likely (vulnerable) to rupture, while the calcified type has the lowest stress values and, consequently, is expected to remain stable. The results could be used in the plaque vulnerability anticipation and have clinical implications in interventions and surgeries, including balloon angioplasty, bypass, and stenting.

An investigation of cerebral bridging veins rupture due to head trauma.

Subdural hematoma (SDH) is common abnormality that is caused by the rupture of cerebral bridge veins (BVs). It occurs in more than 30% of severe head injuries. The purpose of this research was to develop a numerical model to examine the effects of brain atrophy and age on the rupture of bridging veins in subdural hematoma. Three types of models were developed to simulate subdural hematoma, namely global solid, global FSI, and local solid models. In the next step, a head impact with the head injury criterion (HIC) value of 744 was applied as a loading condition to global models. For the global solid models, we measured the relative displacement between the skull and brain. We extracted the pressure distribution from the global FSI models. The data were used as boundary conditions on the local models to evaluate the damage to the cerebral bridge veins precisely The results showed that the relative displacement was greater in the atrophied model compared to the healthy one (2.64 and 2.20 mm, respectively). In addition, the pressure value was higher in atrophied models. In the healthy local model, the maximum strain on BVs was around 1.38, while in the atrophied model, it was 2.77. The head impact, which had a HIC value of 744, did not cause serious injury to a human with a healthy brain, but it caused severe damage to an atrophied brain. The degeneration of the brain and intracranial space changes are two important factors for the movement of the brain and its vulnerability to impact.

Biomechanical analysis of tracheal stent during cough reflex.

Tracheal stenting is a common method which is widely used to cure different tracheal disorders including airways stenosis, chronic coughs, and accidents. In this study, we aimed to analyze the reaction of the trachea wall to exhale in three phases of light, moderate, and vigorous activities at air flows of 15 L/min (light), 26 L/min (medium), and 30 L/min (vigorous). Fluid structure interaction (FSI) was used for the numerical analysis using computed tomography (CT) images. The flow was assumed incompressible and turbulent. The stent is silicone with a Young's modulus equal to 1 MPa, Poisson's ratio 0.28, and density of 2330 kg/m3. The stent length was 60 mm and fix support boundary condition was applied for all inputs and outputs. Numerical simulation was performed using ANSYS software. The induced stresses, strains, wall deformation, flow pressure, and the flow velocity were obtained. The results showed that the stent prevented the local deformation of the wall of trachea and it reduced the induced strain in the position. But the stenting could lead to stress concentration. Finally, the stent prevented the damage to the trachea muscles during coughs in row.

A color spectrographic phonocardiography (CSP) applied to the detection and characterization of heart murmurs: preliminary results.

BACKGROUND: Although cardiac auscultation remains important to detect abnormal sounds and murmurs indicative of cardiac pathology, the application of electronic methods remains seldom used in everyday clinical practice. In this report we provide preliminary data showing how the phonocardiogram can be analyzed using color spectrographic techniques and discuss how such information may be of future value for noninvasive cardiac monitoring. METHODS: We digitally recorded the phonocardiogram using a high-speed USB interface and the program Gold Wave http://www.goldwave.com in 55 infants and adults with cardiac structural disease as well as from normal individuals and individuals with innocent murmurs. Color spectrographic analysis of the signal was performed using Spectrogram (Version 16) as a well as custom MATLAB code. RESULTS: Our preliminary data is presented as a series of seven cases. CONCLUSIONS: We expect the application of spectrographic techniques to phonocardiography to grow substantially as ongoing research demonstrates its utility in various clinical settings. Our evaluation of a simple, low-cost phonocardiographic recording and analysis system to assist in determining the characteristic features of heart murmurs shows promise in helping distinguish innocent systolic murmurs from pathological murmurs in children and is expected to useful in other clinical settings as well.

Digital Subtraction Phonocardiography (DSP) applied to the detection and characterization of heart murmurs.

BACKGROUND: During the cardiac cycle, the heart normally produces repeatable physiological sounds. However, under pathologic conditions, such as with heart valve stenosis or a ventricular septal defect, blood flow turbulence leads to the production of additional sounds, called murmurs. Murmurs are random in nature, while the underlying heart sounds are not (being deterministic). INNOVATION: We show that a new analytical technique, which we call Digital Subtraction Phonocardiography (DSP), can be used to separate the random murmur component of the phonocardiogram from the underlying deterministic heart sounds. METHODS: We digitally recorded the phonocardiogram from the anterior chest wall in 60 infants and adults using a high-speed USB interface and the program Gold Wave http://www.goldwave.com. The recordings included individuals with cardiac structural disease as well as recordings from normal individuals and from individuals with innocent heart murmurs. Digital Subtraction Analysis of the signal was performed using a custom computer program called Murmurgram. In essence, this program subtracts the recorded sound from two adjacent cardiac cycles to produce a difference signal, herein called a "murmurgram". Other software used included Spectrogram (Version 16), GoldWave (Version 5.55) as well as custom MATLAB code. RESULTS: Our preliminary data is presented as a series of eight cases. These cases show how advanced signal processing techniques can be used to separate heart sounds from murmurs. Note that these results are preliminary in that normal ranges for obtained test results have not yet been established. CONCLUSIONS: Cardiac murmurs can be separated from underlying deterministic heart sounds using DSP. DSP has the potential to become a reliable and economical new diagnostic approach to screening for structural heart disease. However, DSP must be further evaluated in a large series of patients with well-characterized pathology to determine its clinical potential.

A comparative finite element simulation of locking compression plate materials for tibial fracture treatment.

The locking compression plate (LCP) system has several advantages in fracture fixation combining angular stability with the use of locking screws with traditional fixation techniques. However, the system is complex and requiring careful attention to biomechanical principles and good surgical technique. Due to the set of complicate stresses and strains in the LCP system after implantation, the material, which is being used here, is deemed important. However, so far the materials have been limited to the stainless steel (SS) or titanium (Ti). This study was therefore aimed at investigate the biomechanical performance of the internal tibial locked plates at different material properties, including SS, Ti, carbon/polyether ether ketone (PEEK) composite, in treating medial tibial fracture using patient-specific finite element (FE) model of the human tibia. The carbon/PEEK composite materials were used at three different fiber plies configurations. Simulated loading was applied at 60:40 ratios on the medial:lateral aspect. The model was fixed distally in all degrees of freedom. The results revealed the highest stress (307.10 MPa) and the lowest strain (0.14%) at Ti LCP system. The carbon/PEEK LCP system at configuration I and III showed low stress (∼60 MPa) and high strain (0.70%), which are suitable points for designing of an internal LCP system. On the other hand, the highest value of stress in callus region was 4.78 MPa (Carbon PEEK/Configuration I) and the strain variations of callus region were between 1.46% and 3.82% among all materials. These results implied the advantage of carbon/PEEK composite materials in LCP system as they can tolerate higher strains at lower stresses.

A finite element study of fatigue load effects on total hip joint prosthesis.

The main goal of this study was to perform a fatigue analysis and compare the results for different materials. A 72 years old patient was chosen and his hip radiographic/CT scan images were used to construct the geometry. We used four different material including Titanium, Titanium alloy, Cobalt-Chrome, and Stainless steel. The material characteristics of these prostheses were extracted from the literature. All models were exported to ANSYS software for mathematical analysis and the Von-Mises criteria, deformations, and the fatigue life were calculated for each material. Our findings showed that titanium prosthesis tolerated the lowest stress (i.e., 591 MPa for static, and 687 MPa for fatigue loading) and highest safety factor (i.e., 1.54). We found out that Titanium material could be used as the most appropriate one for the hip prosthesis due to lower stress concentration and longer life compared to other materials.

Biomechanical role of posterior cruciate ligament in total knee arthroplasty: A finite element analysis.

BACKGROUND AND OBJECTIVE: The knee joint is a complex structure which is vulnerable to injury due to various types of loadings as a consequence of walking, running, stair climbing, etc. Total knee arthroplasty (TKA) is a widely used and successful orthopedic procedure which during that the posterior cruciate ligament (PCL) can either be retained or substituted. Different surgical techniques suggest retention or sacrifice of the PCL in TKA for the treatment of osteoarthritis which may alter the post-op outcomes. The objective of this study was to evaluate the biomechanical role of PCL after TKA surgery using finite element (FE) modeling. METHODS: A three-dimensional (3D) FE model of the prosthetic knee was developed and its validity was compared to available studies in literature. Further, the effect of the retention or removing of the PCL as well as its degradation (i.e. variation in mechanical properties) and angle on knee biomechanics were evaluated during a weight-bearing squatting movement. RESULTS: The validity of the intact model were confirmed. The results revealed higher stresses in the PCL and tibial insert at higher femoral flexion angles. In addition, the effect of variations in the stiffness of the PCL was found to be negligible at lower while considerable at higher femoral flexion angles. The variations in the elevation angle of the PCL from 89° to 83° at the critical femoral angles of 60° and 120° showed the highest von Mises stresses in the tibial insert. CONCLUSIONS: The results have implications not only for understanding the stresses in the prosthetic knee model under squat movement but also for providing comprehensive information about the effects of variations in the PCL stiffness and balancing on the induced stresses of the PCL and tibial insert.

In silico investigation of sneezing in a full real human upper airway using computational fluid dynamics method.

BACKGROUND AND OBJECTIVE: Sneezing is one of the most critical conditions that can occur in the human upper airway. As some reports confirm the injury to the human upper respiratory airway while sneezing. Therefore, the accurate study of the distribution of pressure and velocity in this case is of great importance. METHODS: In the present study, using a real human upper airway model, the pressure and velocity of the airflow generated in the tract during the sneezing have been investigated. Also, considering the results from a spirometer device as a boundary condition in the simulation process, the calculations have become reliable. RESULTS: According to the results, during sneezing, taking into account that the average outlet flow rate from the mouth is 4.79 L/s, the velocity of outlet airflow from the mouth and nose reaches 5.3 and 8.4 m/s, respectively. These values were 11.5 and 19, respectively, when the desired maximum flow rate was 10.58 L/s. It can be concluded that the increasing of trachea flow rate, leads to higher percentage of the outlet flow rate from the nose . The highest average pressure and velocity have been occurred in the trachea. Among other salient results of this report, increased average static pressure of larynx to approximately 10 kPa can be pointed which indicates that this area is critical so that the thyroid cartilage defect is likely to occur. It is also noteworthy that the increase of speed at nasopharynx is up to 125 m/s so that the cross-section changing in this area leads the fluid acts as a jet flow. Due to the specific geometry of the nasal cavity, some streams similar to poor shocks are formed, these shocks get stronger by increasing of the flow rate. The thyroid cartilage and nasal cavity are exposed to maximum static pressure extremums, respectively. CONCLUSIONS: We introduced a model simulating a normal sneezing for two cases using a healthy 30-year-old male person. We believe that the model should be applied for different persons and an atlas of data could be obtained from different cases. This may help the medical system to have more data about the sneezing process.

Fluid-structure interaction assessment of blood flow hemodynamics and leaflet stress during mitral regurgitation.

The aim of this study is to simulate the Mitral Regurgitation (MR) disease progression from mild to severe intensity. A Fluid Structure Interaction (FSI) model was developed to extract the hemodynamic parameters of blood flow in mitral regurgitation (MR) during systole. A two-dimensional (2D) geometry of the mitral valve was built based on the data resulting from Magnetic Resonance Imaging (MRI) dimensional measurements. The leaflets were assumed to be elastic. Using COMSOL software, the hemodynamic parameters of blood flow including velocity, pressure, and Von Mises stress contours were obtained by moving arbitrary Lagrange-Euler mesh. The results were obtained for normal and MR cases. They showed the effects of the abnormal distance between the leaflets on the amount of returned flow. Furthermore, the deformation of the leaflets was measured during systole. The results were found to be consistent with the relevant literature.