Search

Background The objective of this study was to develop a method to evaluate the effects of an aortic dissection on hemodynamic parameters by conducting a comparison with that of a healthy (nondissected) aorta. Open-source software will be implemented, no proprietary software/application will be used to ensure accessorily and repeatability, in all the data analysis and processing. Computed tomography (CT) images of aortic dissection are used for the model geometry segmentation. Boundary conditions from literature are implemented to computational fluid dynamics (CFD) to analyze the hemodynamic parameters. Methods A numerical simulation model was created by obtaining accurate 3-dimensional geometries of aortae from CT images. In this study, CT images of 8 cases of aortic dissection (Stanford type-A and type-B) and 3 cases of healthy aortae are used for the actual aorta model geometry segmentation. These models were exported into an open-source CFD software, OpenFOAM, where a simplified pulsating flow was simulated by controlling the flow pressure. Ten cycles of the pulsatile flow (0.50 sec/cycle) conditions, totaling 5 sec, were calculated. Results The pressure distribution, wall shear stress (WSS) and flow velocity streamlines within the aorta and the false lumen were calculated and visualized. It was found that the flow velocity and WSS had a high correlation in high WSS areas of the intermittent layer between the true and false lumen. Most of the Stanford type-A dissections in the study showed high WSS, over 38 Pa, at the systole phase. This indicates that the arterial walls in type-A dissections are more likely to be damaged with pulsatile flow. Conclusions Using CFD to estimate localized high WSS areas may help in deciding to treat a type-A or B dissection with a stent graft to prevent a potential rupture.

Coronavirus and its spread all over the world have been the most challenging crisis in 2020. Hospitals are categorized among the most vulnerable centers due to their presumably highest traffic of this virus. In this study, centrifugal isolation of coronavirus is successfully deployed for purifying hospitals’ air using air conditioners and ducts, suggesting an efficient setup. Numerical simulations have been used to testify the proposed setup due to the complexities of using experimental investigation such as high cost and clinical hazards of the airborne SARS-CoV-2 in the air. Results show that a 20-cm pipe with an inlet velocity of 4 m/s constitutes the best choice for the separation and purification of air from the virus. The proposed scalable method also efficiently separates larger particles, but it can separate smaller particles too. Numerical results also suggest installing the air purifying system on the floor of the hospitals’ room for maximum efficiency.

Tracking and understanding data quality, analysis and reproducibility are critical concerns in the biological sciences. This is especially true in genomics where next generation sequencing (NGS) based technologies such as ChIP-seq, RNA-seq and ATAC-seq are generating a flood of genome-scale data. However, such data are usually processed with automated tools and pipelines, generating tabular outputs and static visualisations. Interpretation is normally made at a high level without the ability to visualise the underlying data in detail. Conventional genome browsers are limited to browsing single locations and do not allow for interactions with the dataset as a whole. Multi Locus View (MLV), a web-based tool, has been developed to allow users to fluidly interact with genomics datasets at multiple scales. The user is able to browse the raw data, cluster, and combine the data with other analysis and annotate the data. User datasets can then be shared with other users or made public for quick assessment from the academic community. MLV is publically available at https://mlv.molbiol.ox.ac.uk., Martin Sergeant et al. develop a web-based tool called Multi Locus View that allows researchers to interact with genomics datasets at multiple scales. Users can browse, annotate, combine or analyse raw data on this interface. The user datasets generated in this tool can also be made public. This tool is expected to increase public access to the entirety of genomic data.

Non-anatomical placement may occur during the surgical implantation of the meniscal implant, and its influence on the resulting biomechanics of the knee joint has not been systematically studied. The purpose of this study was to evaluate the biomechanical effects of non-anatomical placement of the meniscal implant on the knee joint during a complete walking cycle. Three-dimensional finite element (FE) analyses of the knee joint were performed, based on the model developed from magnetic resonance images and the loading conditions derived from the gait pattern of a healthy male subject, for the following physiological conditions: (i) knee joint with intact native meniscus, (ii) medial meniscectomized knee joint, (iii) knee joint with anatomically placed meniscal implant, and (iv) knee joint with the meniscal implant placed in four different in vitro determined non-anatomical locations. While the native menisci were modeled using the nonlinear hyperelastic Holzapfel-Gasser-Ogden (HGO) constitutive model, the meniscal implant was modeled using the isotropic hyperelastic neo-Hookean model. Placement of the meniscal implant in the non-anatomical lateral-posterior and lateral-anterior locations significantly increased the peak contact pressure in the medial compartment. Placement of the meniscal implant in non-anatomical locations significantly altered the tibial rotational kinematics and increased the total force acting at the meniscal horns. Results suggest that placement of the meniscal implant in non-anatomical locations may restrain its ability to be chondroprotective and may initiate or accelerate cartilage degeneration. In conclusion, clinicians should endeavor to place the implant as closest as possible to the anatomical location to restore the normal knee biomechanics.

Atherosclerosis is a progressive disease that can significantly reduce blood supply to vital organs, being one of the main causes of death worldwide. In this work, a numerical and experimental study in 3D printed stenotic coronary arteries, considering both steady and pulsatile blood flow conditions, is presented. The results revealed that a degree of stenosis superior to 50% creates disturbed flows downstream of the contraction, with an accented increase in the wall shear stress measurements at the stenosis throat. Finally, the multiphase mixture was investigated and compared with a single-phase modelling, and only slight differences were observed right after the stenosis throat.

Purpose: this study aims to quantify the patient-specific hemodynamics of complex conduit routing configurations of coronary artery bypass grafting (CABG) operation which are specifically suitable for off-pump surgeries. Coronary perfusion efficacy and local hemodynamics of multiple left internal mammary artery (LIMA) with sequential and end-to-side anastomosis are investigated. Using a full anatomical model comprised of aortic arch and coronary artery branches the optimum perfusion configuration in multi-vessel coronary artery stenosis is desired. Methodology: two clinically relevant CABG configurations are created using a virtual surgical planning tool where for each configuration set, the stenosis level, anastomosis distance and angle were varied. A non-Newtonian computational fluid dynamics solver in OpenFOAM incorporated with resistance boundary conditions representing the coronary perfusion physiology was developed. The numerical accuracy is verified and results agreed well with a validated commercial cardiovascular flow solver and experiments. For segmental performance analysis, new coronary perfusion indices to quantify deviation from the healthy scenario were introduced. Results: the first simulation configuration set;-a CABG targeting two stenos sites on the left anterior descending artery (LAD), the LIMA graft was capable of 31 mL/min blood supply for all the parametric cases and uphold the healthy LAD perfusion in agreement with the clinical experience. In the second end-to-side anastomosed graft configuration set;-the radial artery graft anastomosed to LIMA, a maximum of 64 mL/min flow rate in LIMA was observed. However, except LAD, the obtuse marginal (OM) and second marginal artery (m2) suffered poor perfusion. In the first set, average wall shear stress (WSS) were in the range of 4 to 35 dyns/cm(2)for in LAD. Nevertheless, for second configuration sets the WSS values were higher as the LIMA could not supply enough blood to OM and m2. Conclusion: the virtual surgical configurations have the potential to improve the quality of operation by providing quantitative surgical insight. The degree of stenosis is a critical factor in terms of coronary perfusion and WSS. The sequential anastomosis can be done safely if the anastomosis angle is less than 90 degrees regardless of degree of stenosis. The smaller proposed perfusion index value,O(0.04 - 0) x 10(2), enable us to quantify the post-op hemodynamic performance by comparing with the ideal healthy physiological flow., Scientific and Technological Research Council of Turkey (TÜBİTAK); 1003 Priority-Research Program Grant; European Union (EU); Horizon 2020; European Research Council (ERC) Proof of Concept Grant; KidsSurgicalPlan

Transcatheter aortic valve replacement (TAVR) strongly depends on the calcification patterns, which may lead to a malapposition of the stented valve and complication onsets in terms of structure kinematics and paravalvular leakage (PVL). From one anatomical-resembling model of the aortic root, six configurations with different calcific deposits were built. TAVR fluid-structure interaction simulations predicted different outcomes for the different calcifications patterns in terms of the final valve configuration in the implantation site and the PVL estimations. In particular models with deposits along the cups coaptation resulted in mild PVL, while those with deposits along the attachment line in moderate PVL.

Mechanical heart valves (MHVs) are widely used to replace diseased heart valves, but it may suffer from cavitation due to the rapid closing velocity of the leaflets, resulting in the damage of red blood cells and platelets. The aim of this study is to apply computational fluid dynamics (CFD) method to investigate the cavitation in bileaflets mechanical heart valves (BMHVs) and discuss the effects of the conduit and leaflet geometries on cavitation intensity. Firstly, CFD method together with moving-grid technology were applied and validated by comparing with experimental results obtained from other literature. Then the leaflets movement and the flow rate of BMHVs with different conduit geometries and leaflet geometries are compared. At last, the duration time of the saturated vapor pressure and the closing velocity of leaflets at the instant of valve closure were used to represent the cavitation intensity. Larger closing velocity of leaflets at the instant of valve closure means higher cavitation intensity. For BMHVs with different conduit geometries, the conduit with Valsalva sinuses has the maximum cavitation intensity and the straight conduit has the minimum cavitation intensity, but the leaflets cannot reach the fully opened state in a straight conduit. For BMHVs with different leaflet geometries, in order to minimize the cavitation intensity, the leaflets are better to have a large thickness and a small rotational radius. CFD method is a promising method to deal with cavitation in BMHVs, and the closing velocity of leaflets has the same trend with the cavitation intensity. By using CFD method, the effects of the conduit geometry and the leaflet geometry on cavitaion in BMHVs are found out.

Atherosclerotic plaque in the femoral is the leading cause of peripheral artery disease (PAD), the worse consequence of which may lead to ulceration and gangrene of the feet. Numerical studies on fluid-structure interactions (FSI) of atherosclerotic femoral arteries enable quantitative analysis of biomechanical features in arteries. This study aims to investigate the hemodynamic performance and its interaction with femoral arterial wall based on the patient-specific model with multiple plaques (calcified and lipid plaques). Three types of models, calcification-only, lipid-only and calcification-lipid models, are established. Hyperelastic material coefficients of the human femoral arteries obtained from experimental studies are employed for all simulations. Oscillation of WSS is observed in the healthy downstream region in the lipid-only model. The pressure around the plaques in the two-plaque model is lower than that in the corresponding one-plaque models due to the reduction of blood flow domain, which consequently diminishes the loading forces on both plaques. Therefore, we found that stress acting on the plaques in the two-plaque model is lower than that in the corresponding one-plaque models. This finding implies that the lipid plaque, accompanied by the calcified plaque around, might reduce its risk of rupture due to the reduced the stress acting on it.

This work aimed to determine the influence of strengthening or weakening of the transversus abdominis (TrA) muscle on loads in the lumbar spine using musculoskeletal modelling methods. The input kinematic data of two positions (a standing position and a position during a sitting-down task) were angles in the elbow joint (0°;4°), shoulder joint (0°;3°), hip joint (0°;75°) and knee joint (0°;69°) as well as the torso tilt angle (0°;32°). It was shown that a change in the TrA physiological cross section area (PCSA) has a crucial impact on lumbar spine loads (2xTrA PCSA causes a reduction in the force in joint L5-S1 by 11% for a standing position and by 25% for a sitting-down position) and value of intra-abdominal pressure.

Atherosclerosis preferentially occurs near the junction of branching vessels, where blood recirculation tends to occur (Malek et al. in J Am Med Assoc 282(21):2035-2042, 1999, https://doi.org/10.1001/jama.282.21.2035 ). For decades, CFD has been used to predict flow patterns such as separation and recirculation zones in hemodynamic models, but those predictions have rarely been validated with experimental data. In the context of verification and validation (VV), we first conduct a CFD benchmark calculation that reproduces the vortex detection experiments of Karino and Goldsmith (1980) with idealised branching blood vessels (Karino and Goldsmith in Trans. Am. Soc. Artif. Internal Organs 26:500-506, 1980). The critical conditions for the formation of recirculation vortices, the so-called critical Reynolds numbers, are the main parameters for comparison with the experimental data to demonstrate the credibility of the CFD workflow. We then characterise the wall shear stresses and develop a surrogate model for the size of formed vortices.An automated parametric study generating more than 12,000 CFD simulations was performed, sweeping the geometries and flow conditions found in the experiments by Karino and Goldsmith. The flow conditions were restricted to steady-state laminar flow, with a range of inflow Reynolds numbers up to 350, with various flow ratios between the main branch outlet and side branch outlet. The side branch diameter was scaled relative to the main branch diameter, ranging from 1.05/3 to 3/3; and the branching angles ranged in size from [Formula: see text] to [Formula: see text]. Recirculation vortices were detected by the inversion of the velocity vector at certain locations, as well as by the inversion of the wall shear stress (WSS) vector.The CFD simulations demonstrated good agreement with the experimental data on the critical Reynolds numbers. The spatial distributions of WSS on each branch were analysed to identify potential regions of disease. Once a vortex is formed, the size of the vortex increases by the square root of the Reynolds number. The CFD data was fitted to a surrogate model that accurately predicts the vortex size without the need to run computationally more expensive CFD simulations.This benchmark study validates the CFD simulation of vortex detection in idealised branching vessels under comprehensive flow conditions. This work also proposes a surrogate model for the size of the vortex, which could reduce the computational requirements in the studies related to branching vessels and complex vascular systems.

Stenosis in the bifurcated coronary artery has attracted wide attention among the researchers. Many investigations have been carried out by means of Computational Fluid Dynamics (CFD) to better understand the physical mechanism inside the stenotic bifurcated artery. However, the main focus of the existing publications is limited to the variation of hemodynamic parameters affected by the stenosis and bifurcation structure. The present study aims to make further evaluations of stenosis development and diagnostic lesion assessments based on the critical values of hemodynamic and practical diagnostic parameters. The bifurcated coronary artery with initial stenosis source has been studied in 2 D unsteady model. Different locations of initial stenosis source have been found to greatly affect the orientation of the further stenosis development. In addition, different stenosis severities (diameter stenosis: 30%, 50% and 70%), different stenosis lengths (3 mm, 6 mm and 9 mm) and different pulse rates (75 bpm, 100 bpm and 120 bpm) as controlling parameters have been investigated.

This paper is focused on specification of conditions at failure in bovine cortical bone. Both experimental and analytical studies are conducted. The experimental part includes a series of novel direct shear tests which examine the sensitivity of shear strength to the applied normal stress for different orientations of the sample microstructure. These experiments are supplemented by standard axial compression and tension tests in order to define and quantify a general form of failure criterion. The analytical part examines two different methodologies, viz. critical plane approach and microstructure tensor approach, for defining the anisotropic strength criterion. A procedure for identification of material parameters is outlined which is based on the results of the performed material tests.

Transcatheter aortic valve replacement (TAVR) is a minimally invasive procedure that provides an effective alternative to open-heart surgical valve replacement for treating advanced calcific aortic valve disease patients. However, complications, such as valve durability, device migration, paravalvular leakage (PVL), and thrombogenicity may lead to increased overall post-TAVR morbidity and mortality. A series of numerical studies involving a self-expandable TAVR valve were performed to evaluate these complications. Structural studies were performed with finite element (FE) analysis, followed by computational fluid dynamics (CFD) simulations, and fluid-structure interaction (FSI) analysis. The FE analysis was utilized to study the effect of TAVR valve implantation depth on valve anchorage in the Living Heart Human Model, which is capable of simulating beating heart during repeated cardiac cycles. The TAVR deployment cases where no valve migration was observed were then used to calculate the post-deployment thrombogenic potential via CFD simulations. FSI analysis followed to further assess the post-deployment TAVR hemodynamic performance for different implantation depths. The deployed valves PVL, geometric and effective orifice areas, and the leaflets structural and flow stress magnitudes were compared to determine the device optimal landing zone. The combined structural and hemodynamic analysis indicated that with the TAVR valve deployed at an aft ventricle position an optimal performance was achieved in the specific anatomy studied. Given the TAVR’s rapid expansion to younger lower-risk patients, the comprehensive numerical methodology proposed here can potentially be used as a predictive tool for both procedural planning and valve design optimization to minimize the reported complications.

Multidimensional correlation magnetic resonance imaging (MRI) is an emerging imaging modality that is capable of disentangling highly heterogeneous and opaque systems according to chemical and physical interactions of water within them. Using this approach, the conventional three dimensional MR scalar images are replaced with spatially resolved multidimensional spectra. The ensuing abundance in microstructural and chemical information is a blessing that incorporates a real challenge: how does one distill and refine it into images while retaining its significant components? In this paper we introduce a general framework that preserves the spectral information from spatially resolved multidimensional data. Equal weight is given to significant spectral components at the single voxel level, resulting in a summarized image spectrum. This spectrum is then used to define spectral regions of interest that are utilized to reconstruct images of sub-voxel components. Using numerical simulations we first show that, contrary to the conventional approach, the proposed framework preserves spectral resolution, and in turn, sensitivity and specificity of the reconstructed images. The retained spectral resolution allows, for the first time, to observe an array of distinct $${T}_{1}$$ T 1 −$${T}_{2}$$ T 2 −$$\langle D\rangle $$ ⟨ D ⟩ components images of the human brain. The robustly generated images of sub-voxel components overcome the limited spatial resolution of MRI, thus advancing multidimensional correlation MRI to fulfilling its full potential.

Fracture healing and growth of the bones are highly related to the stress level. Numerical analysis of stresses is the most effective means to determine the stress level, but it usually requires sufficient resolution to ensure an accurate description of geometry features of bones. In this paper, high-fidelity smoothed finite element method (S-FEM) skull models are created using computed tomography (CT) and micro-computed tomography (μCT) images of a juvenile pig skull. The material properties of the heterogeneous bone are modeled by a varying distribution of Young's modulus mapped to each element and smoothing domain to accurately capture the high heterogeneity. Different types of S-FEM models, including node-based, edge-based, and face-based, are developed for this high-fidelity modeling work. It is found that S-FEM has higher accuracy, in terms of displacements, stresses, and strain energy, compared to the traditional finite element method (FEM). Graphical abstract.

Loop-mediated isothermal amplification or LAMP has been identified to be an efficient technology for point-of-care diagnostics. Paper-based LAMP technique has tremendous potential in replacing the existing tube-based technology as the manufacturing cost of a paper-based device is comparatively lower and easy-to-use. LAMP-based paper diagnostic device for Mycobacterium tuberculosis (MTB) detection is of extreme importance as it will help in early and rapid diagnosis of the affected patients. The fabrication of these devices requires assessment of design parameters on the extent of LAMP amplification reaction. Hence, CFD studies would be extremely beneficial from the design perspective. The current work presents an insight into the CFD simulations for LAMP amplification reaction on a porous paper membrane (nitrocellulose membrane). The convection-diffusion-reaction model is solved on a COMSOL Multiphysics 5.0 platform. Studies on effect of pore size, aspect ratio and initial DNA concentration on the extent of DNA amplification reaction have been carried out. The current paper-based technique is effective in detecting a minimum of 5 copies of DNA contrasting the previous semi-quantitative technique which demonstrated the detection of minimum 98 copies. Overall, the simulation results displayed almost 96% enhancement in the DNA amplification rate on paper membrane. Graphical abstract Graphical abstract for the computational study of DNA amplification reaction via LAMP technique on a porous paper membrane.

The sacrum is a large trilateral bone located at the base of the vertebral column serving to transfer the body weight from the trunk to the pelvis and lower extremities. Over the years, an abundance of sacral anatomical divergences has been reported, including numerical and/or morphological variations of sacral entities. The majority of these anatomical alternations has been incidentally identified during radiological investigations, surgical procedures or discovered in anatomical, anthropological and forensic research studies. Throughout international literature, however, there is a scarcity of an integrative recording of all known anatomical variations of the sacrum in a single study. This constitutes the objective of the present paper: to provide an exhaustive systematic review of the relevant literature, as well as to thoroughly describe all the recognized deviations of the sacrum structure, while highlighting the aspects of their clinical significance.

In-stent restenosis (ISR) is related to local haemodynamics in the arteries after stent intervention. However, the haemodynamics of stents implanted into tapered vessels is rarely studied and remains unclear. This study aimed to study the haemodynamic performance of a stent in a tapered artery to reveal the haemodynamic differences between tapered and cylindrical stents after stent implantation and guide the stent selection for the treatment of coronary artery stenosis. Cylindrical and tapered stents were implanted into the tapered arteries. A model of a cylindrical stent implanted into a cylindrical artery was established as the contrast model. Using the finite element method, the flow velocity and wall shear stress distribution of the three models were compared. At t1, t2, t3 and t4, the flow rate of the tapered artery with tapered stents (TT) after the implantation increased by 8.59, 3.80, 12.81 and 3.66%, respectively. In addition, the wall shear stress in the tapered arteries of TT was 23.48, 36.67, 13.00 and 8.06% higher than that of the tapered arteries with cylindrical stents (TC). The implantation of a tapered stent in the tapered artery can effectively improve intravascular haemodynamics. The tapered stent allows the tapered artery to obtain better haemodynamics and reduces the probability of ISR.

PURPOSE: Patient-specific models of the heart can be used to improve the diagnosis of cardiac diseases, but practical application of these models can be impeded by the computational costs and numerical uncertainties of fitting mechanistic models to clinical measurements from individual patients. Reliable and efficient tuning of these models within clinically appropriate error bounds is a requirement for practical deployment in the time-constrained environment of the clinic. METHODS: We developed an optimization framework to tune parameters of patient-specific mechanistic models using routinely-acquired non-invasive patient data more efficiently than manual methods. We employ a hybrid particle swarm and pattern search optimization algorithm, but the framework can be readily adapted to use other optimization algorithms. RESULTS: We apply the proposed framework to tune full-cycle lumped parameter circulatory models using clinical data. We show that our framework can be easily adapted to optimize cross-species models by tuning the parameters of the same circulation model to four canine subjects. CONCLUSIONS: This work will facilitate the use of biomechanics and circulatory cardiac models in both clinical and research environments by ameliorating the tedious process of manually fitting the parameters.

More than eighty percent of pancreatic cancer involves ductal adenocarcinoma with an abundant desmoplastic extracellular matrix surrounding the solid tumor entity. This aberrant tumor microenvironment facilitates a strong resistance of pancreatic cancer to medication. Although various therapeutic strategies have been reported to be effective in mice with pancreatic cancer, they still need to be tested quantitatively in wider animal-based experiments before being applied as therapies. To aid the design of experiments, we develop a cell-based mathematical model to describe cancer progression under therapy with a specific application to pancreatic cancer. The displacement of cells is simulated by solving a large system of stochastic differential equations with the Euler–Maruyama method. We consider treatment with the PEGylated drug PEGPH20 that breaks down hyaluronan in desmoplastic stroma followed by administration of the chemotherapy drug gemcitabine to inhibit the proliferation of cancer cells. Modeling the effects of PEGPH20 + gemcitabine concentrations is based on Green’s fundamental solutions of the reaction–diffusion equation. Moreover, Monte Carlo simulations are performed to quantitatively investigate uncertainties in the input parameters as well as predictions for the likelihood of success of cancer therapy. Our simplified model is able to simulate cancer progression and evaluate treatments to inhibit the progression of cancer. Electronic supplementary material The online version of this article (10.1007/s10237-019-01219-0) contains supplementary material, which is available to authorized users.

Column selection and experimental design are recognized to be wear-and-tear processes to obtain a good fingerprinting result with comprehensive two dimensional gas chromatography (GC × GC). The processes involve a large number of experiments for analysis of each sample due to the optimized conditions depending on the selected column sets. In this study, the analytical solutions of time summation model combined with temperature dependent linear solvation energy relationship (LSER) were derived for constant flow separation in GC × GC. The derived equations allow calculation of analyte retention time in first and second dimensional separation (1tR and 2tR), under temperature program separation employing any column combination with known LSER database. As a result, optimization software was developed enabling simulation of several hundred thousand GC × GC results (fingerprints) for separation of each sample. Good correlations between our predicted results and the results obtained with the previously established numerical approach, were obtained with the R2 of 1.000 and 0.998 for simulation of 1tR and 2tR, respectively. The developed approaches were further applied to simulation of 96,000 individual fingerprints in temperature programmed GC × GC, with the focus on application of 16 column sets including non-ionic liquid and ionic liquid (IL) stationary phases for separation of 678 model compounds. These approaches resulted in the computational time of 1 day compared with 1 year provided by the numerical method. Best column sets and experimental conditions (secondary column lengths and temperature programs) could then be extracted according to maximizing number of separated peaks in separation, which represents the quality of each fingerprinting.

Background Validating new techniques for fetal cardiovascular magnetic resonance (CMR) is challenging due to random fetal movement that precludes repeat measurements. Consequently, fetal CMR development has been largely performed using physical phantoms or postnatal volunteers. In this work, we present an open-source simulation designed to aid in the development and validation of new approaches for fetal CMR. Our approach, fetal extended Cardiac-Torso cardiovascular magnetic resonance imaging (Fetal XCMR), builds on established methods for simulating CMR acquisitions but is tailored toward the dynamic physiology of the fetal heart and body. We present comparisons between the Fetal XCMR phantom and data acquired in utero, resulting in image quality, anatomy, tissue signals and contrast. Methods Existing extended Cardiac-Torso models are modified to create maternal and fetal anatomy, combined according to simulated motion, mapped to CMR contrast, and converted to CMR data. To provide a comparison between the proposed simulation and experimental fetal CMR images acquired in utero, images from a typical scan of a pregnant woman are included and simulated acquisitions were generated using matching CMR parameters, motion and noise levels. Three reconstruction (static, real-time, and CINE), and two motion estimation methods (translational motion, fetal heart rate) from data acquired in transverse, sagittal, coronal, and short-axis planes of the fetal heart were performed to compare to in utero acquisitions and demonstrate feasibility of the proposed simulation framework. Results Overall, CMR contrast, morphologies, and relative proportions of the maternal and fetal anatomy are well represented by the Fetal XCMR images when comparing the simulation to static images acquired in utero. Additionally, visualization of maternal respiratory and fetal cardiac motion is comparable between Fetal XCMR and in utero real-time images. Finally, high quality CINE image reconstructions provide excellent delineation of fetal cardiac anatomy and temporal dynamics for both data types. Conclusion The fetal CMR phantom provides a new method for evaluating fetal CMR acquisition and reconstruction methods by simulating the underlying anatomy and physiology. As the field of fetal CMR continues to grow, new methods will become available and require careful validation. The fetal CMR phantom is therefore a powerful and convenient tool in the continued development of fetal cardiac imaging. Electronic supplementary material The online version of this article (10.1186/s12968-019-0539-2) contains supplementary material, which is available to authorized users.

As a response to the Tokyo Electric Power Company's Fukushima Daiichi nuclear power plant accident in 2011, seven TEPCO workers whose exposure doses were expected to be >250 mSv (a tentative dose limit stipulated by the Japanese central authority) attended Japan's National Institute for Radiological Sciences for additional internal dose measurements. The National Institute for Radiological Sciences examination revealed that these workers' internal doses came mainly from their intake of the radionuclide I during emergency operations. In this study, we performed numerical simulations based on individual volume-pixel (voxel) phantoms of six of the seven workers for a more sophisticated evaluation of their internal doses, taking into account the individual thyroid size and other specific parameters. The voxel phantoms were created from magnetic resonance imaging scan images. As a result, the individual thyroid volumes ranged from 6.5 to 28.2 cm and were considerably smaller than the reference value (~20 cm) adopted in the International Commission on Radiation Protection's dosimetric model for four of the six subjects. Compared to the original estimates of the thyroid absorbed dose, our preliminary evaluation revealed values that were increased by approximately 3-fold or decreased by 30% at maximum. A wide difference in the individual thyroid size would be one of the significant modifiers in the current dose estimation of subjects of the ongoing epidemiological study project. The present simulations also provided evidence that the direct thyroid measurements by the National Institute for Radiological Sciences to determine the workers' I thyroid contents were sufficiently accurate.

T-wave in electrocardiogram (ECG) is a vital wave component and has potential of diagnosing various cardiac disorders. The present work proposes a novel technique for T-wave peak detection using minimal pre-processing and simple root mean square based decision rule. The technique uses a two-stage median filter and a Savitzky–Golay smoothing filter for pre-processing. P-QRS-complex is removed from the filtered ECG, and T-wave is left as the most prominent wave segment, which can be detected using a root mean square based adaptive threshold. An RR-interval based T-wave peak correction strategy has been proposed which can handle the challenges of morphological variations in the T-wave, thus increases the detection accuracy. The proposed technique has been substantiated on a standard QT-database. The detection sensitivity = 97.01%, positive predictivity = 99.61%, detection error rate = 3.36%, and accuracy = 96.66% have been achieved. A T-wave detection technique requiring minimal pre-processing and with simple decision rule has been designed. The noticeably high positive predictivity rate of the proposed technique shows its efficiency to detect T-wave peak.

Hyaline cartilage undergoes many substantial age-related physiochemical and biomechanical changes that reduce its ability to overcome the effects of mechanical stress and injury. In quest of therapeutic options, magnetic stimulation and electrical stimulation (ES) have been proposed for improving tissue engineering approaches for the repair of articular cartilage. The aim of this study is to summarize in silico investigations involving induced electrical properties of cartilage tissue due to various biophysical stimuli along their respective mathematical descriptions. Based on these, a preliminary numerical study involving electromechanical transduction in bovine cartilage tissue has been carried out using an open source finite element computational software. The simulation results have been compared to experimental results from the literature. This study serves as a basis for further in silico studies to better understand the behavior of hyaline cartilage tissue due to ES and to find an optimal stimulation protocol for the cartilage regeneration. Moreover, it provides an overview of the basic models along with mathematical description and scope for future research regarding electrical behavior of the cartilage tissue using open source software. Impact Statement The presented research summarizes the basic models with mathematical description regarding electrical behavior of the cartilage tissue. A preliminary numerical study involving electromechanical transduction in bovine cartilage tissue sample has been carried out using an open source finite element software. This research will provide scope for future research regarding electrical behavior of the cartilage tissue using open source software.

Balloon-occluded transarterial chemoembolisation (B-TACE) is an intraarterial transcatheter treatment for liver cancer. In B-TACE, an artery-occluding microballoon catheter occludes an artery and promotes collateral circulation for drug delivery to tumours. This paper presents a methodology for analysing the haemodynamics during B-TACE, by combining zero-dimensional and three-dimensional modelling tools. As a proof of concept, we apply the methodology to a patient-specific hepatic artery geometry and analyse two catheter locations. Results show that the blood flow redistribution can be predicted in this proof-of-concept study, suggesting that this approach could potentially be used to optimise catheter location.

The number of hip resurfacing arthroplasty procedures has declined dramatically in recent years, for reasons related to the survival rate. Some studies suggest that metal particles are the main critical problem, but do not specify the effect of femoral position on the failure rate. The present study aims to analyze whether the positioning of the resurfacing head implant is important in the distribution of bone strains and in the risk of fracture of the femur. Three in vitro experimental models received the Birmingham hip resurfacing implant to replicate the total hip joint. The resurfacing head of the implanted models was placed in three different offset positions: in a positive offset, with the same femoral head center and in a negative offset. The numerical models were validated by correlating numerical and experimental results. Comparing experimental results from the implanted and intact femurs highlights a strain increase of up to 48% in the proximal medial femur region for positive offset and up to 18% in the neutral position. A reduction of 72% for negative offset (valgus position) was also measured experimentally. A significant change in strain distributions was observed with a resurfacing hip system and increased risk of neck fracture was found using the resurfacing head in positive offset. The iliac bone presents a high decrease in strains that will induce bone loss in the long term. Among the offset positions tested, results suggest that the negative offset (valgus position) and the natural position are the best equilibrated for better long-term results.

The long and binding treatment of tuberculosis (TB) at least 6–8 months for the new cases, the partial immunity given by BCG vaccine, the loss of immunity after a few years doing that strategy of TB control via vaccination and treatment of infectious are not sufficient to eradicate TB. TB is an infectious disease caused by the bacillus Mycobacterium tuberculosis. Adults are principally attacked. In this work, we assess the impact of vaccination in the spread of TB via a deterministic epidemic model (SV ELI) (Susceptible, Vaccinated, Early latent, Late latent, Infectious). Using the Lyapunov–Lasalle method, we analyse the stability of epidemic system (SV ELI) around the equilibriums (disease-free and endemic). The global asymptotic stability of the unique endemic equilibrium whenever $ \mathcal R_{0} \gt 1 $ is proved, where $ \mathcal R_{0} $ is the reproduction number. We prove also that when $ \mathcal R_{0} $ is less than 1, TB can be eradicated. Numerical simulations, using some TB data found in the literature in relation with Cameroon, are conducted to approve analytic results, and to show that vaccination coverage is not sufficient to control TB, effective contact rate has a high impact in the spread of TB.

Maintaining sustainable ecosystems are important for all the inhabitants of earth. Also, an important component of sustainable ecosystems is the maintenance of healthy coexistence of consumers and their resources which can include diseases in the species involved. We formulate a model, where the resources are plants, to explore how consumer-resource coexistence could of itself limit the spread of infectious diseases. The important mathematical features of the model are discussed using the basic reproduction number and the consumption number. The results show an association between species coexistence and a decrease in ecosystem resource disease. The possible importance of these results are discussed.

Over the last few decades, a plethora of tools has been developed for neuroscientists to interface with the brain. Implementing these tools requires precisely removing sections of the skull to access the brain. These delicate cranial microsurgical procedures need to be performed on the sub-millimeter thick bone without damaging the underlying tissue and therefore, require significant training. Automating some of these procedures would not only enable more precise microsurgical operations, but also facilitate widespread use of advanced neurotechnologies. Here, we introduce the “Craniobot”, a cranial microsurgery platform that combines automated skull surface profiling with a computer numerical controlled (CNC) milling machine to perform a variety of cranial microsurgical procedures on mice. The Craniobot utilizes a low-force contact sensor to profile the skull surface and uses this information to perform precise milling operations within minutes. We have used the Craniobot to perform intact skull thinning and open small to large craniotomies over the dorsal cortex.

BACKGROUND Although many studies have reported that the presence of minor or major ST-T change of electrocardiography (ECG) was associated with a risk of cardiovascular events, it is not clear whether there is a difference in the prognostic power depending on the summation of ST-T area (ST-Tarea) assessed by a quantitative method. METHODS Electrocardiograms were performed in 834 clinical patients with one or more cardiovascular risks. ST-Tarea was assessed as the area enclosed by the baseline from the end of the QRS complex to the end of the ST-T segment using a computerized quantitative method. We used the lower magnitude of ST-Tarea in the V5 or V6 lead for the analysis. RESULTS After a mean follow-up 8.4 ± 2.9 years (7,001 person-years), there were 92 cardiovascular events. With adjustment for covariates, the results from Cox proportional hazards models (Model 1) suggested that the lowest quartile of ST-Tarea was associated with a higher risk for cardiovascular outcome compared with the remaining quartile groups (hazard ratio, 2.08; 95% confidence interval, 1.36–3.16, P < 0.01). Even when adding the ECG left ventricular hypertrophy by Cornell voltage (Model 2) and Cornell product (Model 3) to Model 1, the significance remained (both P < 0.01). When we used ST-Tarea as a continuous variable substitute for the lowest quartile of ST-Tarea, these associations were similar in all models (all P < 0.01). CONCLUSION The lower summations of ST-T area assessed by a computerized quantitative method were associated with increased risk of cardiovascular disease incidence in a clinical population.

Central systolic blood pressure (cSBP, the peak of the central waveform) is usually regarded as the determinant of peripheral systolic blood pressure with amplification of peripheral systolic blood pressure (pSBP) measured with reference to cSBP. However, the earlier portion of the central waveform up to the first systolic shoulder (P1) may be the major determinant of pSBP. We performed in silico simulation studies and examined previously acquired experimental data (n = 131) in which peripheral and central blood pressure waveforms had been acquired both invasively and noninvasively to examine the determinants of pSBP. Measurements were made at baseline and during perturbation of hemodynamics by inotropic and vasoactive drugs. In silico simulations using a central-to-peripheral transfer function demonstrated that pSBP is dependent on P1 and the rate of change (dP/dt) of central pressure up to the time of P1 but not cSBP. In computational simulations, peripheral reflection in the radial artery was closely related to dP/dt, and 97% of the variability in amplification as measured with reference to P1 was explained by dP/dt. In vivo, amplification of pSBP over P1 was correlated with dP/dt (R > 0.75, P < 0.0001 for all data sets), and P1 and dP/dt were independently correlated with pSBP, explaining 90% of the variability in pSBP. We conclude that P1 and dP/dt are major determinants of pSBP and that pSBP and cSBP are, in part, determined by different cardiac, central, and peripheral vascular properties.NEW & NOTEWORTHY Peripheral systolic BP is determined mainly by the first shoulder and the rate of rise of the central systolic blood pressure waveform rather than the peak of this waveform (central systolic BP). Peripheral and central systolic blood pressure are determined by different cardiac and vascular properties.

In response to the soaring needs of human mobility data, especially during disaster events such as the COVID-19 pandemic, and the associated big data challenges, we develop a scalable online platform for extracting, analyzing, and sharing multi-source multi-scale human mobility flows. Within the platform, an origin-destination-time (ODT) data model is proposed to work with scalable query engines to handle heterogenous mobility data in large volumes with extensive spatial coverage, which allows for efficient extraction, query, and aggregation of billion-level origin-destination (OD) flows in parallel at the server-side. An interactive spatial web portal, ODT Flow Explorer, is developed to allow users to explore multi-source mobility datasets with user-defined spatiotemporal scales. To promote reproducibility and replicability, we further develop ODT Flow REST APIs that provide researchers with the flexibility to access the data programmatically via workflows, codes, and programs. Demonstrations are provided to illustrate the potential of the APIs integrating with scientific workflows and with the Jupyter Notebook environment. We believe the platform coupled with the derived multi-scale mobility data can assist human mobility monitoring and analysis during disaster events such as the ongoing COVID-19 pandemic and benefit both scientific communities and the general public in understanding human mobility dynamics.

Three-D head geometrical models of eight healthy subjects and 11 hydrocephalus patients were built using their CINE phase-contrast MRI data and used for computer simulations under three different inlet/outlet boundary conditions (BCs). The maximum cerebrospinal fluid (CSF) pressure and the ventricular system volume were more effective and accurate than the other parameters in evaluating the patients’ conditions. In constant CSF pressure, the computational patient models were 18.5% more sensitive to CSF volume changes in the ventricular system under BC “C”. Pulsatile CSF flow rate diagrams were used for inlet and outlet BCs of BC “C”. BC “C” was suggested to evaluate the intracranial compliance of the hydrocephalus patients. The results suggested using the computational fluid dynamic (CFD) method and the fully coupled fluid-structure interaction (FSI) method for the CSF dynamic analysis in patients with external and internal hydrocephalus, respectively., Seifollah Gholampour et al. develop a computational model to examine the flow of cerebrospinal fluid (CSF) in hydrocephalus patients and healthy controls, and simulate how different biophysical parameters can influence CSF dynamics in the brain. Ultimately, their results could be used to better examine the CSF dynamics in a healthy or hydrocephalus brain, without the need for invasive procedures.

CFD has emerged as a promising diagnostic tool for clinical trials, with tremendous potential. However, for real clinical applications to be useful, overall statistical findings from large population samples (e.g., multiple cases and models) are needed. Fully resolved solutions are not a priority, but rather rapid solutions with fast turn-around times are desired. This leads to the issue of what are the minimum modelling criteria for achieving adequate accuracy in respiratory flows for large-scale clinical applications, with a view to rapid turnaround times. This study simulated a highly-resolved solution using the large eddy simulation (LES) method as a reference case for comparison with lower resolution models that included larger time steps and no turbulence modelling. Differences in solutions were quantified by pressure loss, flow resistance, unsteadiness, turbulence intensity, and hysteresis effects from multiple cycles. The results demonstrated that sufficient accuracy could be achieved with lower resolution models if the mean flow was considered. Furthermore, to achieve an established transient result unaffected by the initial start-up quiescent effects, the results need to be taken from at least the second respiration cycle. It was also found that the exhalation phase exhibited strong turbulence. The results are expected to provide guidance for future modelling efforts for clinical and engineering applications requiring large numbers of cases using simplified modelling approaches.

This work deals with the finite element (FE) implementation of a biphasic poroelastic formulation specifically developed to address the intricate behaviour of the Intervertebral Disc (IVD) and other highly hydrated soft tissues. This formulation is implemented in custom FE solver V-Biomech, being the validation performed with a lumbar IVD model, which was compared against the analogous FE model of Williams et al. and the experiments of Tyrrell et al. Good agreement with these benchmarks was achieved, meaning that V-Biomech and its novel poroelastic formulation are a viable alternative for simulation of biphasic soft tissues., The first author would like to express his gratitude to "" (Portugal) for PhD grant SFRH/BD/63882/2009 and LAETA project UIDB/50022/2020 (through IDMEC). Both authors would like to thank Professors Jacques Huyghe (University of Limerick, Republic of Ireland) and Damien Lacroix (University of Sheffield, United Kingdom) for the useful discussions on this topic.

In this paper we address the possibility to perform imaging of two samples within the same acquisition time using coupled ceramic resonators and one transmit/receive channel. We theoretically and experimentally compare the operation of our ceramic dual‐resonator probe with a wire‐wound solenoid probe, which is the standard probe used in ultrahigh‐field magnetic resonance microscopy. We show that due to the low‐loss ceramics used to fabricate the resonators, and a favorable distribution of the electric field within the conducting sample, a dual probe, which contains two samples, achieves an SNR enhancement by a factor close to the square root of 2 compared with a solenoid optimized for one sample., The possibility of simultaneous imaging of two samples is theoretically and experimentally studied using coupled resonators. Our ceramic probe is compared with a compact solenoid probe, the standard probe in ultrahigh‐field magnetic resonance microscopy. Due to low‐loss ceramics and a beneficial distribution of the electric field in a conducting sample, a dual probe containing two samples achieves an SNR enhancement by a factor close to the square root of 2 compared with a solenoid optimized for one sample.

The abdominal aorta is the largest artery in the abdominal cavity that supplies blood flows to vital organs through the complex visceral arterial branches, including the celiac trunk (the liver, stomach, spleen, etc.), the renal arteries (the kidneys) and the superior and inferior mesenteric arteries (the small and large intestine, pancreas, etc.). An accurate simulation of blood flows in this network of arteries is important for the understanding of the hemodynamics in various organs of healthy and diseased patients, but the computational cost is very high. As a result, most researchers choose to focus on a portion of the artery or use a low-dimensional approximation of the artery. In the present work, we introduce a parallel algorithm for the modeling of pulsatile flows in the abdominal aorta with branches to the primary organs, and an organ-based two-level method for calculating the resistances for the outflow boundary conditions. With this highly parallel approach, the simulation of the blood flow for a cardiac cycle of the anatomically detailed aorta can be obtained within a few hours, and the blood distribution to organs including liver, spleen and kidneys are also computed with certain accuracy. Moreover, we discuss the significant hemodynamic differences resulted from the influence of the peripheral branches. In addition, we examine the accuracy of the results with respect to the mesh size and time-step size and show the high parallel scalability of the proposed algorithm with up to 3000 processor cores.

Background The responsibility for helping patients understand potential health benefits and risks, especially regarding screening tests, falls largely to general practitioners (GPs). The Berlin Numeracy Test (BNT) specifically measures risk literacy (i.e., the ability to understand different aspects of statistical numeracy associated with accurate interpretation of information about risks). This study explored the association between risk literacy levels and clinical experience in GPs vs. medical students. Additionally, the effect of GP risk literacy on evaluation of the predictive value of screening tests was examined. Methods The participants were 84 GPs and 92 third-year medical students who completed the BNT (total score range 0–4 points). The GPs received an additional case scenario on mammography screening as a simple measure of performance in applying numeracy skills. Results Despite having an average of 25.9 years of clinical experience, GPs scored no better than medical students on risk literacy (GPs: 2.33 points, 95% confidence interval [CI] 2.08–2.59; students: 2.34, 95% CI 2.07–2.61; P = .983). Of all GPs, 71.6% (n = 58) greatly overestimated the real predictive value. Conclusions In this study, we found no difference in risk literacy between current students and current GPs. GPs lack risk literacy and consequently do not fully understand numeric estimates of probability in routine screening procedures.

Several deep-learning models have been proposed to shorten MRI scan time. Prior deep-learning models that utilize real-valued kernels have limited capability to learn rich representations of complex MRI data. In this work, we utilize a complex-valued convolutional network (ℂNet) for fast reconstruction of highly under-sampled MRI data and evaluate its ability to rapidly reconstruct 3D late gadolinium enhancement (LGE) data. ℂNet preserves the complex nature and optimal combination of real and imaginary components of MRI data throughout the reconstruction process by utilizing complex-valued convolution, novel radial batch normalization, and complex activation function layers in a U-Net architecture. A prospectively under-sampled 3D LGE cardiac MRI dataset of 219 patients (17 003 images) at acceleration rates R = 3 through R = 5 was used to evaluate ℂNet. The dataset was further retrospectively under-sampled to a maximum of R = 8 to simulate higher acceleration rates. We created three reconstructions of the 3D LGE dataset using (1) ℂNet, (2) a compressed-sensing-based low-dimensional-structure self-learning and thresholding algorithm (LOST), and (3) a real-valued U-Net (realNet) with the same number of parameters as ℂNet. LOST-reconstructed data were considered the reference for training and evaluation of all models. The reconstructed images were quantitatively evaluated using mean-squared error (MSE) and the structural similarity index measure (SSIM), and subjectively evaluated by three independent readers. Quantitatively, ℂNet-reconstructed images had significantly improved MSE and SSIM values compared with realNet (MSE, 0.077 versus 0.091; SSIM, 0.876 versus 0.733, respectively; p < 0.01). Subjective quality assessment showed that ℂNet-reconstructed image quality was similar to that of compressed sensing and significantly better than that of realNet. ℂNet reconstruction was also more than 300 times faster than compressed sensing. Retrospective under-sampled images demonstrate the potential of ℂNet at higher acceleration rates. ℂNet enables fast reconstruction of highly accelerated 3D MRI with superior performance to real-valued networks, and achieves faster reconstruction than compressed sensing.

PURPOSE: Modeling of hemolysis due to fluid stresses faces significant methodological challenges, particularly in geometries with turbulence or complex flow patterns. It is currently unclear how existing phenomenological blood-damage models based on laminar viscous stresses can be implemented into turbulent computational fluid dynamics simulations. The aim of this work is to generalize the existing laminar models to turbulent flows based on first principles, and validate this generalization with existing experimental data. METHODS: A novel analytical and numerical framework for the simulation of flow-induced hemolysis based on the intermittency-corrected turbulent viscous shear stress (ICTVSS) is introduced. The proposed large-eddy simulation framework is able to seamlessly transition from laminar to turbulent conditions in a single flow domain by linking laminar shear stresses to dissipation of mechanical energy, accounting for intermittency in turbulent dissipation, and relying on existing power-law hemolysis models. Simulations are run to reproduce previously published hemolysis data with bovine blood in a benchmark geometry. Two sets of experimental data are relied upon to tune power-law parameters and justify that tuning. The first presents hemolysis measurements in a simple laminar flow, and the second is hemolysis in turbulent flow through the FDA benchmark nozzle. Validation is performed by simulation of blood injected into a turbulent jet of phosphate-buffered saline, with modifications made to account for the local concentration of blood. RESULTS: Hemolysis predictions are found to be very sensitive to power-law parameters in the turbulent case, though a set of parameters is presented that both matches the turbulent data and is well-justified by the laminar data. The model is shown to be able to predict the general behavior of hemolysis in a second turbulent case. Results suggest that wall shear may play a dominant role in most cases. CONCLUSION: The ICTVSS framework of generalizing laminar power-law models to turbulent flows shows promise, but would benefit from further numerical validation and carefully designed experiments.

The present study aims to compare spinal stability after two different minimally invasive techniques, the lateral lumbar interbody fusion (LLIF) and the transforaminal lumbar interbody fusion (TLIF) approaches. Two nonlinear three-dimensional finite element (FE) models of the L4-L5 functional spinal unit (FSU) were subjected to the loads that usually act on the lumbar spine. Findings show that the LLIF approach yields better results for torsion load case, due to the larger surface area of the implant. For extension, flexion and lateral bending loads, the TLIF approach presents smaller displacements probably due to the anterior placement of the cage and to the smaller damaged area of the annulus fibrosus.

The purpose of this study is to investigate the hemodynamic significance of various degrees of coronary area of stenosis (AS) and multiple sequential stenoses (MSS) in normal and hypertension pressure conditions. MSS in a single branch coronary artery pose challenges to determine the physiological assessment in the prevalent invasive intervention. The hemodynamic parameters of each stenosis are influenced by other stenoses in the single branch of MSS coronary artery. In this study, we entirely use open source tools and techniques for coronary computed tomography angiography (CCTA) image segmentation, 3D reconstruction, grid generation and hemodynamic simulations. The results yield different hemodynamic parameters such as velocity magnitude, mean arterial pressure difference, flow-pressure linear relation, wall shear stress (WSS) and eventually virtual fractional flow reserve (vFFR) allowing for the prediction and the assessment of lumen area severity conditions in MSS coronaries.

We propose a model of a joint spread of heroin use and HIV infection. The unique disease-free equilibrium always exists and it is stable if the basic reproduction numbers of heroin use and HIV infection are both less than 1. The semi-trivial equilibrium of HIV infection (heroin use) exists if the basic reproduction number of HIV infection (heroin use) is larger than 1 and it is locally stable if and only if the invasion number of heroin use (HIV infection) is less than 1. When both semi-trivial equilibria lose their stability, a coexistence equilibrium occurs, which may not be unique. We compare the model to US data on heroin use and HIV transmission. We conclude that the two diseases in the US are in a coexistence regime. Elasticities of the invasion numbers suggest two foci for control measures: targeting the drug abuse epidemic and reducing HIV risk in drug-users.

Many researchers from a variety of fields including computer science, network science and mathematics have focused on how to contain the outbreaks of Internet misinformation that threaten social systems and undermine societal health. Most research on this topic treats the connections among individuals as static, but these connections change in time, and thus social networks are also temporal networks. Currently there is no theoretical approach to the problem of containing misinformation outbreaks in temporal networks. We thus propose a misinformation spreading model for temporal networks and describe it using a new theoretical approach. We propose a heuristic-containing (HC) strategy based on optimizing final outbreak size that outperforms simplified strategies such as those that are random-containing (RC) and targeted-containing (TC). We verify the effectiveness of our HC strategy on both artificial and real-world networks by performing extensive numerical simulations and theoretical analyses. We find that the HC strategy greatly increases the outbreak threshold and decreases the final outbreak threshold., Comment: 22 pages, 9 figures

Most fish swim with body undulations that result from fluid–structure interactions between the fish’s internal tissues and the surrounding water. Gaining insight into these complex fluid–structure interactions is essential to understand how fish swim. To this end, we developed a dedicated experimental–numerical inverse dynamics approach to calculate the lateral bending moment distributions for a large-amplitude undulatory swimmer that moves freely in three-dimensional space. We combined automated motion tracking from multiple synchronised high-speed video sequences, computation of fluid dynamic stresses on the swimmer’s body from computational fluid dynamics, and bending moment calculations using these stresses as input for a novel beam model of the body. The bending moment, which represent the system’s net actuation, varies over time and along the fish’s central axis due to muscle actions, passive tissues, inertia, and fluid dynamics. Our three-dimensional analysis of 113 swimming events of zebrafish larvae ranging in age from 3 to 12 days after fertilisation shows that these bending moment patterns are not only relatively simple but also strikingly similar throughout early development and from fast starts to periodic swimming. This suggests that fish larvae may produce and adjust swimming movements relatively simply, yet effectively, while restructuring their neuromuscular control system throughout their rapid development., Most fish swim with body undulations that result from fluid-structure interactions between the fish's internal tissues and the surrounding water. A novel combined experimental-numerical approach shows that free-swimming larval fish generate complex swimming motions using simple actuation patterns.

In this paper, a numerical and experimental study of the shock absorption properties of bike helmets is presented. Laboratory compression and tensile tests were carried out on samples of expanded polystyrene (EPS) and polycarbonate (PC), respectively constituting the internal shock absorption layer and the external hard shell of composite helmets. The measured responses of the two materials were then exploited to calibrate the relevant elasto-plastic constitutive models, adopted in full-scale finite element analyses of a helmet subject to standardized impacts. The simulations allowed assessing the time evolution of the acceleration measured inside the headform (according e.g., to EN 1078) and the failure mechanisms of the helmet, if any, as induced by the localization of plastic deformations.

Numerical assessment of the pressure drop across an aortic coarctation using CFD is a promising approach to replace invasive catheter-based measurements. The aim of this study was to investigate and quantify the uncertainty of numerical calculation of the pressure drop introduced during two essential steps of medical image processing: segmentation of the patient-specific geometry and measurement of patient-specific flow rates from 4D-flow-MRI. Based on the baseline segmentation, geometries with different stenosis diameters were generated for a sample of ten patients. The pressure drop generated by these geometries was calculated for different volume flow rates using computational fluid dynamics. Based on these simulations, a second order polynomial fit was calculated. Based on these polynomial fits an uncertainty of pressure drop calculation was quantified. The calculated pressure drop values varied strongly between the patients. In four patients, pressure drops above and below the clinical threshold of 20 mmHg were found. The median standard deviation of the pressure drop was 2.3 mmHg. The sensitivity of the pressure drop toward changes in the volume flow rate or the stenosis geometry varied between patients. The uncertainty of numerical pressure drop calculation introduced by uncertainties during image segmentation and measurement of volume flow rates was comparable to the uncertainty of pressure drop measurements using invasive catheterization. However, in some patients this uncertainty would have led to different treatment decision. Therefore, patient-specific uncertainty assessment might help to better understand the reliability of a numerically calculated biomarker as the pressure drop across an aortic coarctation.

In this study, laryngeal flow fields are investigated and compared in normal larynx and models of larynx with unilateral vocal fold paralysis (UVFP). In paralytic models, three fixed initial glottal gaps are considered to understand the positive or probable negative impacts of surgical operation on unilaterally paralytic larynx, by which the paralyzed vocal fold is brought closer to the mid-plane. Various features of the flow fields have been discussed in detail including glottal gap width, glottal flow rate, glottal exit pressure pattern and glottal jet evolution. The numerical solution of fluid-structure interaction is carried out using ANSYS, and the results confirm some of the favorable effects of surgery on the patient's larynx. It is also shown that by tightening the glottal gap, some of the problems caused by the presence of a motionless vocal fold, such as leakage through glottal gap in the closure phase resulting in breathy voice can be moderated, although some of the symptoms of this disorder remain relatively unchanged.