Penetrating and blunt trauma to the neck: clinical presentation, assessment and emergency management.

Penetrating and blunt trauma to the neck: clinical presentation, assessment ana emergency management. In Belgium, and even in Western Europe, penetrating and blunt injury to the neck is relatively uncommon in both the civilian and military populations. Pre-hospital and emergency assessment and management will therefore always prove challenging, as individual exposure to this specific type of injury remains low. Historically, the neck has been divided into three anatomical zones with specific landmarks to guide the diagnostic and therapeutic approach to penetrating neck injuries. Most penetrating injuries need to be explored surgically, although with the advent of multi-detector computed tomographic angiography (MDCTA), which yields high diagnostic sensitivity, this inflexible approach has recently changed to a more targeted management, based on clinical, radiographic and, if deemed necessary, endoscopic findings. However, some authors have addressed their concern about this novel, 'no-zone' approach, since the risk of missing less apparent aerodigestive tract injuries may increase. It is recommended, therefore, that all patients with penetrating neck injuries be closely observed, irrespective of the initial findings. The incidence of blunt neck injury is much lower, and this makes risk assessment and management even more difficult in comparison with penetrating injuries. Again, MDCTA is most often the first diagnostic tool if a blunt neck injury is suspected, due to its good sensitivity for blunt cerebrovascular injuries (BCVI) as well as for aerodigestive tract injuries. Specific patterns of injury and unexpected neurological and neuro-radiological findings in trauma patients should always warrant further investigation. Despite ongoing debate, systemic anticoagulation is recommended for most BCVI, sometimes combined with endovascular treatment. Aerodigestive tract injuries may present dramatically, but are often more subtle, making the diagnosis more difficult than other types of neck injuries. Treatment may be conservative if damage is minimal, but surgery is warranted in all other cases.

Low temperature deposition of 2D WS2 layers from WF6 and H2S precursors: impact of reducing agents.

We demonstrate the impact of reducing agents for Chemical Vapor Deposition (CVD) and Atomic Layer Deposition (ALD) of WS2 from WF6 and H2S precursors. Nanocrystalline WS2 layers with a two-dimensional structure can be obtained at low deposition temperatures (300-450 °C) without using a template or anneal.

Influence of valve size, orientation and downstream geometry of an aortic BMHV on leaflet motion and clinically used valve performance parameters.

The aim of this study was to reconcile some of our own previous work and the work of others to generate a physiologically realistic numerical simulation environment that allows to virtually assess the performance of BMHVs. The model incorporates: (i) a left ventricular deformable model to generate a physiological inflow to the aortic valve; (ii) a patient-specific aortic geometry (root, arch and descending aorta); (iii) physiological pressure and flow boundary conditions. We particularly studied the influence of downstream geometry, valve size and orientation on leaflet kinematics and functional indices used in clinical routine. Compared to the straight tube geometry, the patient-specific aorta leads to a significant asynchronous movement of the valve, especially during the closing of the valve. The anterior leaflet starts to close first, impacts the casing at the closed position and remains in this position. At the same time, the posterior leaflet impacts the pivoting mechanisms at the fully open position. At the end of systole, this leaflet subsequently accelerates to the closed position, impacting the casing with an angular velocity of approximately -477 rad/s. The valve size greatly influences the transvalvular pressure gradient (TPG), but does not change the overall leaflet kinematics. This is in contrast to changes in valve orientation, where changing valve orientation induces large differences in leaflet kinematics, but the TPG remains approximately the same.

Functional respiratory imaging as a tool to assess upper airway patency in children with obstructive sleep apnea.

OBJECTIVE: We aim to investigate if anatomical and functional properties of the upper airway using computerized 3D models derived from computed tomography (CT) scans better predict obstructive sleep apnea (OSA) severity than standard clinical markers. METHODS: Consecutive children with suspected OSA underwent polysomnography, clinical assessment of upper airway patency, and a CT scan while awake. A three-dimensional (3D) reconstruction of the pharyngeal airway was built from these images, and computational fluid dynamics modeling of low inspiratory flow was performed using open-source software. RESULTS: Thirty-three children were included (23 boys; mean age, was 6.0±3.2y). OSA was diagnosed in 23 patients. Children with OSA had a significantly lower volume of the overlap region between tonsils and the adenoids (median volume, 1408 mm compared to 2173 mm; p=0.04), a lower mean cross-sectional area at this location (median volume, 69.3mm(2) compared to 114.3mm2; p=0.04), and a lower minimal cross-sectional area (median volume, 17.9 mm2 compared to 25.9 mm2; p=0.05). Various significant correlations were found between several imaging parameters and the severity of OSA, most pronounced for upper airway conductance (r=-0.46) (p<0.01) for correlation between upper airway conductance and the apnea-hypopnea index. No differences or significant correlations were observed with clinical parameters of upper airway patency. Preliminary data after treatment showed that none of the patients with residual OSA had their smallest cross-sectional area located in segment 3, and this frequency was significantly lower than in their peers whose sleep study normalized (64%; p=0.05). CONCLUSION: Functional imaging parameters are highly correlated with OSA severity and are a more powerful correlate than clinical scores of upper airway patency. Preliminary data also showed that we could identify differences in the upper airway of those subjects who did not benefit from a local upper airway treatment.

Haemodynamic impact of stent-vessel (mal)apposition following carotid artery stenting: mind the gaps!

Carotid artery stenting (CAS) has emerged as a minimally invasive alternative to endarterectomy but its use in clinical treatment is limited due to the post-stenting complications. Haemodynamic actors, related to blood flow in the stented vessel, have been suggested to play a role in the endothelium response to stenting, including adverse reactions such as in-stent restenosis and late thrombosis. Accessing the flow-related shear forces acting on the endothelium in vivo requires space and time resolutions which are currently not achievable with non-invasive clinical imaging techniques but can be obtained from image-based computational analysis. In this study, we present a framework for accurate determination of the wall shear stress (WSS) in a mildly stenosed carotid artery after the implantation of a stent, resembling the commercially available Acculink (Abbott Laboratories, Abbott Park, Illinois, USA). Starting from angiographic CT images of the vessel lumen and a micro-CT scan of the stent, a finite element analysis is carried out in order to deploy the stent in the vessel, reproducing CAS in silico. Then, based on the post-stenting anatomy, the vessel is perfused using a set of boundary conditions: total pressure is applied at the inlet, and impedances that are assumed to be insensitive to the presence of the stent are imposed at the outlets. Evaluation of the CAS outcome from a geometrical and haemodynamic perspective shows the presence of atheroprone regions (low time-average WSS, high relative residence time) colocalised with stent malapposition and stent strut interconnections. Stent struts remain unapposed in the ostium of the external carotid artery disturbing the flow and generating abnormal shear forces, which could trigger thromboembolic events.

A computational exploration of helical arterio-venous graft designs.

Although arterio-venous grafts (AVGs) are the second best option as long-term vascular access for hemodialysis, they suffer from complications caused by intimal hyperplasia, mainly located in vessel regions of low and oscillating wall shear stress. However, certain flow patterns in the bulk may reduce these unfavorable hemodynamic conditions. We therefore studied, with computational fluid dynamics (CFD), the impact of a helical AVG design on the occurrence of (un)favorable hemodynamic conditions at the venous anastomosis. Six CFD-models of an AVG in closed-loop configuration were constructed: one conventional straight graft, and five helical designed grafts with a pitch of 105 mm down to 35 mm. At the venous anastomosis, disturbed shear was assessed by quantifying the area with unfavorable conditions, and by analyzing averaged values in a case-specific patch. The bulk hemodynamics were assessed by analyzing the kinetic helicity in and the pressure drop over the graft. The most helical design scores best, being instrumental to suppress disturbed shear in the venous segment. There is, however, no trivial relationship between the number of helix turns of the graft and disturbed shear in the venous segment, when a realistic closed-loop AVG model is investigated. Bulk flow investigation showed a marked increase of helicity intensity in, and a moderate pressure drop over the AVG by introducing a lower pitch. At the venous anastomosis, unfavorable hemodynamic conditions can be reduced by introducing a helical design. However, due to the complex flow conditions, the optimal helical design for an AVG cannot be derived without studying case by case.

Computational models for wall shear stress estimation in scaffolds: a comparative study of two complete geometries.

Fluid mechanical stimuli are known to upregulate cell differentiation and matrix formation. Since wall shear stress plays an important role various studies tried to estimate the scaffold fluid dynamic environment. However, because of the geometrical complexity, nearly all studies created their CFD model based on a submodel of the entire scaffold assuming that the model covers heterogeneity sufficiently. However to the authors' knowledge no study exist providing guidelines in this matter. In a previous study we demonstrated that submodels are influenced by the boundary conditions, inevitable when flow channels are chopped off. For the current study we therefore developed µCT based models of two complete scaffold geometries (one titanium and one hydroxyapatite). Imposing a 0.04 ml/min flow rate resulted in a surface area averaged wall shear stress of 1.41 mPa for titanium and 1.09 mPa for hydroxyapatite. In order to get insight in required model size we subdivided the domain in regions of different size. From our results we propose a model size between 6 and 10 times the average pore size. The wall shears stress should be calculated on a region at least one pore size away from the boundaries. These guidelines could be of use for computationally more costly simulations where it is not possible to simulate the complete scaffold domain.

Experimental validation of a pulse wave propagation model for predicting hemodynamics after vascular access surgery.

Hemodialysis patients require a vascular access that is, preferably, surgically created by connecting an artery and vein in the arm, i.e. an arteriovenous fistula (AVF). The site for AVF creation is chosen by the surgeon based on preoperative diagnostics, but AVFs are still compromised by flow-associated complications. Previously, it was shown that a computational 1D-model is able to describe pressure and flow after AVF surgery. However, predicted flows differed from measurements in 4/10 patients. Differences can be attributed to inaccuracies in Doppler measurements and input data, to neglecting physiological mechanisms or to an incomplete physical description of the pulse wave propagation after AVF surgery. The physical description can be checked by validating against an experimental setup consisting of silicone tubes mimicking the aorta and arm vasculature both before and after AVF surgery, which is the aim of the current study. In such an analysis, the output uncertainty resulting from measurement uncertainty in model input should be quantified. The computational model was fed by geometrical and mechanical properties collected from the setup. Pressure and flow waveforms were simulated and compared with experimental waveforms. The precision of the simulations was determined by performing a Monte Carlo study. It was concluded that the computational model was able to simulate mean pressures and flows accurately, whereas simulated waveforms were less attenuated than experimental ones, likely resulting from neglecting viscoelasticity. Furthermore, it was found that in the analysis output uncertainties, resulting from input uncertainties, cannot be neglected and should thus be considered.

How can cells sense the elasticity of a substrate? An analysis using a cell tensegrity model.

A eukaryotic cell attaches and spreads on substrates, whether it is the extracellular matrix naturally produced by the cell itself, or artificial materials, such as tissue-engineered scaffolds. Attachment and spreading require the cell to apply forces in the nN range to the substrate via adhesion sites, and these forces are balanced by the elastic response of the substrate. This mechanical interaction is one determinant of cell morphology and, ultimately, cell phenotype. In this paper we use a finite element model of a cell, with a tensegrity structure to model the cytoskeleton of actin filaments and microtubules, to explore the way cells sense the stiffness of the substrate and thereby adapt to it. To support the computational results, an analytical 1D model is developed for comparison. We find that (i) the tensegrity hypothesis of the cytoskeleton is sufficient to explain the matrix-elasticity sensing, (ii) cell sensitivity is not constant but has a bell-shaped distribution over the physiological matrix-elasticity range, and (iii) the position of the sensitivity peak over the matrix-elasticity range depends on the cytoskeletal structure and in particular on the F-actin organisation. Our model suggests that F-actin reorganisation observed in mesenchymal stem cells (MSCs) in response to change of matrix elasticity is a structural-remodelling process that shifts the sensitivity peak towards the new value of matrix elasticity. This finding discloses a potential regulatory role of scaffold stiffness for cell differentiation.

Flip angle optimization for dynamic contrast-enhanced MRI-studies with spoiled gradient echo pulse sequences.

Spoiled gradient echo pulse (SPGRE) sequences are commonly used in dynamic contrast-enhanced MRI (DCE-MRI) studies to measure the contrast agent concentration in a tissue of interest over time. However, due to improper tuning of the SPGRE parameters, concentration uncertainty can be very high, even at high signal-to-noise ratio in the MR measurement. In this work, an optimization procedure is proposed for selecting the optimal value of the SPGRE-flip angle FA(opt), given the expected concentration range. The optimization condition ensures that every concentration in the assumed range has the lowest possible uncertainty. By decoupling the R(1)- and R*(2)-effects caused by the presence of the contrast agent, a contour plot has been generated from which FA(opt) can be read off for any study design. Investigation of ten recent DCE-MRI studies showed that improper flip angle selection unnecessarily increases the concentration uncertainty, up to 742% and 72% on average for the typical physiological concentration ranges of 0-2 mM in tumour tissue and 0-10 mM in blood, respectively. Simulations show that the reduced noise levels on the concentration curves, observed at the optimal flip angle, effectively increase the precision of the kinetic parameters estimates (up to 82% for K(trans), 82% for ν(e) and 92% for ν(p) in the case of an individually measured arterial input function (AIF), up to 53% for K(trans), 59% for ν(e) and 67% for ν(p) in the case of a standard AIF). In vivo experiments confirm the potential of flip angle optimization to increase the reproducibility of the kinetic parameter estimates.

Patient-specific computational haemodynamics: generation of structured and conformal hexahedral meshes from triangulated surfaces of vascular bifurcations.

Measuring the blood flow is still limited by current imaging technologies and is generally overcome using computational fluid dynamics (CFD) which, because of the complex geometry of blood vessels, has widely relied on tetrahedral meshes. Hexahedral meshes offer more accurate results with lower-density meshes and faster computation as compared to tetrahedral meshes, but their use is limited by the far more complex mesh generation. We present a robust methodology for conformal and structured hexahedral mesh generation - applicable to complex arterial geometries as bifurcating vessels - starting from triangulated surfaces. Cutting planes are used to slice the lumen surface and to construct longitudinal Bezier splines. Afterwards, an isoparametric transformation is used to map a parametrically defined quadrilateral surface mesh into the vessel volume, resulting in stacks of sections which can then be used for sweeping. Being robust and open source based, this methodology may improve the current standard in patient-specific mesh generation and enhance the reliability of CFD to patient-specific haemodynamics.

In vitro, in vivo and numerical assessment of the working principle of the truCCOMS continuous cardiac output catheter system.

The truCCOMS cardiac output monitor system provides a continuous and instantaneous measurement of cardiac output, derived from the amount of energy required for heating a filament to maintain a fixed 2 degrees C blood temperature difference between two thermistors located distally on a pulmonary artery catheter. Clinical studies, however, reported relatively poor accuracy of the cardiac output estimation, possibly due to linearly assumed power-cardiac output relationship used for calibration of the catheters. We experimentally studied the shape of the truCCOMS calibration relationship (i) in a hydraulic bench model of the right heart and (ii) in vivo intact animal model. The results showed a nonlinear relationship between the power input into the heating element and the cardiac output; which could satisfactorily be described with an exponential relationship. Comparison of the performance of the same catheters in vitro and in vivo showed that the in vitro determined calibration relationship should not be used for in vivo measurements. Finally, we also simulated the working principle of the catheter using a simplified numerical model of the blood flow and heat transfer around the catheter. The computed results also suggested a pronounced nonlinear relationship between power and cardiac output in pulsatile conditions. We conclude that the observed over- and underestimation of high- and low flows, respectively, by the current truCCOMS system is likely to arise from its linear calibration relationship. An appropriate calibration scheme accounting for the intrinsic nonlinear power-cardiac output relationship and the difference between in vitro and in vivo conditions should improve the clinical performance of the system.

An assessment of ductus venosus tapering and wave transmission from the fetal heart.

Pressure and flow pulsations in the fetal heart propagate through the precordial vein and the ductus venosus (DV) but are normally not transmitted into the umbilical vein (UV). Pulsations in the umbilical vein do occur, however, in early pregnancy and in pathological conditions. Such transmission into the umbilical vein is not well understood. In particular, the effect of the impedance changes in the DV due to its tapered geometry is not known. This paper presents a mathematical model that we developed to study the transmission of pulsations, originating in the fetal heart, through the DV to the umbilical vein. In our model, the tapered geometry of the DV was found to be of minor importance and the only effective reflection site in the DV appears to be at the DV inlet. Differences between the DV inlet and outlet flow were also found to be minor for medium to large umbilical vein-DV diameter ratios. Finally, the results of a previously proposed lumped model were found to agree well with the present model of the DV-umbilical vein bifurcation.

Determining carotid artery pressure from scaled diameter waveforms: comparison and validation of calibration techniques in 2026 subjects.

Calibrated diameter distension waveforms could provide an alternative for local arterial pressure assessment more widely applicable than applanation tonometry. We compared linearly and exponentially calibrated carotid diameter waveforms to tonometry readings. Local carotid pressures measured by tonometry and diameter waveforms measured by ultrasound were obtained in 2026 subjects participating in the Asklepios study protocol. Diameter waveforms were calibrated using a linear and an exponential calibration scheme and compared to measured tonometry waveforms by examining the mean root-mean-squared error (RMSE), carotid systolic blood pressure (SBPcar) and augmentation index (AIx) of calibrated and measured pressures. Mean RMSE was 5.2(3.3) mmHg (mean(stdev)) for linear and 4.6(3.6) mmHg for exponential calibration. Linear calibration yielded an underestimation of SBPcar by 6.4(4.1) mmHg which was strongly correlated to values of brachial pulse pressure (PPbra) (R = 0.4, P < 0.05). Exponential calibration underestimated true SBPcar by 1.9(3.9) mmHg, independent of PPbra. AIx was overestimated by linear calibration by 1.9(10.1)%, the difference significantly increasing with increasing AIx (R = 0.25, P < 0.001) and by exponential calibration by 5.4(10.6)%, independently of the value of AIx. Properly calibrated diameter waveforms offer a viable alternative for local pressure estimation at the carotid artery. Compared to linear calibration, exponential calibration significantly improves the pressure estimation.

Three- and four-element Windkessel models: assessment of their fitting performance in a large cohort of healthy middle-aged individuals.

Lumped-parameter models are used to estimate the global arterial properties by fitting the model to measured (aortic) pressure and flow. Different model configurations coexist, and it is still an open question as to which model optimally reflects the arterial tree and leads to correct estimates of arterial properties. An assessment was made of the performance of (a) the three-element Windkessel model (WK3) consisting of vascular resistance R, total arterial compliance C, and characteristic impedance Zc; (b) a four-element model with an inertance element L placed in parallel with Zc (WK4-p); and (c) a four-element model with L placed in series with Zc (WK4-s). Models were fitted to data measured non-invasively in 2404 healthy subjects, aged between 35 and 55 years. It was found that model performance segregated into two groups. In a group containing 20 per cent of the dataset (characterized by low blood pressure and wave reflection) the WK4-p model outperformed the other models, with model behaviour as envisioned by its promoters. In these cases, the WK3 and WK4-s models led to increased overestimation of total arterial compliance and underestimation of characteristic impedance. However, in about 80 per cent of the cases, the WK4-p model showed a behaviour that was very similar to that of the WK3 and WK4-s models. Here, the WK4-s model yielded the best quality of fit, although model parameters reached physically impossible values for L in about 12 per cent of all cases. The debate about which lumped-parameter model is the better approximation of the arterial tree is therefore still not fully resolved.

Experimental and numerical modelling of the ventriculosinus shunt (El-Shafei shunt).

This study assesses malresorptive hydrocephalus treatment by ventriculosinus shunting with the shunt in the antegrade or retrograde position. First, an experimental model of the cerebral ventricles, the arachnoid villi, the cortical veins, and the superior sagittal sinus was built. For this purpose, the compliance of a human cortical vein was measured and then modelled by means of Penrose tubes. The dimensions of the superior sagittal sinus were determined in vivo by measurements on magnetic resonance imaging scans of 21 patients. Second, a numerical model of the cortical veins and the superior sagittal sinus was built. The numerical results were validated with the results from the experimental model. The experimental and numerical pressure difference between the intracranial pressure and the static sinus pressure was small (0-20 Pa) and corresponded to the theoretically expected values. No overdrainage was found in either the antegrade or the retrograde position of the shunt. Blood reflow was only found while mimicking lumbar puncture or changes in position with the experimental model (lowering the intracranial pressure or increasing the sinus pressure rapidly). Optimal results can be obtained with the shunt positioned in the most downstream half of the superior sagittal sinus. The experimental and numerical results confirm the potential of ventriculosinus shunting as therapy for malresorptive hydrocephalus patients. The ventriculosinus shunt thus proves to be a promising technique.

Numerical study of the uniformity of balloon-expandable stent deployment.

Stents are small tubelike structures, implanted in coronary and peripheral arteries to reopen narrowed vessel sections. This endovascular intervention remains suboptimal, as the success rate is limited by restenosis. This renarrowing of a stented vessel is related to the arterial injury caused by stent-artery and balloon-artery interactions, and a local subsequent inflammatory process. Therefore, efforts to optimize the stent deployment remain very meaningful. Several authors have studied with finite element modeling the mechanical behavior of balloon-expandable stents, but none of the proposed models incorporates the folding pattern of the balloon. We developed a numerical model in which the CYPHER stent is combined with a realistic trifolded balloon. In this paper, the impact of several parameters such as balloon length, folding pattern, and relative position of the stent with respect to the balloon catheter on the free stent expansion has been investigated. Quantitative validation of the modeling strategy shows excellent agreement with data provided by the manufacturer and, therefore, the model serves as a solid basis for further investigations. The parametric analyses showed that both the balloon length and the folding pattern have a considerable influence on the uniformity and symmetry of the transient stent expansion. Consequently, this approach can be used to select the most appropriate balloon length and folding pattern for a particular stent design in order to optimize the stent deployment. Furthermore, it was demonstrated that small positioning inaccuracies may change the expansion behavior of a stent. Therefore, the placement of the stent on the balloon catheter should be accurately carried out, again in order to decrease the endothelial damage.

Impact of hemodialysis duration on the removal of uremic retention solutes.

Several studies have stressed the importance of dialysis time in the removal of uremic retention solutes. To further investigate this, nine stable chronic hemodialysis patients were dialyzed for 4, 6, or 8 h processing the same total blood and dialysate volume by the Genius system and high-flux FX80 dialyzers. Inlet blood and outlet dialysate were analyzed for urea, creatinine, phosphorus, and beta2-microglobulin at various times. Total solute removal, dialyzer extraction ratios, and total cleared volumes were significantly larger during prolonged dialysis for urea, creatinine, phosphorus, and beta2-microglobulin. Reduction ratios increased progressively, except for phosphate and beta2-microglobulin, where the ratios remained constant after 2 h. In contrast, no significant difference was found for the reduction ratios of all solutes and Kt/V(urea) between the three different sessions. With longer dialyses, solutes are efficiently removed from the deeper compartments of the patient's body. Our study shows that care must be taken when using Kt/Vurea or reduction ratios as the only parameters to quantify dialysis adequacy.

PIV validation of blood-heart valve leaflet interaction modelling.

The aim of this study was to validate the 2D computational fluid dynamics (CFD) results of a moving heart valve based on a fluid-structure interaction (FSI) algorithm with experimental measurements. Firstly, a pulsatile laminar flow through a monoleaflet valve model with a stiff leaflet was visualized by means of Particle Image Velocimetry (PIV). The inflow data sets were applied to a CFD simulation including blood-leaflet interaction. The measurement section with a fixed leaflet was enclosed into a standard mock loop in series with a Harvard Apparatus Pulsatile Blood Pump, a compliance chamber and a reservoir. Standard 2D PIV measurements were made at a frequency of 60 bpm. Average velocity magnitude results of 36 phase-locked measurements were evaluated at every 10 degrees of the pump cycle. For the CFD flow simulation, a commercially available package from Fluent Inc. was used in combination with inhouse developed FSI code based on the Arbitrary Lagrangian-Eulerian (ALE) method. Then the CFD code was applied to the leaflet to quantify the shear stress on it. Generally, the CFD results are in agreement with the PIV evaluated data in major flow regions, thereby validating the FSI simulation of a monoleaflet valve with a flexible leaflet. The applicability of the new CFD code for quantifying the shear stress on a flexible leaflet is thus demonstrated.