Finite Element Based Nonlinear Normalization of Human Lumbar Intervertebral Disc Stiffness to Account for Its Morphology

Disc degeneration, usually associated with low back pain and changes of intervertebral stiffness, represents a major health issue. As the intervertebral disc (IVD) morphology influences its stiffness, the link between mechanical properties and degenerative grade is partially lost without an efficient normalization of the stiffness with respect to the morphology. Moreover, although the behavior of soft tissues is highly nonlinear, only linear normalization protocols have been defined so far for the disc stiffness. Thus, the aim of this work is to propose a nonlinear normalization based on finite elements (FE) simulations and evaluate its impact on the stiffness of human anatomical specimens of lumbar IVD. First, a parameter study involving simulations of biomechanical tests (compression, flexion/extension, bilateral torsion and bending) on 20 FE models of IVDs with various dimensions was carried out to evaluate the effect of the disc's geometry on its compliance and establish stiffness/morphology relations necessary to the nonlinear normalization. The computed stiffness was then normalized by height (H), cross-sectional area (CSA), polar moment of inertia (J) or moments of inertia (Ixx, Iyy) to quantify the effect of both linear and nonlinear normalizations. In the second part of the study, T1-weighted MRI images were acquired to determine H, CSA, J, Ixx and Iyy of 14 human lumbar IVDs. Based on the measured morphology and pre-established relation with stiffness, linear and nonlinear normalization routines were then applied to the compliance of the specimens for each quasi-static biomechanical test. The variability of the stiffness prior to and after normalization was assessed via coefficient of variation (CV). The FE study confirmed that larger and thinner IVDs were stiffer while the normalization strongly attenuated the effect of the disc geometry on its stiffness. Yet, notwithstanding the results of the FE study, the experimental stiffness showed consistently higher CV after normalization. Assuming that geometry and material properties affect the mechanical response, they can also compensate for one another. Therefore, the larger CV after normalization can be interpreted as a strong variability of the material properties, previously hidden by the geometry's own influence. In conclusion, a new normalization protocol for the intervertebral disc stiffness in compression, flexion, extension, bilateral torsion and bending was proposed, with the possible use of MRI and FE to acquire the discs' anatomy and determine the nonlinear relations between stiffness and morphology. Such protocol may be useful to relate the disc's mechanical properties to its degree of degeneration.

Sudden death in a case of sickle cell anemia: post-mortem computed tomography and autopsy correlation from a radiologist's perspective

Sickle cell anemia (SCA) is a hemolytic disease characterized by the production of abnormal hemoglobin chains and distorted red blood cell morphology or sickling. "Sickle cell crisis" includes vaso-occlusive crisis, a plastic crisis, sequestration crisis, haemolytic crisis and often culminating in serious complications, organ damage and even sudden death. Post-mortem computed tomography (PMCT) findings of sickle cell disease have never been reported in literature. This case of sudden death from acute hemolytic crisis in SCA where post-mortem computed tomography (PMCT) and autopsy findings complemented each other, both revealing findings invisible to the other and both crucial to the case.

Serial [18F]-fluoromisonidazole PET during radiochemotherapy for locally advanced head and neck cancer and its correlation with outcome.

PURPOSE The aim was to assess changes of tumour hypoxia during primary radiochemotherapy (RCT) for head and neck cancer (HNC) and to evaluate their relationship with treatment outcome. MATERIAL AND METHODS Hypoxia was assessed by FMISO-PET in weeks 0, 2 and 5 of RCT. The tumour volume (TV) was determined using FDG-PET/MRI/CT co-registered images. The level of hypoxia was quantified on FMISO-PET as TBRmax (SUVmaxTV/SUVmean background). The hypoxic subvolume (HSV) was defined as TV that showed FMISO uptake ⩾1.4 times blood pool activity. RESULTS Sixteen consecutive patients (T3-4, N+, M0) were included (mean follow-up 31, median 44months). Mean TBRmax decreased significantly (p<0.05) from 1.94 to 1.57 (week 2) and 1.27 (week 5). Mean HSV in week 2 and week 5 (HSV2=5.8ml, HSV3=0.3ml) were significantly (p<0.05) smaller than at baseline (HSV1=15.8ml). Kaplan-Meier plots of local recurrence free survival stratified at the median TBRmax showed superior local control for less hypoxic tumours, the difference being significant at baseline and after 2weeks (p=0.031, p=0.016). CONCLUSIONS FMISO-PET documented that in most HNC reoxygenation starts early during RCT and is correlated with better outcome.

Oxygen isotopes of lake marl at Gerzensee and Leysin (Switzerland), covering the Younger Dryas and two minor oscillations, and their correlation to the GRIP ice core

The ratio of oxygen isotopes is a temperature proxy both in precipitation and in the calcite of lacustrine sediments. The very similar oxygen-isotope records from Greenland ice cores and European lake sediments during the Last Glacial Termination suggest that the drastic climatic changes occurred quasi-simultaneously on an extra-regional, probably hemispheric scale. In order to study temporal relations of the different parameters recorded in lake sediments, for example biotic response times to rapid climatic changes, a precise chronology is required. In unlaminated lake sediments there is not yet available a method to provide a high-resolution chronology, especially for periods with radiocarbon plateaux. Alternatively, an indirect time scale can be constructed by linking the lake stratigraphy with other well-dated climate records. New oxygen-isotope records from Gerzensee and Leysin, with an estimated sampling resolution of between 15 and 40 years, match the Greenlandic isotope record in many details. Under the assumption that the main variations in temperature and thus in oxygen isotopes occurred about simultaneously in Greenland and Switzerland, we have assigned a time scale to the lake sediments of Gerzensee and Leysin by wiggle-matching their stable-isotope records with those of Greenland ice cores, which are among the best dated climatic archives. We estimate a precision of 20 to 100 years during the Last Glacial Termination.

Late Cretaceous-Cenozoic magnetostratigraphic and biostratigraphic correlation of DSDP Hole 73-519 in the South Atlantic

DSDP Leg 73 sediment cores allow direct calibrations of magnetostratigraphy and biostratigraphy for much of the latest Cretaceous to Cenozoic in the mid-latitude South Atlantic Ocean. A complete record of the Cenozoic was not obtained, however, because strong dissolution, poor core recovery and intense core disturbance have masked the biostratigraphy or magnetostratigraphy over some intervals of all recovered sections. DSDP Leg 73 results show the following correlations: Early/middle Miocene in Chron 16 Oligocene/Miocene within Subchron C6CN Eocene/Oligocene within Subchron C13R Middle/late Eocene top of Chron C17 Early/late Paleocene top of Subchron C27N Cretaceous/Tertiary within Subchron C29R

Column Generation Algorithms for Nonlinear Optimization II: Numerical Investigations

García et al. present a class of column generation (CG) algorithms for nonlinear programs. Its main motivation from a theoretical viewpoint is that under some circumstances, finite convergence can be achieved, in much the same way as for the classic simplicial decomposition method; the main practical motivation is that within the class there are certain nonlinear column generation problems that can accelerate the convergence of a solution approach which generates a sequence of feasible points. This algorithm can, for example, accelerate simplicial decomposition schemes by making the subproblems nonlinear. This paper complements the theoretical study on the asymptotic and finite convergence of these methods given in [1] with an experimental study focused on their computational efficiency. Three types of numerical experiments are conducted. The first group of test problems has been designed to study the parameters involved in these methods. The second group has been designed to investigate the role and the computation of the prolongation of the generated columns to the relative boundary. The last one has been designed to carry out a more complete investigation of the difference in computational efficiency between linear and nonlinear column generation approaches. In order to carry out this investigation, we consider two types of test problems: the first one is the nonlinear, capacitated single-commodity network flow problem of which several large-scale instances with varied degrees of nonlinearity and total capacity are constructed and investigated, and the second one is a combined traffic assignment model

Anticorrelation between Ion Acceleration and Nonlinear Coherent Structures from Laser-Underdense Plasma Interaction

In laser-plasma experiments, we observed that ion acceleration from the Coulomb explosion of the plasma channel bored by the laser, is prevented when multiple plasma instabilities such as filamentation and hosing, and nonlinear coherent structures (vortices/post-solitons) appear in the wake of an ultrashort laser pulse. The tailoring of the longitudinal plasma density ramp allows us to control the onset of these insabilities. We deduced that the laser pulse is depleted into these structures in our conditions, when a plasma at about 10% of the critical density exhibits a gradient on the order of 250 {\mu}m (gaussian fit), thus hindering the acceleration. A promising experimental setup with a long pulse is demonstrated enabling the excitation of an isolated coherent structure for polarimetric measurements and, in further perspectives, parametric studies of ion plasma acceleration efficiency.

Estimation of a probability in inverse binomial sampling under normalized linear-linear and inverse-linear loss

Sequential estimation of the success probability p in inverse binomial sampling is considered in this paper. For any estimator pˆ , its quality is measured by the risk associated with normalized loss functions of linear-linear or inverse-linear form. These functions are possibly asymmetric, with arbitrary slope parameters a and b for pˆ

and pˆ>p , respectively. Interest in these functions is motivated by their significance and potential uses, which are briefly discussed. Estimators are given for which the risk has an asymptotic value as p→0, and which guarantee that, for any p∈(0,1), the risk is lower than its asymptotic value. This allows selecting the required number of successes, r, to meet a prescribed quality irrespective of the unknown p. In addition, the proposed estimators are shown to be approximately minimax when a/b does not deviate too much from 1, and asymptotically minimax as r→∞ when a=b.

Comparing linear and non linear wind flow models

Assessing wind conditions on complex terrain has become a hard task as terrain complexity increases. That is why there is a need to extrapolate in a reliable manner some wind parameters that determine wind farms viability such as annual average wind speed at all hub heights as well as turbulence intensities. The development of these tasks began in the early 90´s with the widely used linear model WAsP and WAsP Engineering especially designed for simple terrain with remarkable results on them but not so good on complex orographies. Simultaneously non-linearized Navier Stokes solvers have been rapidly developed in the last decade through CFD (Computational Fluid Dynamics) codes allowing simulating atmospheric boundary layer flows over steep complex terrain more accurately reducing uncertainties. This paper describes the features of these models by validating them through meteorological masts installed in a highly complex terrain. The study compares the results of the mentioned models in terms of wind speed and turbulence intensity.