Computer simulation-based nanothermal field and tissue damage analysis for cardiac tumor ablation

Radiofrequency ablation is a nominally invasive technique to eradicate cancerous or non-cancerous cells by heating. However, it is still hampered to acquire a successful cell destruction process due to inappropriate RF intensities that will not entirely obliterate tumorous tissues, causing in treatment failure. In this study, we are acquainted with a nanoassisted RF ablation procedure of cardiac tumor to provide better outcomes for long-term survival rate without any recurrences. A three-dimensional thermo-electric energy model is employed to investigate nanothermal field and ablation efficiency into the left atrium tumor. The cell death model is adopted to quantify the degree of tissue injury while injecting the Fe3O4 nanoparticles concentrations up to 20% into the target tissue. The results reveal that when nanothermal field extents as a function of tissue depth (10 mm) from the electrode tip, the increasing thermal rates were approximately 0.54362%, 3.17039%, and 7.27397% for the particle concentration levels of 7%, 10%, and 15% compared with no-particle case. In the 7% Fe3O4 nanoparticles, 100% fractional damage index is achieved after ablation time of 18 s whereas tissue annihilation approach proceeds longer to complete for no-particle case. The outcomes indicate that injecting nanoparticles may lessen ablation time in surgeries and prevent damage to adjacent healthy tissue.

Surrogate metamodels from digital image correlation for testing high-performance composite vessels

In this work, a workflow has been developed for the generation of surrogate metamodels to predict and evaluate failure with a confidence above 95 % in initial service conditions of high-performance cylindrical vessels manufactured in composites by Roll Wrapping technology. Currently, there is no specific testing standardization for this type of vessel and to fill this gap probabilistic numerical models were developed, performed by the Finite Element Method, fed with the material characteristics obtained experimentally by 2D digital image correlation from flat specimens. From the initial numerical model, a surrogate metamodel was generated by stochastic approximations. Once the metamodels were obtained by robust engineering, an experimental ring-ring tensile test was developed under service conditions and deformations were measured by high-precision 3D digital image correlation. Parametric and robust tests showed that the results of the metamodel did not show statistically significant differences, with errors in the rupture part of less than 2 % with respect to the results obtained in the test, being proposed as a basis for new test procedures.

Evaluating the effect of functionally graded materials on bone remodeling around dental implants

OBJECTIVES

This study evaluates the potential for osseointegration and remodeling of customized dental implants made from Titanium-Hydroxyapatite Functionally Graded Material (Ti-HAP FGM) with optimized geometry, using the finite element method (FEM).

METHODS

The study utilized CT scan images to model and assemble various geometrical designs of dental implants in a mandibular slice. The mechanical properties of Ti-HAP FGMs were computed by varying volume fractions (VF) of hydroxyapatite (0-20%), and a bone remodeling algorithm was used to evaluate the biomechanical characteristics of the ultimate bone configuration in the peri-implant tissue.

RESULTS

The findings of the FEA reveal that osseointegration improves with changes in the density and mechanical properties of the bone surrounding Ti-HAP implants, which are influenced by the varying VF of hydroxyapatite in the FGM.

SIGNIFICANCE

Increasing the hydroxyapatite fraction improves osseointegration, and appropriate length and diameter selection of Ti-HAP dental implants contribute to their stability and longevity.

Three-dimensional theoretical model for effectively describing the effect of craniomaxillofacial structural factors on loading situation in the temporomandibular joint

BACKGROUND

Temporomandibular joint (TMJ) overloading is considered a primary cause of temporomandibular joint disorders (TMD). Accordingly, craniomaxillofacial structural parameters affect the loading situation in the TMJ. However, no effective method exists for quantitatively measuring the loading variation in human TMJs. Clinical statistics, which draws from general rules from large amounts of clinical data, cannot entry for exploring the underlying biomechanical mechanism in craniomaxillofacial system. The finite element method (FEM) is an effective tool for analyze the stress and load on TMJs for several cases in a short period of time; however, it is difficult to generalize general patterns through calculations between different cases due to the different geometric characteristics and occlusal contacts between each case.

METHODS

(1) This study included 88 subjects with 176 unilateral data to measure angle (α) of the distance to the plane of occlusion. The bone destruction score was evaluated for clinical statistics. To rule out effects of the potential factors and ensure the generality of the study, one participant with no obvious bone destruction was selected as the standard case for establishing the three-dimensional (3D) theoretical model and FEM. (2) Three groups of forces, including biting, muscles and joint reaction forces on mandible, were adopted to establish a 3D theoretical model. (3) By modifying the sagittal α and coronal three types of deviation angle (φ) of the original model, nine candidate models were obtained for the FEM studies.

RESULTS

(1) The static equilibrium equations, were used to establish a 3D theoretical model for describing the loading of the TMJ. The theoretical model was validated by monotonously modifying the structural parameter in comparison to two-dimensional theoretical models reported previously; (2) The force on the TMJ gradually decreased with α, and this trend was validated by both clinic statistics and FEM results; (3) The effects of the three types of deviation angle were different. The results of the case where only rotating biting forces were considered was consistent with clinical statistics, indicating that the side with lower α experiences higher TMJ load. (4) Changing the unilateral proportionality coefficients of biting and muscle force produced opposite effects, wherein the effects of the muscle force were stronger than those of the biting forces.

CONCLUSIONS

A negative correlation was observed between the joint load and α. Among the three types of asymmetric deformities, occlusal deviations were the primary factors leading to TMD. Unilateral occlusion can result in a greater load on the ipsilateral joint and should be avoided when using the side corresponding to the TMD. This study provides a theoretical basis for the biomechanical mechanism of TMD and also enables the targeted mitigation and treatment of TMD through structural modification.

Enhanced mixing efficiency and reduced droplet size with novel droplet generators

Nowadays, droplet microfluidics has become widely utilized for high-throughput assays. Efficient mixing is crucial for initiating biochemical reactions in many applications. Rapid mixing during droplet formation eliminates the need for incorporating micromixers, which can complicate the chip design. Furthermore, immediate mixing of substances upon contact can significantly improve the consistency of chemical reactions and resulting products. This study introduces three innovative designs for droplet generators that achieve efficient mixing and produce small droplets. The T-cross and cross-T geometries combine cross and T junction mixing mechanisms, resulting in improved mixing efficiency. Numerical simulations were conducted to compare these novel geometries with traditional T and cross junctions in terms of mixing index, droplet diameter, and eccentricity. The cross-T geometry exhibited the highest mixing index and produced the smallest droplets. For the flow rate ratio of 0.5, this geometry offered a 10% increase in the mixing index and a decrease in the droplet diameter by 10% compared to the T junction. While the T junction has the best mixing efficiency among traditional droplet generators, it produces larger droplets, which can increase the risk of contamination due to contact with the microchannel walls. Therefore, the cross-T geometry is highly desirable in most applications due to its production of considerably smaller droplets. The asymmetric cross junction offered a 8% increase in mixing index and around 2% decrease in droplet diameter compared to the conventional cross junction in flow rate ratio of 0.5. All novel geometries demonstrated comparable mixing efficiency to the T junction. The cross junction exhibited the lowest mixing efficiency and produced larger droplets compared to the cross-T geometry (around 1%). Thus, the novel geometries, particularly the cross-T geometry, are a favorable choice for applications where both high mixing efficiency and small droplet sizes are important.

Investigating the effects of solution viscosity on the stability and success rate of SECCM imaging

Due to the capability of simultaneously detecting the morphology and electrochemical information of samples and limiting the electrochemical reaction to a range approximately the size of the inner diameter of the pipette tip opening, scanning electrochemical cell microscopy (SECCM) enables higher precision local electrochemical measurement and surface material delivery and has been demonstrating unique advantages and broad application prospects. However, the meniscus droplet at the pipette tip of SECCM is equivalent to the opening radius of the pipette tip, which is usually tens of nanometers to hundreds of nanometers. The tiny meniscus droplet makes it susceptible to evaporation and crystallization, which increases the likelihood of the pipette colliding with the sample during the scanning process, resulting in the failure of scanning. In this paper, the influence of solution viscosity on the shape variation of the droplet at the tip during the movement of the pipette of SECCM was studied by finite element analysis. It is proved that the increase of solution viscosity is helpful in reducing the shape variation of the droplet at the tip during the movement of the pipette. Then scanning experiments were carried out using a flat Au substrate and Au substrates with rounded triangle and rounded rectangular convex structures as the samples. According to the experimental results, increasing solution viscosity improves scanning success rates and scanning quality and effectively lowers the MSE of the scanning results. The experimental results also show that SECCM can image at a higher speed when the solution's viscosity increases since the deformation of the droplet at the tip is less than with a typical solution.

Biomechanical effects of joint disc perforation on the temporomandibular joint: a 3D finite element study

BACKGROUND

Disc perforation (DP) is a severe type of Temporomandibular Disorder (TMD). DP may induce changes in the internal stresses of the temporomandibular joint (TMJ). Herein, this study attempts to investigate the biomechanical effects of different positions and sizes of DP on the TMJ using a biomechanical approach, to explore the mechanical pathogenesis of TMD.

METHODS

Eleven three-dimensional finite element (FE)models of the TMJ were constructed based on CBCT imaging files of a patient with DP on the left side. These models included the disc with anterior displacement and discs with different locations and sizes of perforations on the affected disc. FE methods were conducted on these models.

RESULTS

Anterior displacement of the disc leads to a significant increase in the maxim von Mises stress (MVMS) in both TMJs, with the affected side exhibiting a more pronounced effect. DP occurring at the posterior band and the junction between the disc and the bilaminar region has a greater impact on the MVMS of both TMJs compared to perforations at other locations. As the size of the perforation increases, both sides of the TMJs exhibit an increase in the magnitude of MVMS.

CONCLUSIONS

Unilateral disc anterior displacement results in an increased stress on both TMJs. Unilateral DP further affects the stress on both sides of the TMJs. TMD is a progressive condition, and timely intervention is necessary in the early stages to prevent the worsening of the condition.

Outcome measures for electric field modeling in tES and TMS: A systematic review and large-scale modeling study

BACKGROUND

Electric field (E-field) modeling is a potent tool to estimate the amount of transcranial magnetic and electrical stimulation (TMS and tES, respectively) that reaches the cortex and to address the variable behavioral effects observed in the field. However, outcome measures used to quantify E-fields vary considerably and a thorough comparison is missing.

OBJECTIVES

This two-part study aimed to examine the different outcome measures used to report on tES and TMS induced E-fields, including volume- and surface-level gray matter, region of interest (ROI), whole brain, geometrical, structural, and percentile-based approaches. The study aimed to guide future research in informed selection of appropriate outcome measures.

METHODS

Three electronic databases were searched for tES and/or TMS studies quantifying E-fields. The identified outcome measures were compared across volume- and surface-level E-field data in ten tES and TMS modalities targeting two common targets in 100 healthy individuals.

RESULTS

In the systematic review, we extracted 308 outcome measures from 202 studies that adopted either a gray matter volume-level (n = 197) or surface-level (n = 111) approach. Volume-level results focused on E-field magnitude, while surface-level data encompassed E-field magnitude (n = 64) and normal/tangential E-field components (n = 47). E-fields were extracted in ROIs, such as brain structures and shapes (spheres, hexahedra and cylinders), or the whole brain. Percentiles or mean values were mostly used to quantify E-fields. Our modeling study, which involved 1,000 E-field models and > 1,000,000 extracted E-field values, revealed that different outcome measures yielded distinct E-field values, analyzed different brain regions, and did not always exhibit strong correlations in the same within-subject E-field model.

CONCLUSIONS

Outcome measure selection significantly impacts the locations and intensities of extracted E-field data in both tES and TMS E-field models. The suitability of different outcome measures depends on the target region, TMS/tES modality, individual anatomy, the analyzed E-field component and the research question. To enhance the quality, rigor, and reproducibility in the E-field modeling domain, we suggest standard reporting practices across studies and provide four recommendations.

Image data and computational grids for computing brain shift and solving the electrocorticography forward problem

This article describes the dataset applied in the research reported in NeuroImage article "Patient-specific solution of the electrocorticography forward problem in deforming brain" [1] that is available for download from the Zenodo data repository (https://zenodo.org/record/7687631) [2]. Preoperative structural and diffusion-weighted magnetic resonance (MR) and postoperative computed tomography (CT) images of a 12-year-old female epilepsy patient under evaluation for surgical intervention were obtained retrospectively from Boston Children's Hospital. We used these images to conduct the analysis at The University of Western Australia's Intelligent Systems for Medicine Laboratory using SlicerCBM [3], our open-source software extension for the 3D Slicer medical imaging platform. As part of the analysis, we processed the images to extract the patient-specific brain geometry; created computational grids, including a tetrahedral grid for the meshless solution of the biomechanical model and a regular hexahedral grid for the finite element solution of the electrocorticography forward problem; predicted the postoperative MRI and DTI that correspond to the brain configuration deformed by the placement of subdural electrodes using biomechanics-based image warping; and solved the patient-specific electrocorticography forward problem to compute the electric potential distribution within the patient's head using the original preoperative and predicted postoperative image data. The well-established and open-source file formats used in this dataset, including Nearly Raw Raster Data (NRRD) files for images, STL files for surface geometry, and Visualization Toolkit (VTK) files for computational grids, allow other research groups to easily reuse the data presented herein to solve the electrocorticography forward problem accounting for the brain shift caused by implantation of subdural grid electrodes.

Biomechanical analysis of a single-level Customized cage screw fixation for Anterior Cervical Discectomy and Fusion in the cervical spine: an in-silico study

The usual surgical choice of customized cage fixation is anterior cervical discectomy and fusion (ACDF) for cervical spondylosis with disc herniation. The safe and successful cage fixation for ACDF surgery benefits patients with cervical disc degenerative disease by easing their discomfort and regaining function. The cage prevents mobility between the vertebrae by using cage fixation to anchor the neighbouring vertebrae. The goal of the current study is to develop a customized cage-screw implant for single-level cage fixation at C4-C5 level of the cervical spine (C2-C7). The Finite Element Analysis (FEA) is performed for the intact and implanted cervical spine and analysed the flexibility, stress of the implant and implant adjacent bone during three physiological loading conditions being analysed. Lower surface of the C7 vertebrae is fixed and 50 N compressive force with 1 Nm moment are applied on the C2 vertebrae for simulated lateral bending, axial rotation and flexion-extension. The flexibility is decreased at single level of fixation (C4-C5 level) by 64% to 86%, as compared to natural cervical spine. The flexibility is increased 3% to 17% at the nearest levels of fixation. The maximum Von Mises stress in PEEK cage varies from 24 to 59 MPa and for Ti-6Al-4V screw the stress varies from 84MPa to 121MPa which are far below the yield stress of PEEK (95MPa) and Ti-6Al-4V (750MPa).

Choosing the optimal mandible position for inferior alveolar nerve block (IANB) using finite element analysis

BACKGROUND

One of the most popular methods of local anesthesia in dentistry, inferior alveolar nerve block (IANB) involves the blockade of the inferior alveolar nerve (IAN) and lingual nerve (LN) in the pterygomandibular space. Despite the large number of works describing the contents of this space, the spatial displacements of the anatomical structures of this area at different positions of the mandible have not been sufficiently studied. The aim of our study was to study the spatial movements of the IAN and inferior alveolar artery (IAA) at various positions of the mandible using computer simulation and finite element analysis to find the safest way to conduct IANB.

MATERIALS AND METHODS

Reverse engineering was used to create a model of the cranial base and the mandible based on computed tomography (CT) data obtained from patient N (male, 24 years old), the arteries of the head and neck were designed from the data of multiphase angiography of patient M (female, 61 years old). Masticatory muscles, sphenomandibular ligament, temporomandibular joint and mandibular nerve were modeled in the SolidWorks software package based on an open database of anatomical structures. The finite element grid was generated in the Solidworks software. In the first series of experiments, the displacement of the mandible was modeled along the vertical axis down by 48 mm (maximum opening of the mouth), in the second series, the jaw was displaced vertically by 48 mm with a simultaneous transversal movement of 10 mm, in the third series, the jaw was displaced along the vertical axis down by 34 mm and transversally by 22 mm.

RESULTS

The largest distance between IAN and IAA was noted in the third series of experiments. The distance between the nerve and the vessel was minimal in the first series, with an open mouth without lateral displacements.

CONCLUSION

The generated computer model opens new possibilities for studying the dynamic anatomy of the pterygomandibular space. The results of this study can be used for further experimental and clinical trials to find the safest approach to the implementation of IANB, as well as applied in the practice of the educational process.

The impact of the parameters of the constitutive model on the distribution of strain in the femoral head

The rapid spread of the finite element method has caused that it has become, among other methods, the standard tool for pre-clinical estimates of bone properties. This paper presents an application of this method for the calculation and prediction of strain and stress fields in the femoral head. The aim of the work is to study the influence of the considered anisotropy and heterogeneity of the modeled bone on the mechanical fields during a typical gait cycle. Three material models were tested with different properties of porous bone carried out in literature: a homogeneous isotropic model, a heterogeneous isotropic model, and a heterogeneous anisotropic model. In three cases studied, the elastic properties of the bone were determined basing on the Zysset-Curnier approach. The tensor of elastic constants defining the local properties of porous bone is correlated with a local porosity and a second order fabric tensor describing the bone microstructure. In the calculations, a model of the femoral head generated from high-resolution tomographic scans was used. Experimental data were drawn from publicly available database "Osteoporotic Virtual Physiological Human Project." To realistically reflect the load on the femoral head, main muscles were considered, and their contraction forces were determined based on inverse kinematics. For this purpose, the results from OpenSim packet were used. The simulations demonstrated that differences between the results predicted by these material models are significant. Only the anisotropic model allowed for the plausible distribution of stresses along the main trabecular groups. The outcomes also showed that the precise evaluation of the mechanical fields is critical in the context of bone tissue remodeling under mechanical stimulations.

Quantitative Comparison of Analytical Solution and Finite Element Method for Investigation of Near-infrared Light Propagation in Brain Tissue Model

Introduction

Functional near-infrared spectroscopy (fNIRS) is an imaging method in which a light source and detector are installed on the head; consequently, the re-emission of light from human skin contains information about cerebral hemodynamic alteration. The spatial probability distribution profile of photons penetrating tissue at a source spot, scattering into the tissue, and being released at an appropriate detector position, represents the spatial sensitivity.

Methods

Modeling light propagation in a human head is essential for quantitative near-infrared spectroscopy and optical imaging. The specific form of the distribution of light is obtained using the theory of perturbation. An analytical solution of the perturbative diffusion equation (DE) and finite element method (FEM) in a Slab media (similar to the human head) makes it possible to study light propagation due to absorption and scattering of brain tissue.

Results

The simulation result indicates that sensitivity is slowly decreasing in the deep area, and the sensitivity below the source and detector is the highest. The depth sensitivity and computation time of both analytical and FEM methods are compared. The simulation time of the analytical approach is four times larger than the FEM.

Conclusion

In this paper, an analytical solution and the performance of FEM methods when applied to the diffusion equation for heterogeneous media with a single spherical defect are compared. The depth sensitivity along with the computation time of simulation has been investigated for both methods. For simple and Slab modes of the human brain, the analytical solution is the right candidate. Whenever the brain model is sophisticated, it is possible to use FEM methods, but it costs a higher computation time.

Highlights

Analytical and finite element method (FEM) depth sensitivity are almost the same.FEM requires more computation time, but can handle complicated head models.The analytical method is proposed for the first step and simple head models.

Plain Language Summary

The functional near-infrared spectroscopy (fNIRS) is a type of neuromonitoring that uses near-infrared light to measure brain activity indirectly and is similar to electroencephalography (EEG). A single-channel fNIRS system contains a near-infrared light source, which emits near-infrared light (NIR), and a detector is placed near the source. A light intensity change received by detectors indicates brain activity when NIR light penetrates into the gray matter. It is necessary to have a prior understanding of light penetration depth in order to measure brain activity more accurately. fNIRS can be better understood, optimized, and investigated through modeling light propagation in brain tissue. In order to study light in tissues, analytical and numerical methods can be used. In this work, we compared these two approaches quantitatively in a simple slab medium. We concluded that the numerical method takes too much time to calculate the results, but it can be applied to complicated head models. The results of these studies provide researchers with new insights into the modeling and simulation of fNIRS and diffuse optical tomography.

A Fluid-Structure Interaction Analysis of Blood Clot Motion in a Branch of Pulmonary Arteries

INTRODUCTION

Pulmonary embolism (PE) is one of the most prevalent diseases amid hospitalized patients taking many people's lives annually. This phenomenon, however, has not been investigated via numerical simulations.

METHODS

In this study, an image-based model of pulmonary arteries has been constructed from a 44-year-old man's computed tomography images. The fluid-structure interaction method was used to simulate the motion of the blood clot. In this regard, Navier-Stokes equations, as the governing equations, have been solved in an arbitrary Lagrangian-Eulerian (ALE) formulation.

RESULTS

According to our results, the velocity of visco-hyperelastic model of the emboli was relatively higher than the emboli with hyperelastic model, despite their similar behavioral pattern. The stresses on the clot were also investigated and showed that the blood clot continuously sustained stresses greater than 165 Pa over an about 0.01 s period, which can cause platelets to leak and make the clot grow or tear apart.

CONCLUSIONS

It could be concluded that in silico analysis of the cardiovascular diseases initiated from clot motion in blood flow is a valuable tool for a better understanding of these phenomena.

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.

Can individually targeted and optimized multi-channel tDCS outperform standard bipolar tDCS in stimulating the primary somatosensory cortex?

BACKGROUND

Transcranial direct current stimulation (tDCS) has emerged as a non-invasive neuro-modulation technique. Most studies show that anodal tDCS increases cortical excitability, however, with variable outcomes. Previously, we have shown in computer simulations that our multi-channel tDCS (mc-tDCS) approach, the distributed constrained maximum intensity (D-CMI) method can potentially lead to better controlled tDCS results due to the improved directionality of the injected current at the target side for individually optimized D-CMI montages.

OBJECTIVE

In this study, we test the application of the D-CMI approach in an experimental study to stimulate the somatosensory P20/N20 target source in Brodmann area 3b and compare it with standard bipolar tDCS and sham conditions.

METHODS

We applied anodal D-CMI, the standard bipolar and D-CMI based Sham tDCS for 10 min to target the 20 ms post-stimulus somatosensory P20/N20 target brain source in Brodmann area 3b reconstructed using combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis in realistic head models with calibrated skull conductivity in a group-study with 13 subjects. Finger-stimulated somatosensory evoked fields (SEF) were recorded and the component at 20 ms post-stimulus (M20) was analyzed before and after the application of the three tDCS conditions in order to read out the stimulation effect on Brodmann area 3b.

RESULTS

Analysis of the finger stimulated SEF M20 peak before (baseline) and after tDCS shows a significant increase in source amplitude in Brodmann area 3b for D-CMI (6-16 min after tDCS), while no significant effects are found for standard bipolar (6-16 min after tDCS) and sham (6-16 min after tDCS) stimulation conditions. For the later time courses (16-26 and 27-37 min post-stimulation), we found a significant decrease in M20 peak source amplitude for standard bipolar and sham tDCS, while there was no effect for D-CMI.

CONCLUSION

Our results indicate that targeted and optimized, and thereby highly individualized, mc-tDCS can outperform standard bipolar stimulation and lead to better control over stimulation outcomes with, however, a considerable amount of additional work compared to standard bipolar tDCS.

Condyle bone destruction: the association between temporomandibular joint vibration and finite element analysis

OBJECTIVE

The objective of the study is to investigate the stress distributions of condyle and articular disc at different mandibular plane angles and the association between the temporomandibular joint (TMJ) vibration and anatomical relationship.

SETTING AND SAMPLE POPULATION

Data from 195 untreated patients were analyzed. Patients were, respectively, divided into 3 groups, based on cone-beam computed tomography (CBCT): Group 0 presenting normal condyles, Group 1 presenting mild bone change, and Group 2 presenting severe bone change; based on magnetic resonance imaging (MRI): normal disc position (N), disc displacement with reduction (DDR) and disc displacement without reduction (DDNR); based on the mandibular plane angle: low, normal and high.

METHODS

(1) According to peak frequency, average amplitude, and occurrence phase, association was assessed between TMJ condition and TMJ vibration; (2) A three-dimensional finite element model of masticatory system was established and the relationship between mandibular plane angle and condyle bone destruction was described.

RESULTS

The average amplitude of TMJ vibration negatively correlated with pathological condition of the articular disc (p < 0.01). The angle of SN-MP was strongly relevant to bone destruction of condyle (p < 0.01), and the biomechanical analysis showed that with the increase of SN-MP angle, the area of stress concentration in the posterior slope of condyle rose.

CONCLUSIONS

The average amplitude of TMJ vibration decreased with the pathological position state of articular disc, and condyle bone destruction was positively associated with SN-MP angle.

Probabilistic maps for deep brain stimulation - Impact of methodological differences

BACKGROUND

Group analysis of patients with deep brain stimulation (DBS) has the potential to help understand and optimize the treatment of patients with movement disorders. Probabilistic stimulation maps (PSM) are commonly used to analyze the correlation between tissue stimulation and symptomatic effect but are applied with different methodological variations.

OBJECTIVE

To compute a group-specific MRI template and PSMs for investigating the impact of PSM model parameters.

METHODS

Improvement and occurrence of dizziness in 68 essential tremor patients implanted in caudal zona incerta were analyzed. The input data includes the best parameters for each electrode contact (screening), and the clinically used settings. Patient-specific electric field simulations (n = 488) were computed for all DBS settings. The electric fields were transformed to a group-specific MRI template for analysis and visualization. The different comparisons were based on PSMs representing occurrence (N-map), mean improvement (M-map), weighted mean improvement (wM-map), and voxel-wise t-statistics (p-map). These maps were used to investigate the impact from input data (clinical/screening settings), clustering methods, sampling resolution, and weighting function.

RESULTS

Screening or clinical settings showed the largest impacts on the PSMs. The average differences of wM-maps were 12.4 and 18.2% points for the left and right sides respectively. Extracting clusters based on wM-map or p-map showed notable variation in volumes, while positioning was similar. The impact on the PSMs was small from weighting functions, except for a clear shift in the positioning of the wM-map clusters.

CONCLUSION

The distribution of the input data and the clustering method are most important to consider when creating PSMs for studying the relationship between anatomy and DBS outcome.

Euler-Lagrange numerical simulation of improved magnetic drug delivery in a three-dimensional CT-based carotid artery bifurcation

BACKGROUND AND OBJECTIVE

Magnetic drug targeting (MDT) is a promising method to improve the therapy efficiency for cardiovascular diseases (CVDs) and cancers. In MDT, therapeutic agents are bonded to superparamagnetic iron oxide nanoparticle (SPION) cores and then are guided toward the damaged tissue through a magnetic field. Fundamentally, it's vital to steer the SPIONs to the desired location to increase the capture efficiency at the target lesion. Hence, the present study aims to enhance the drug delivery to the desired branch in a carotid bifurcation. Besides, it is tried to decrement the particles' entry to the unwanted outlet by using four different magnet configurations (with a maximum magnetic flux density of 0.4 T) implanted adjacent to the artery wall. Also, the effect of particles' diameter -ranging from 200 nm to 2 µm- on the drug delivery performance is studied in the four cases.

METHODS

The Eulerian-Lagrangian approach with one-way coupling is employed for numerical simulation of the problem using the finite element method (FEM). The dominant forces acting on particles are drag and magnetophoretic. A computed tomography (CT) model of the carotid bifurcation is adopted to have a 3D realistic geometry. The blood flow is considered to be laminar, incompressible, pulsatile, and non-Newtonian. Boundary conditions are applied using the three-element Windkessel equation.

RESULTS

Results are presented in terms of velocity, pressure, magnetic field flux density, wall shear stress, and streamlines. Also, the number of particles in each direction is presented for the four studied cases. The results show that using proper magnets configurations makes it possible to guide more particles to the desired branch (up to 4%) while preventing particles from entering the unwanted branch (up to 13%). By defining connectivity between oscillatory shear index (OSI) value and magnetic drug delivery efficacy, it becomes clear that places with lower OSI values are more proper to place the magnets than areas with higher OSI values.

CONCLUSIONS

It was observed that increasing the diameter of particles does not necessarily result in a higher drug delivery efficiency. The configuration of the magnets and the size of particles are the main affecting parameters that should be chosen precisely to meet the most efficient drug delivery performance. Magnetic drug targeting (MDT) is a promising method to improve the therapy efficiency for cardiovascular diseases (CVDs) and cancers. Fundamentally, it's vital to steer the superparamagnetic iron oxide nanoparticles (SPIONs) to the target lesion location to increase the capture efficiency. Hence, the present study aims to enhance the drug delivery to the desired branch in a 3D carotid bifurcation. Besides, it is tried to decrement the particles' entry to the unwanted outlet by using four different magnet configurations implanted adjacent to the artery wall. The Eulerian-Lagrangian approach with one-way coupling is employed for numerical simulation of the problem using the finite element method (FEM). The dominant forces acting on particles are drag and magnetophoretic. The blood flow is laminar, incompressible, pulsatile, and non-Newtonian. The results show that it is possible to guide more particles to the desired branch (up to 4%) while preventing particles from entering the unwanted branch (up to 13%). By defining connectivity between oscillatory shear index (OSI) value and magnetic drug delivery efficacy, it becomes clear that places with lower OSI values are more proper to place the magnets than areas with higher OSI values.

Multipair transcranial temporal interference stimulation for improved focalized stimulation of deep brain regions: A simulation study

Temporal interference stimulation (TIS) has been proved to be effective in stimulating deep brain regions while avoiding the stimulation of neocortical regions in animal experiments. In the traditional TIS, two alternating currents are injected with different frequencies via two electrode pairs attached to the scalp. In the human brain, however, it is difficult to achieve a focal stimulation of deep brain structures due to the high complexity of human brain structures. In this study, we hypothesized that the use of multiple electrode pairs may contribute to the more focalized delivery of temporal interference (TI) currents to the target site in the deep area of the brain. Based on this hypothesis, we proposed a novel multipair TIS method that employs more than two electrode pairs for improved focalized stimulation of the deep brain region (in this study, the head of the right hippocampus). Three realistic finite element models were used to validate the feasibility of the proposed multipair TIS. Additional electrode pairs were sequentially added to the conventional two-electrode pairs with the aim of maximizing the delivery of TI currents to the target while minimizing TI currents in the neocortical regions. The results confirmed that the multipair TIS provides better focalized stimulation than the conventional two-pair TIS for all three head models. It is expected that the proposed multipair TIS can be used to enhance the effectiveness of noninvasive deep brain stimulation.

Finite element contact stress analysis of the temporomandibular joints of patients with temporomandibular disorders under mastication

BACKGROUND AND OBJECTIVES

Temporomandibular disorders (TMD) represent a group of diseases occurred in the temporomandibular joint (TMJ) and its surrounding tissues. In epidemiological studies, up to 75% of adults have shown at least one sign of temporomandibular disorders during their examinations. The incongruous biomechanical environment in the TMJ is the main pathogenic factor of TMD. This study attempts to determine the mechanical differences in different groups of TMD patients through biomechanics and to explain the mechanical pathogenesis of TMD according to various cases.

METHODS

Eleven control subjects and eleven TMD patients were selected and divided into three groups: the control group, bilateral TMD group, and unilateral TMD group. The contact stresses of the articular discs, condyles and temporal bones were analyzed using finite element method and statistical analysis of variance.

RESULTS

The results indicated that the contact stresses in the joints with TMD were significantly greater in the Bi-Group (Bilateral TMD patients) compared to the Control-Group. The TMD side always exhibited greater stresses in the Uni-Group [Unilateral TMD patients) under various conditions (clenching on the TMD side or asymptomatic side). The greatest stress of all the groups occurred at the contralateral side with TMD when clenching on the asymptomatic side.

CONCLUSIONS

Excessive protection would lead to greater stress on the affected side and increased TMD risk on the asymptomatic side. Clinically, the abnormal stress distributions of the disc represented poor buffering and articular clicking. The asymmetric distributions of the articular fossa manifested the deviation of mouth opening or inconsistent TMJ loading.

Reconstructing subcortical and cortical somatosensory activity via the RAMUS inverse source analysis technique using median nerve SEP data

This study concerns reconstructing brain activity at various depths based on non-invasive EEG (electroencephalography) scalp measurements. We aimed at demonstrating the potential of the RAMUS (randomized multiresolution scanning) technique in localizing weakly distinguishable far-field sources in combination with coinciding cortical activity. As we have shown earlier theoretically and through simulations, RAMUS is a novel mathematical method that by employing the multigrid concept, allows marginalizing noise and depth bias effects and thus enables the recovery of both cortical and subcortical brain activity. To show this capability with experimental data, we examined the 14-30 ms post-stimulus somatosensory evoked potential (SEP) responses of human median nerve stimulation in three healthy adult subjects. We aim at reconstructing the different response components by evaluating a RAMUS-based estimate for the primary current density in the nervous tissue. We present source reconstructions obtained with RAMUS and compare them with the literature knowledge of the SEP components and the outcome of the unit-noise gain beamformer (UGNB) and standardized low-resolution brain electromagnetic tomography (sLORETA). We also analyzed the effect of the iterative alternating sequential technique, the optimization technique of RAMUS, compared to the classical minimum norm estimation (MNE) technique. Matching with our previous numerical studies, the current results suggest that RAMUS could have the potential to enhance the detection of simultaneous deep and cortical components and the distinction between the evoked sulcal and gyral activity.

Individually optimized multi-channel tDCS for targeting somatosensory cortex

OBJECTIVE

Transcranial direct current stimulation (tDCS) is a non-invasive neuro-modulation technique that delivers current through the scalp by a pair of patch electrodes (2-Patch). This study proposes a new multi-channel tDCS (mc-tDCS) optimization method, the distributed constrained maximum intensity (D-CMI) approach. For targeting the P20/N20 somatosensory source at Brodmann area 3b, an integrated combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis is used with individualized skull conductivity calibrated realistic head modeling.

METHODS

Simulated electric fields (EF) for our new D-CMI method and the already known maximum intensity (MI), alternating direction method of multipliers (ADMM) and 2-Patch methods were produced and compared for the individualized P20/N20 somatosensory target for 10 subjects.

RESULTS

D-CMI and MI showed highest intensities parallel to the P20/N20 target compared to ADMM and 2-Patch, with ADMM achieving highest focality. D-CMI showed a slight reduction in intensity compared to MI while reducing side effects and skin level sensations by current distribution over multiple stimulation electrodes.

CONCLUSION

Individualized D-CMI montages are preferred for our follow up somatosensory experiment to provide a good balance between high current intensities at the target and reduced side effects and skin sensations.

SIGNIFICANCE

An integrated combined MEG and EEG source analysis with D-CMI montages for mc-tDCS stimulation potentially can improve control, reproducibility and reduce sensitivity differences between sham and real stimulations.

Biomechanical analysis of the cervical spine segment as a method for studying the functional and dynamic anatomy of the human neck

BACKGROUND

Traditionally, dynamic and functional anatomy, in particular the dynamic anatomy of the neck, is studied on cadaveric material. However, the development of in vivo visualization technologies and in silico modeling has made it possible to expand these possibilities. Despite significant progress in the study of dynamic and functional anatomy of the neck by means of in silico methods, the issues of validating the developed models and taking into account the pronounced nonlinearity of soft tissues as well as local anisotropy remain open. The aim of this study was to develop a virtual dynamic anatomical model of the human neck and reproduce the dynamic processes in the cervical spine from this model using the finite element method.

MATERIALS AND METHODS

Reverse engineering was used to generate a dynamic anatomical model of the neck from CT data (both male, 24 and 22 years old). Two segments of the cervical spine (C3-C5, C2-T1) were isolated from the resulting model for finite element analysis. Finite element mesh generation and contact interactions were performed using the HyperMesh software (Altair Engineering Inc, Troy, Michigan, USA). The anisotropic hyperelastic Holzapfel-Gasser-Ogden model was used to describe the material behavior of the fibrous rings of the disc. Material modeling and finite element analysis were performed using Abaqus CAE 6.14 software (Simulia, Johnston, Rhode Island, USA).

RESULTS

A technique for creating a virtual dynamic anatomical model of the neck was elaborated and implemented. The model includes 79 major anatomical structures of the neck segmented from radiological data. A finite element analysis of the cervical spine was performed. The results of finite element analysis of the C3-C5 segment under axial load were compared with in vitro data. The proposed model shows nonlinear deformation of the disc under static loading; the model predicted displacement values agree well with the experimental ones. The displacement of the С3-С5 central vertebra with an axial load of 800 N reaches a value of 0.65 mm. For the segment C2-T1, data on intradiscal pressure, stress plots and displacements during flexion were obtained. The maximum stress value of 10.036 MPa is observed in the C3-C4 disc.

CONCLUSION

Simulation results using the proposed methodology are in good agreement with experimental data. The generated biomechanical models allow describing dynamic phenomena in the cervical spine and obtaining a wide range of quantitative properties of anatomical objects, which are otherwise inaccessible to classical methods for studying dynamic and functional anatomy.

Hemodynamic alternations following stent deployment and post-dilation in a heavily calcified coronary artery: In silico and ex-vivo approaches

In this work, hemodynamic alterations in a patient-specific, heavily calcified coronary artery following stent deployment and post-dilations are quantified using in silico and ex-vivo approaches. Three-dimensional artery models were reconstructed from OCT images. Stent deployment and post-dilation with various inflation pressures were performed through both the finite element method (FEM) and ex vivo experiments. Results from FEM agreed very well with the ex-vivo measurements, interms of lumen areas, stent underexpansion, and strut malapposition. In addition, computational fluid dynamics (CFD) simulations were performed to delineate the hemodynamic alterations after stent deployment and post-dilations. A pressure time history at the inlet and a lumped parameter model (LPM) at the outlet were adopted to mimic the aortic pressure and the distal arterial tree, respectively. The pressure drop across the lesion, pertaining to the clinical measure of instantaneous wave-free flow ratio (iFR), was investigated. Results have shown that post-dilations are necessary for the lumen gain as well as the hemodynamic restoration towards hemostasis. Malapposed struts induced much higher shear rate, flow disturbances and lower time-averaged wall shear stress (TAWSS) around struts. Post-dilations mitigated the strut malapposition, and thus the shear rate. Moreover, stenting induced larger area of low TAWSS (<0.4 Pa) and lager volume of high shear rate (>2000 s^-1), indicating higher risks of in-stent restenosis (ISR) and stent thrombosis (ST), respectively. Oscillatory shear index (OSI) and relative residence time (RRT) indicated the wall regions more prone to ISR are located near the malapposed stent struts.

Polymer-Based Biocompatible Packaging for Implantable Devices: Packaging Method, Materials, and Reliability Simulation

Polymer materials attract more and more interests for a biocompatible package of novel implantable medical devices. Medical implants need to be packaged in a biocompatible way to minimize FBR (Foreign Body Reaction) of the implant. One of the most advanced implantable devices is neural prosthesis device, which consists of polymeric neural electrode and silicon neural signal processing integrated circuit (IC). The overall neural interface system should be packaged in a biocompatible way to be implanted in a patient. The biocompatible packaging is being mainly achieved in two approaches; (1) polymer encapsulation of conventional package based on die attach, wire bond, solder bump, etc. (2) chip-level integrated interconnect, which integrates Si chip with metal thin film deposition through sacrificial release technique. The polymer encapsulation must cover different materials, creating a multitude of interface, which is of much importance in long-term reliability of the implanted biocompatible package. Another failure mode is bio-fluid penetration through the polymer encapsulation layer. To prevent bio-fluid leakage, a diffusion barrier is frequently added to the polymer packaging layer. Such a diffusion barrier is also used in polymer-based neural electrodes. This review paper presents the summary of biocompatible packaging techniques, packaging materials focusing on encapsulation polymer materials and diffusion barrier, and a FEM-based modeling and simulation to study the biocompatible package reliability.

Mechanical performances of balloon post-dilation for improving stent expansion in calcified coronary artery: Computational and experimental investigations

Stent deployment in a calcified coronary artery is often associated with suboptimal outcomes such as stent underexpansion and malapposition. Post-dilation after stent deployment is commonly used for optimal stent implantation. There is no guideline for choosing the post-dilation balloon diameter and inflation pressure. In this work, ex-vivo/in-silico experiments were performed to investigate the efficacy of post-dilation balloon diameter and inflation pressure in improving the stent expansion in a calcified lesion. Post-dilations with three balloon diameters (3 mm, 3.5 mm, and 4 mm) were performed. For each balloon diameter, three inflation pressures (10 atm, 20 atm, and 30 atm) were sequentially applied. In ex-vivo experiments, optical coherence tomography images were acquired during the stenting procedure, i.e., pre- and post-deployment of 3 mm diameter stent, as well as after each post-dilation. The results from in-silico experiments were compared with ex-vivo experiments in terms of lumen area. In addition, stretch ratio analysis was developed to predict the stent-induced lumen area, along with the strain analysis and the in-silico experiments. Results have shown that target lumen area could be achieved with an oversized nominal balloon diameter of +0.5 mm (i.e., 0.5 mm greater than reference lumen diameter) at an inflation pressure of 20 atm. After each post-dilation, fibrotic tissue demonstrated a larger strain, contributing to improved lumen gain. However, minimal changes were observed in calcification. Moreover, a strong correlation (R² = 0.95) between the stretch ratio of fibrotic tissue and lumen area after each post-dilation was observed. This indicated that the morphology of the fibrotic tissue could be a potential marker to predict the lumen gain. The detailed mechanistic quantifications of a single lesion cannot be generalized to all clinical cases. However, this work could be used to provide a fundamental understanding of the post-dilations, to develop experimental protocols for producing generalized guidelines, and to exploit their potential for optimal pre- and post-stent strategies.

Analysis and improvement of positioning reliability and accuracy of theta pipette configuration for scanning ion conductance microscopy

Scanning ion conductance microscopy (SICM) as an emerging non-contact scanning probe microscopy technique and featuring its strong in-situ detectability for soft and viscous samples, is increasingly used in biomedical and materials related studies. In SICM measurements, employing theta pipette as SICM probe to scan sample is an effective method to extend the applications of SICM for multi-parameter measurement. There are two crucial but still unclear issues that influence the reliability and accuracy of the usage of theta pipette in the SICM measurements, which are the safe feedback threshold and the horizontal measurement offset. In this work, aiming at the theta pipette configuration of SICM, we systematically investigated the two issues of the theta pipette by both finite element method (FEM) simulation and SICM experiments. The FEM analysis results show that the safe feedback threshold of the one side barrel of the theta pipette is above 99.5%, and the horizontal measurement offset is ~0.53 times of the inner radius of the probe tip. Based on this, we proposed an improved scanning method used by the theta pipette to solve the reliability and accuracy problems caused by the feedback threshold too close to the reference current (100%) and the measurement offset error at the tip radius level. Then through testing the polydimethylsiloxane (PDMS) samples with different embossed patterns with the improved method of SICM, we can conclude that the improved method can enhance the scanning reliability by adding the double barrels approaching process and increase the positioning accuracy by compensating an offset distance. The theoretical analysis and the improved scanning method in this work demonstrate more property and usage details of the theta pipette, and further improve the reliability and accuracy of the diversified multifunctional applications of the theta pipette for SICM to meet the increasingly complex and precise research needs.

Thermo-mechano-chemical modeling and computation of thermosetting polymers used in post-installed fastening systems in concrete structures

As thermoset polymers find frequent implementation in engineering design, their application in structural engineering is rather limited. One key reason relies on the ongoing curing process in typical applications such as post-installed adhesive anchors, joints by structural elements or surface-mounted laminates glued by adhesive polymers. Mechanochemistry including curing and aging under thermal as well as mechanical loading causes a multiphysics problem to be discussed. For restricting the variety of material models based on empirical observations, we aim at a thermodynamically sound strategy for modeling thermosets. By providing a careful analysis and clearly identifying the assumptions and simplifications, we present the general framework for modeling and computational implementation of thermo-mechano-chemical processes by using open-source codes.

Bio-inspired wideband antenna for wireless applications based on perturbation technique

The evolution of advancement in communication technologies and ever-increasing demand by users for compact communication devices has necessitated a shift in the design approach to achieve antenna structures that are compact and robust. Owing to the diverse communication requirements, antenna systems operating across wide bands have become a necessity. An antenna that is capable of working effectively in several bands is called wideband antenna. In this work, a bio-inspired microstrip antenna (Bi-MPA) for wideband application is proposed and simulated. The radiating patch of the proposed Bi-MPA is the shape of Carica Papaya leaf. The structure was realized through the perturbation of the circular shape patch. The proposed antenna has an impedance bandwidth of 4.3 GHz (1.9 GHz–6.2 GHz) at a return loss of 10 dB while it exhibits a narrow band at 7.2 GHz (6.99–7.44 GHz) and 9.3 GHz (9.15–9.35 GHz) bands. The gain of the proposed antenna is between 2.60 dB and 10.22 dB and the radiation pattern is quasi-omnidirectional. The proposed Bi-MPA is compact and suitable for global system for mobile communication (GSM1900), Universal Mobile Telecommunication System (UMTS), Wireless Local Area Network (WLAN), Long Term Evolution (LTE2300 and LTE2600), Worldwide Interoperability for Microwave Access (WiMAX), C-band, X-band, and sub6 GHz fifth-generation (5G) band. Our contribution to the scientific community in this work is that we have proposed a single antenna structure that is suitable for communication in all the bands mentioned in order to ensure compactness in the mobile devices as compared to base station antennas.

Modeling radiofrequency responses of realistic multi-electrode leads containing helical and straight wires

Purpose

To present a modeling workflow for the evaluation of a lead electromagnetic model (LEM) consisting of a transfer function (TF) and a calibration factor. The LEM represents an analytical relationship between the RF response of a lead and the incident electromagnetic field. The study also highlights the importance of including key geometric details of the lead and the electrode when modeling multi-electrode leads.

Methods

The electrical and thermal responses of multi-electrode leads with helical and straight wires were investigated using 3D electromagnetic (EM) and thermal co-simulations. The net dissipated power (P) around each lead electrode and the net temperature increase at the electrodes (ΔT) were obtained for a set of incident EM fields with different spatial distributions. A reciprocity approach was used to determine a TF for each electrode based on the results of the computational model. The evaluation of the calibration factors and the TF validation were performed using the linear regression of P versus the LEM predictions.

Results

P and ΔT were investigated for four multi-electrode leads and four single-electrode leads containing either helical or straight wires. All electrodes of the multi-electrode lead were found to be points of high power deposition and temperature rise. The LEMs for the individual electrodes varied substantially. A significant dependence of the calibration factors on the surrounding tissue medium was also found. Finally, the model showed that the TF, the calibration factor, P and ΔT for multi-electrode leads differ significantly from those for single-electrode leads.

Conclusion

These results highlight the need to evaluate a LEM for each electrode of a multi-electrode lead as well as for each possible surrounding medium. It is also shown that the results derived from simulations based on simplified single-electrode leads can significantly mislead multi-electrode lead analyses.

Modeling radio-frequency energy-induced heating due to the presence of transcranial electric stimulation setup at 3T

PURPOSE

The purpose of the present study was to develop a numerical workflow for simulating temperature increase in a high-resolution human head and torso model positioned in a whole-body magnetic resonance imaging (MRI) radio-frequency (RF) coil in the presence of a transcranial electric stimulation (tES) setup.

METHODS

A customized human head and torso model was developed from medical image data. Power deposition and temperature rise (ΔT) were evaluated with the model positioned in a whole-body birdcage RF coil in the presence of a tES setup. Multiphysics modeling at 3T (123.2 MHz) on unstructured meshes was based on RF circuit, 3D electromagnetic, and thermal co-simulations. ΔT was obtained for (1) a set of electrical and thermal properties assigned to the scalp region, (2) a set of electrical properties of the gel used to ensure proper electrical contact between the tES electrodes and the scalp, (3) a set of electrical conductivity values of skin tissue, (4) four gel patch shapes, and (5) three electrode shapes.

RESULTS

Significant dependence of power deposition and ΔT on the skin's electrical properties and electrode and gel patch geometries was observed. Differences in maximum ΔT (> 100%) and its location were observed when comparing the results from a model using realistic human tissue properties and one with an external container made of acrylic material. The electrical and thermal properties of the phantom container material also significantly (> 250%) impacted the ΔT results.

CONCLUSION

Simulation results predicted that the electrode and gel geometries, skin electrical conductivity, and position of the temperature sensors have a significant impact on the estimated temperature rise. Therefore, these factors must be considered for reliable assessment of ΔT in subjects undergoing an MRI examination in the presence of a tES setup.

Comparison of magnetic field distributions generated by various permanent magnets for transcranial static magnetic stimulation: A simulation study

Recent experimental studies have shown that static magnetic field can be effective in modulating human brain functions. Following this discovery, a new noninvasive brain stimulation technique was developed: the transcranial static magnetic stimulation (tSMS). Various types of permanent magnets have been used in previous experimental studies, with the aim of validating the effectiveness of tSMS; nevertheless, the spatial distributions of magnetic field generated by these permanent magnets have not been fully investigated. In this study, we compared the distributions of magnetic field on the human cortical surface generated by five different cylindrical magnets (of various dimensions), using the finite element method. Our simulation results demonstrated that the magnitude of magnetic flux density induced in the cortical grey matter of the human brain is proportional to the volume of permanent magnets used, while the magnetic field gradient is not necessarily proportional to the volume of the magnets. Additionally, we showed that the use of magnets with internal holes might not be advantageous. The differences in magnetic field properties induced by various types of permanent magnets suggested that their careful selection, based on magnetic field simulations, might be necessary to increase the effectiveness of tSMS.

Finite element analysis of occlusal splint therapy in patients with bruxism

Background

Bruxism is among the habits considered generally as contributory factors for temporomandibular joint (TMJ) disorders and its etiology is still controversial.

Methods

Three-dimensional models of maxilla and mandible and teeth of 37 patients and 36 control subjects were created using in-vivo image data. The maximum values of stress and deformation were calculated in 21 patients six months after using a splint and compared with those in the initial conditions.

Results

The maximum stresses in the jaw bone and head of mandible were respectively 4.4 and 4.1 times higher in patients than in control subjects. Similar values for deformation were 5.8 and 4.9, respectively. The maximum stress in the jaw bone and head of mandible decreased six months after splint application by up to 71.0 and 72.8%, respectively. Similar values for the maximum deformation were 80.7 and 78.7%, respectively. Following the occlusal splint therapy, the approximation of maximum deformation to the relevant values in control subjects was about 2.6 times the approximation of maximum stress to the relevant values in control subjects. The maximum stress and maximum deformation occurred in all cases in the head of the mandible and the splint had the highest effectiveness in jaw bone adjacent to the molar teeth.

Conclusions

Splint acts as a stress relaxer and dissipates the extra stresses generated as well as the TMJ deformation and deviations due to bruxism. The splint also makes the bilateral and simultaneous loading possible and helps with the treatment of this disorder through regulation of bruxism by creating a biomechanical equilibrium between the physiological loading and the generated stress.

Enhancement of local surface plasmon resonance (LSPR) effect by biocompatible metal clustering based on ZnO nanorods in Raman measurements

The development of size-selective and non-destructive detection techniques for nanosized biomarkers has many reasons, including the study of living cells and diagnostic applications. We present an approach for Raman signal enhancement on biocompatible sensing chips based on surface enhancement Raman spectroscopy (SERS). A sensing chip was fabricated by forming a ZnO-based nanorod structure so that the Raman enhancement occurred at a gap of several tens to several hundred nanometers. The effect of coffee-ring formation was eliminated by introducing the porous ZnO nanorods for the bio-liquid sample. A peculiarity of this approach is that the gold sputtered on the ZnO nanorods initially grows at their heads forming clusters, as confirmed by secondary electron microscopy. This clustering was verified by finite element analysis to be the main factor for enhancement of local surface plasmon resonance (LSPR). This clustering property and the ability to adjust the size of the nanorods enabled the signal acquisition points to be refined using confocal based Raman spectroscopy, which could be applied directly to the sensor chip based on the optimization process in this experiment. It was demonstrated by using common cancer cell lines that cell growth was high on these gold-clad ZnO nanorod-based surface-enhanced Raman substrates. The porosity of the sensing chip, the improved structure for signal enhancement, and the cell assay make these gold-coated ZnO nanorods substrates promising biosensing chips with excellent potential for detecting nanometric biomarkers secreted by cells.

Inherent physiological artifacts in EEG during tDCS

Online imaging and neuromodulation is invalid if stimulation distorts measurements beyond the point of accurate measurement. In theory, combining transcranial Direct Current Stimulation (tDCS) with electroencephalography (EEG) is compelling, as both use non-invasive electrodes and image-guided dose can be informed by the reciprocity principle. To distinguish real changes in EEG from stimulation artifacts, prior studies applied conventional signal processing techniques (e.g. high-pass filtering, ICA). Here, we address the assumptions underlying the suitability of these approaches. We distinguish physiological artifacts - defined as artifacts resulting from interactions between the stimulation induced voltage and the body and so inherent regardless of tDCS or EEG hardware performance - from methodology-related artifacts - arising from non-ideal experimental conditions or non-ideal stimulation and recording equipment performance. Critically, we identify inherent physiological artifacts which are present in all online EEG-tDCS: 1) cardiac distortion and 2) ocular motor distortion. In conjunction, non-inherent physiological artifacts which can be minimized in most experimental conditions include: 1) motion and 2) myogenic distortion. Artifact dynamics were analyzed for varying stimulation parameters (montage, polarity, current) and stimulation hardware. Together with concurrent physiological monitoring (ECG, respiration, ocular, EMG, head motion), and current flow modeling, each physiological artifact was explained by biological source-specific body impedance changes, leading to incremental changes in scalp DC voltage that are significantly larger than real neural signals. Because these artifacts modulate the DC voltage and scale with applied current, they are dose specific such that their contamination cannot be accounted for by conventional experimental controls (e.g. differing stimulation montage or current as a control). Moreover, because the EEG artifacts introduced by physiologic processes during tDCS are high dimensional (as indicated by Generalized Singular Value Decomposition- GSVD), non-stationary, and overlap highly with neurogenic frequencies, these artifacts cannot be easily removed with conventional signal processing techniques. Spatial filtering techniques (GSVD) suggest that the removal of physiological artifacts would significantly degrade signal integrity. Physiological artifacts, as defined here, would emerge only during tDCS, thus processing techniques typically applied to EEG in the absence of tDCS would not be suitable for artifact removal during tDCS. All concurrent EEG-tDCS must account for physiological artifacts that are a) present regardless of equipment used, and b) broadband and confound a broad range of experiments (e.g. oscillatory activity and event related potentials). Removal of these artifacts requires the recognition of their non-stationary, physiology-specific dynamics, and individualized nature. We present a broad taxonomy of artifacts (non/stimulation related), and suggest possible approaches and challenges to denoising online EEG-tDCS stimulation artifacts.

RF-induced heating in tissue near bilateral DBS implants during MRI at 1.5 T and 3T: The role of surgical lead management

Access to MRI is limited for patients with deep brain stimulation (DBS) implants due to safety hazards, including radiofrequency (RF) heating of tissue surrounding the leads. Computational models provide an exquisite tool to explore the multi-variate problem of RF heating and help better understand the interaction of electromagnetic fields and biological tissues. This paper presents a computational approach to assess RF-induced heating, in terms of specific absorption rate (SAR) in the tissue, around the tip of bilateral DBS leads during MRI at 64MHz/1.5 T and 127 MHz/3T. Patient-specific realistic lead models were constructed from post-operative CT images of nine patients operated for sub-thalamic nucleus DBS. Finite element method was applied to calculate the SAR at the tip of left and right DBS contact electrodes. Both transmit head coils and transmit body coils were analyzed. We found a substantial difference between the SAR and temperature rise at the tip of right and left DBS leads, with the lead contralateral to the implanted pulse generator (IPG) exhibiting up to 7 times higher SAR in simulations, and up to 10 times higher temperature rise during measurements. The orientation of incident electric field with respect to lead trajectories was explored and a metric to predict local SAR amplification was introduced. Modification of the lead trajectory was shown to substantially reduce the heating in phantom experiments using both conductive wires and commercially available DBS leads. Finally, the surgical feasibility of implementing the modified trajectories was demonstrated in a patient operated for bilateral DBS.

A realistic, accurate and fast source modeling approach for the EEG forward problem

The aim of this paper is to advance electroencephalography (EEG) source analysis using finite element method (FEM) head volume conductor models that go beyond the standard three compartment (skin, skull, brain) approach and take brain tissue inhomogeneity (gray and white matter and cerebrospinal fluid) into account. The new approach should enable accurate EEG forward modeling in the thin human cortical structures and, more specifically, in the especially thin cortices in children brain research or in pathological applications. The source model should thus be focal enough to be usable in the thin cortices, but should on the other side be more realistic than the current standard mathematical point dipole. Furthermore, it should be numerically accurate and computationally fast. We propose to achieve the best balance between these demands with a current preserving (divergence conforming) dipolar source model. We develop and investigate a varying number of current preserving source basis elements n (n=1,…,n=5). For validation, we conducted numerical experiments within a multi-layered spherical domain, where an analytical solution exists. We show that the accuracy increases along with the number of basis elements, while focality decreases. The results suggest that the best balance between accuracy and focality in thin cortices is achieved with n=4 (or in extreme cases even n=3) basis functions, while in thicker cortices n=5 is recommended to obtain the highest accuracy. We also compare the current preserving approach to two further FEM source modeling techniques, namely partial integration and St. Venant, and show that the best current preserving source model outperforms the competing methods with regard to overall balance. For all tested approaches, FEM transfer matrices enable high computational speed. We implemented the new EEG forward modeling approaches into the open source duneuro library for forward modeling in bioelectromagnetism to enable its broader use by the brain research community. This library is build upon the DUNE framework for parallel finite elements simulations and integrates with high-level toolboxes like FieldTrip. Additionally, an inversion test has been implemented using the realistic head model to demonstrate and compare the differences between the aforementioned source models.

A neural controller for online laser power adjustment during the heat therapy process in the presence of nanoparticles

The present research evaluated the efficiency of a control approach to control the temperature of a breast tumor mass in the presence of nanoparticles exposed to laser radiation. However, if the radiation is carried out in open loop manner it may result in excessive temperature rise healthy cells that exist in the vicinity of tumor's cells. This may lead to the death of healthy cells. So, using closed loop control methods is necessary to guarantee the preservation of healthy cells during the period of radiation. Therefore, in this study, an artificial neural network was trained as a controller. In other words, the trained neural network adjusted the laser power over a period of time in such a way that the temperature in the center of the tumor reached the desired level with an appropriate temporal behavior. The difference between the real temperature of the tumor and the desired temperature of it is the controller input, while the controller output determined the amount of laser power. The simulation studies were carried out using an appropriate physiological model in the presence of nanoparticles. First, Schrödinger equations were solved followed by the effective mass equation. Afterward the optimum number of nanoparticles to be used in the IR field was calculated. Next, the important electro-optical features related to the nanostructure, such as the absorption continuum and reflection continuum had been calculated. The neural network proposed controller was then evaluated through other simulation studies in the tumor mass model. The results showed a promising performance by the trained artificial neural network in adjusting radiated laser power for the desired temperature increase in the center of a tumor mass.

Construction and modeling of a reconfigurable MRI coil for lowering SAR in patients with deep brain stimulation implants

Post-operative MRI of patients with deep brain simulation (DBS) implants is useful to assess complications and diagnose comorbidities, however more than one third of medical centers do not perform MRIs on this patient population due to stringent safety restrictions and liability risks. A new system of reconfigurable magnetic resonance imaging head coil composed of a rotatable linearly-polarized birdcage transmitter and a close-fitting 32-channel receive array is presented for low-SAR imaging of patients with DBS implants. The novel system works by generating a region with low electric field magnitude and steering it to coincide with the DBS lead trajectory. We demonstrate that the new coil system substantially reduces the SAR amplification around DBS electrodes compared to commercially available circularly polarized coils in a cohort of 9 patient-derived realistic DBS lead trajectories. We also show that the optimal coil configuration can be reliably identified from the image artifact on B₁⁺ field maps. Our preliminary results suggest that such a system may provide a viable solution for high-resolution imaging of DBS patients in the future. More data is needed to quantify safety limits and recommend imaging protocols before the novel coil system can be used on patients with DBS implants.

COMETS2: An advanced MATLAB toolbox for the numerical analysis of electric fields generated by transcranial direct current stimulation

BACKGROUND

Since there is no way to measure electric current generated by transcranial direct current stimulation (tDCS) inside the human head through in vivo experiments, numerical analysis based on the finite element method has been widely used to estimate the electric field inside the head. In 2013, we released a MATLAB toolbox named COMETS, which has been used by a number of groups and has helped researchers to gain insight into the electric field distribution during stimulation. The aim of this study was to develop an advanced MATLAB toolbox, named COMETS2, for the numerical analysis of the electric field generated by tDCS.

NEW METHOD

COMETS2 can generate any sizes of rectangular pad electrodes on any positions on the scalp surface. To reduce the large computational burden when repeatedly testing multiple electrode locations and sizes, a new technique to decompose the global stiffness matrix was proposed.

RESULTS

As examples of potential applications, we observed the effects of sizes and displacements of electrodes on the results of electric field analysis. The proposed mesh decomposition method significantly enhanced the overall computational efficiency.

COMPARISON WITH EXISTING METHODS

We implemented an automatic electrode modeler for the first time, and proposed a new technique to enhance the computational efficiency.

CONCLUSIONS

In this paper, an efficient toolbox for tDCS analysis is introduced (freely available at http://www.cometstool.com). It is expected that COMETS2 will be a useful toolbox for researchers who want to benefit from the numerical analysis of electric fields generated by tDCS.

Bacterial gliding fluid dynamics on a layer of non-Newtonian slime: Perturbation and numerical study

Gliding bacteria are an assorted group of rod-shaped prokaryotes that adhere to and glide on certain layers of ooze slime attached to a substratum. Due to the absence of organelles of motility, such as flagella, the gliding motion is caused by the waves moving down the outer surface of these rod-shaped cells. In the present study we employ an undulating surface model to investigate the motility of bacteria on a layer of non-Newtonian slime. The rheological behavior of the slime is characterized by an appropriate constitutive equation, namely the Carreau model. Employing the balances of mass and momentum conservation, the hydrodynamic undulating surface model is transformed into a fourth-order nonlinear differential equation in terms of a stream function under the long wavelength assumption. A perturbation approach is adopted to obtain closed form expressions for stream function, pressure rise per wavelength, forces generated by the organism and power required for propulsion. A numerical technique based on an implicit finite difference scheme is also employed to investigate various features of the model for large values of the rheological parameters of the slime. Verification of the numerical solutions is achieved with a variational finite element method (FEM). The computations demonstrate that the speed of the glider decreases as the rheology of the slime changes from shear-thinning (pseudo-plastic) to shear-thickening (dilatant). Moreover, the viscoelastic nature of the slime tends to increase the swimming speed for the shear-thinning case. The fluid flow in the pumping (generated where the organism is not free to move but instead generates a net fluid flow beneath it) is also investigated in detail. The study is relevant to marine anti-bacterial fouling and medical hygiene biophysics.

Simulation of ultrasonic and EMAT arrays using FEM and FDTD

This paper presents a method which combines electromagnetic simulation and ultrasonic simulation to build EMAT array models. For a specific sensor configuration, Lorentz forces are calculated using the finite element method (FEM), which then can feed through to ultrasonic simulations. The propagation of ultrasound waves is numerically simulated using finite-difference time-domain (FDTD) method to describe their propagation within homogenous medium and their scattering phenomenon by cracks. Radiation pattern obtained with Hilbert transform on time domain waveforms is proposed to characterise the sensor in terms of its beam directivity and field distribution along the steering angle.

The bending strength of tablets with a breaking line--Comparison of the results of an elastic and a "brittle cracking" finite element model with experimental findings

The aim of this work was to ascertain the influence of the position of the breaking line of bevel-edged tablets in a three-point bending test. Two different brands of commercially available, flat-round, bevel-edged tablets with a single central breaking line were studied. Breaking line positions tested, relative to the upper loading roll, were 0°, 22.5°, 45°, 67.5° and 90°. The breaking line faced either up- or downwards during the test. The practical results were compared with FEM results simulating similar test configurations. Tablets failed mainly across the failure plane, resulting in two tablet halves. An exception to this was found for tablets where the breaking line faced down and was positioned at an angle of 22.5° relative to the loading plane. Here the crack followed the breaking line in the centre of the tablets and only diverged towards the loading plane position at the edges of the tablets. The breaking line facing upwards resulted in a significantly higher tensile strength of the tablets compared to it facing downwards. However, with one exception, the orientation of the breaking line relative to the loading plane appeared not to affect the tensile strength values. A fully elastic FEM model indicated that both the position of the breaking line relative to the loading plane and as to whether the breaking line faced up- or downwards during the bending test would result in considerably different failure loads during practical experiments. The results also suggested that regardless of the breaking line position, when it is facing down crack propagation should start at the outer edges propagating towards the midpoint of the discs until failure occurs. Failure should hence always result in equal tablet halves, whereby the failure plane should coincide with the loading plane. Neither predictions fully reflected the practical behaviour of the tablets. Using a brittle cracking FEM model significantly larger tensile stresses for tablets with the breaking line positioned downwards at 0° or 22.5° relative to the loading plane were still predicted, but the differences between model and experimental values was greatly reduced. The remaining differences are more likely due to the inadequacy of the equation available to calculate the experimental tensile strength values. This equation cannot account for the presence of a breaking line and overestimates the thickness of the loading plane by the depth of the breaking line when in 0° or 22.5° position. If the depth of the breaking line is taken into account, the model predictions and the experimental findings are comparable. Also, in the brittle cracking FEM simulations the predicted crack propagation patterns were similar to those found in the experiments, and the model stress distributions across the lower surfaces were much more homogeneous and streamlined parallel to the loading plane. The brittle cracking model hence reflected the practicalities of the bending test more closely. The findings suggested that with the breaking line facing down fracture should always start in the centre of a tablet at its lower surface, initiated by the breaking line. Due to simultaneous development of larger stresses along the y-axis the tablet should still break into two equal halves along the loading plane, unless the position of the breaking line relative to the loading plane was 22.5°. In this case the tablet would fail by a mixed process, whereby failure would occur mainly along the breaking line, but due to simultaneous crack formation at the lower surface close to the bevel edge parallel to the loading plane the final breaking pattern would deviate from the breaking line about half-way from its centre, as seen in the practical experiments.

An optimization algorithm for individualized biomechanical analysis and simulation of tibia fractures

An algorithmic strategy to determine the minimal fusion area of a tibia pseudarthrosis to achieve mechanical stability is presented. For this purpose, a workflow capable for implementation into clinical routine workup of tibia pseudarthrosis was developed using visual computing algorithms for image segmentation, that is a coarsening protocol to reduce computational effort resulting in an individualized volume-mesh based on computed tomography data. An algorithm detecting the minimal amount of fracture union necessary to allow physiological loading without subjecting the implant to stresses and strains that might result in implant failure is developed. The feasibility of the algorithm in terms of computational effort is demonstrated. Numerical finite element simulations show that the minimal fusion area of a tibia pseudarthrosis can be less than 90% of the full circumferential area given a defined maximal von Mises stress in the implant of 80% of the total stress arising in a complete pseudarthrosis of the tibia.

Theoretical investigations into the influence of the position of a breaking line on the tensile failure of flat, round, bevel-edged tablets using finite element methodology (FEM) and its practical relevance for industrial tablet strength testing

Flat, round tablets may have a breaking ("score") line. Pharmacopoeial tablet breaking load tests are diametral in their design, and industrially used breaking load testers often have automatic tablet feeding systems, which position the tablets between the loading platens of the machine with the breaking lines in random orientation to the applied load. The aim of this work was to ascertain the influence of the position of the breaking line in a diametral compression test using finite element methodology (FEM) and to compare the theoretical results with practical findings using commercially produced bevel-edged, scored tablets. Breaking line test positions at an angle of 0°, 22.5°, 45°, 67.5° and 90° relative to the loading plane were studied. FEM results obtained for fully elastic and elasto-plastic tablets were fairly similar, but they highlighted large differences in stress distributions depending on the position of the breaking line. The stress values at failure were predicted to be similar for tablets tested at an angle of 45° or above, whereas at lower test angles the predicted breaking loads were up to three times larger. The stress distributions suggested that not all breaking line angles would result in clean tensile failure. Practical results, however, did not confirm the differences in the predicted breaking loads, but they confirmed differences in the way tablets broke. The results suggest that it is not advisable to convert breaking loads obtained on scored tablets into tablet tensile strength values, and comparisons between different tablets or batches should carefully consider the orientation of the breaking line with respect to the loading plane, as the failure mechanisms appear to vary.

Pressure-volume curves: revisiting the impact of negative turgor during cell collapse by literature review and simulations of cell micromechanics

The Scholander-Hammel pressure chamber has been used in thousands of papers to measure osmotic pressure, πc , turgor pressure, Pt , and bulk modulus of elasticity, ε, of leaf cells by pressure-volume (PV) curve analysis. PV analysis has been questioned in the past. In this paper we use micromechanical analysis of leaf cells to examine the impact on PV curve analysis of negative turgor in living cells (Pt ). Models predict negative Pt (-0.1 to -1.8 MPa) depending on leaf cell size and shape in agreement with experimental values reported by J. J. Oertli. Modeled PV curves have linear regions even when Pt is quite negative, contrary to the arguments of M.T. Tyree. Negative Pt is totally missed by PV curve analysis and results in large errors in derived πc and Pt but smaller errors in ε. A survey of leaf cell sizes vs habitat (arid, temperate, and rainforest), suggests that the majority of published PV curves result in errors of 0.1-1.8 MPa in derived πc and Pt , whereby the error increases with decreasing cell size. We propose that small cell size in leaves is an ecological adaptation that permits plants to endure negative values of water potential with relatively little water loss.

Numerical simulation of the effect of superparamagnetic nanoparticles on microwave rewarming of cryopreserved tissues

In this study, the microwave rewarming process of cryopreserved samples with embedded superparamagnetic (SPM) nanoparticles was numerically simulated. The Finite Element Method (FEM) was used to calculate the coupling of the electromagnetic field and the temperature field in a microwave rewarming system composed of a cylindrical resonant cavity, an antenna source, and a frozen sample phantom with temperature-dependent properties. The heat generated by the sample and the nanoparticles inside the electromagnetic field of the microwave cavity was calculated. The dielectric properties of the biological tissues were approximated using the Debye model, which is applicable at different temperatures. The numerical results showed that, during the rewarming process of the sample phantom without nanoparticles, the rewarming rate was 29.45°C/min and the maximum temperature gradient in the sample was 3.58°C/mm. If nanoparticles were embedded in the sample, and the cavity power was unchanged, the rewarming rate was 47.76°C/min and the maximum temperature gradient in the sample was 1.64°C/mm. In the presence of SPM nanoparticles, the rewarming rate and the maximum temperature gradient were able to reach 20.73°C/min and 0.68°C/mm at the end of the rewarming under the optimized cavity power setting, respectively. The ability to change these temperature behaviors may prevent devitrification and would greatly diminish thermal stress during the rewarming process. The results indicate that the rewarming rate and the uniformity of temperature distribution are increased by nanoparticles. This could be because nanoparticles generated heat in the sample homogeneously and the time-dependent parameters of the sample improved after nanoparticles were homogeneously embedded within it. We were thus able to estimate the positive effect of SPM nanoparticles on microwave rewarming of cryopreserved samples.