Ideal kink instability of a magnetic loop equilibrium

The force-free coronal loop model by Titov & D\'emoulin (1999} is found to be unstable with respect to the ideal kink mode, which suggests this instability as a mechanism for the initiation of flares. The long-wavelength ($m=1$) mode grows for average twists \Phi\ga3.5\pi (at a loop aspect ratio of $\approx$ 5). The threshold of instability increases with increasing major loop radius, primarily because the aspect ratio then also increases. Numerically obtained equilibria at subcritical twist are very close to the approximate analytical equilibrium; they do not show indications of sigmoidal shape. The growth of kink perturbations is eventually slowed down by the surrounding potential field, which varies only slowly with radius in the model. With this field a global eruption is not obtained in the ideal MHD limit. Kink perturbations with a rising loop apex lead to the formation of a vertical current sheet below the apex, which does not occur in the cylindrical approximation.Comment: Astron. Astrophys. Lett., accepte

Financing the new research-based enterprise in New England

Thesis (M.S.)--Massachusetts Institute of Technology, School of Industrial Management, 1963.Includes bibliographical references (leaves 144-147).by Gordon Bruce Baty.M.S

Numerical simulation of two-dimensional single- and multiple-material flow fields

Over the last several years, Sandia National Laboratories has had an interest in developing capabilities to predict the flow fields around vehicles entering or exiting the water at a wide range of speeds. Such prediction schemes have numerous engineering applications in the design of weapon systems. For example, such a scheme could be used to predict the forces and moments experienced by an air-launched anti-submarine weapon on water-entry. Furthermore, a water-exit prediction capability could be used to model the complicated surface closure jet resulting from a missile being shot out of the water. The CCICE (Cell-Centered Implicit Continuous-fluid Eulerian) code developed at Los Alamos National Laboratory (LANL) was chosen to provide the fluid dynamics solver for high speed water-entry and water-exit problems. This implicit time-marching, two-dimensional, conservative, finite-volume code solves the multi-material, compressible, inviscid fluid dynamics equations. The incompressible version of the CCICE code, CCMAC (cell-Centered Marker and Cell), was chosen for low speed water- entry and water-exit problems in order to reduce the computational expense. These codes were chosen to take advantage of certain advances in numerical methods for computational fluid dynamics (CFD) that have taken place at LANL. Notable among these advances is the ability to perform implicit, multi-material, compressible flow simulations, with a fully cell-centered data structure. This means that a single set of control volumes are used, on which a discrete form of the conservation laws is satisfied. This is in control to the more classical staggered mesh methods, in which separate control volumes are defined for mass and momentum. 12 refs

Extragalactic jets with helical magnetic fields: relativistic MHD simulations

Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind the Mach disk, for various jet Lorentz factors and magnetic field helicities. We investigate the influence of helical magnetic fields on jet beam propagation in overdense external medium. We use the AMRVAC code, employing a novel hybrid block-based AMR strategy, to compute ideal plasma dynamics in special relativity. The helicity of the beam magnetic field is effectively transported down the beam, with compression zones in between diagonal internal cross-shocks showing stronger toroidal field regions. In comparison with equivalent low-relativistic jets which get surrounded by cocoons with vortical backflows filled by mainly toroidal field, the high speed jets demonstrate only localized, strong toroidal field zones within the backflow vortical structures. We find evidence for a more poloidal, straight field layer, compressed between jet beam and backflows. This layer decreases the destabilizing influence of the backflow on the jet beam. In all cases, the jet beam contains rich cross-shock patterns, across which part of the kinetic energy gets transferred. For the high speed reference jet considered here, significant jet deceleration only occurs beyond distances exceeding ${\cal O}(100 R_j)$, as the axial flow can reaccelerate downstream to the internal cross-shocks. This reacceleration is magnetically aided, due to field compression across the internal shocks which pinch the flow.Comment: 16 pages, Astronomy and Astrophysics accepted for publicatio

Evidence for a singularity in ideal magnetohydrodynamics: implications for fast reconnection

Numerical evidence for a finite-time singularity in ideal 3D magnetohydrodynamics (MHD) is presented. The simulations start from two interlocking magnetic flux rings with no initial velocity. The magnetic curvature force causes the flux rings to shrink until they come into contact. This produces a current sheet between them. In the ideal compressible calculations, the evidence for a singularity in a finite time $t_c$ is that the peak current density behaves like $|J|_\infty \sim 1/(t_c-t)$ for a range of sound speeds (or plasma betas). For the incompressible calculations consistency with the compressible calculations is noted and evidence is presented that there is convergence to a self-similar state. In the resistive reconnection calculations the magnetic helicity is nearly conserved and energy is dissipated.Comment: 4 pages, 4 figure

A Parametric Study of Erupting Flux Rope Rotation. Modeling the "Cartwheel CME" on 9 April 2008

The rotation of erupting filaments in the solar corona is addressed through a parametric simulation study of unstable, rotating flux ropes in bipolar force-free initial equilibrium. The Lorentz force due to the external shear field component and the relaxation of tension in the twisted field are the major contributors to the rotation in this model, while reconnection with the ambient field is of minor importance. Both major mechanisms writhe the flux rope axis, converting part of the initial twist helicity, and produce rotation profiles which, to a large part, are very similar in a range of shear-twist combinations. A difference lies in the tendency of twist-driven rotation to saturate at lower heights than shear-driven rotation. For parameters characteristic of the source regions of erupting filaments and coronal mass ejections, the shear field is found to be the dominant origin of rotations in the corona and to be required if the rotation reaches angles of order 90 degrees and higher; it dominates even if the twist exceeds the threshold of the helical kink instability. The contributions by shear and twist to the total rotation can be disentangled in the analysis of observations if the rotation and rise profiles are simultaneously compared with model calculations. The resulting twist estimate allows one to judge whether the helical kink instability occurred. This is demonstrated for the erupting prominence in the "Cartwheel CME" on 9 April 2008, which has shown a rotation of \approx 115 degrees up to a height of 1.5 R_sun above the photosphere. Out of a range of initial equilibria which include strongly kink-unstable (twist Phi=5pi), weakly kink-unstable (Phi=3.5pi), and kink-stable (Phi=2.5pi) configurations, only the evolution of the weakly kink-unstable flux rope matches the observations in their entirety.Comment: Solar Physics, submitte

Patient-reported outcomes of periacetabular osteotomy from the prospective ANCHOR cohort study

BACKGROUND: Current literature describing the periacetabular osteotomy (PAO) is mostly limited to retrospective case series. Larger, prospective cohort studies are needed to provide better clinical evidence regarding this procedure. The goals of the current study were to (1) report minimum 2-year patient-reported outcomes (pain, hip function, activity, overall health, and quality of life), (2) investigate preoperative clinical and disease characteristics as predictors of clinical outcomes, and (3) report the rate of early failures and reoperations in patients undergoing contemporary PAO surgery. METHODS: A large, prospective, multicenter cohort of PAO procedures was established, and outcomes at a minimum of 2 years were analyzed. A total of 391 hips were included for analysis (79% of the patients were female, and the average patient age was 25.4 years). Patient-reported outcomes, conversion to total hip replacement, reoperations, and major complications were documented. Variables with a p value of ≤0.10 in the univariate linear regressions were included in the multivariate linear regression. The backward stepwise selection method was used to determine the final risk factors of clinical outcomes. RESULTS: Clinical outcome analysis demonstrated major clinically important improvements in pain, function, quality of life, overall health, and activity level. Increasing age and a body mass index status of overweight or obese were predictive of improved results for certain outcome metrics. Male sex and mild acetabular dysplasia were predictive of lesser improvements in certain outcome measures. Three (0.8%) of the hips underwent early conversion to total hip arthroplasty, 12 (3%) required reoperation, and 26 (7%) experienced a major complication. CONCLUSIONS: This large, prospective cohort study demonstrated the clinical success of contemporary PAO surgery for the treatment of symptomatic acetabular dysplasia. Patient and disease characteristics demonstrated predictive value that should be considered in surgical decision-making. LEVEL OF EVIDENCE: Therapeutic Level IV. See Instructions for Authors for a complete description of levels of evidence

Deterministically Driven Avalanche Models of Solar Flares

We develop and discuss the properties of a new class of lattice-based avalanche models of solar flares. These models are readily amenable to a relatively unambiguous physical interpretation in terms of slow twisting of a coronal loop. They share similarities with other avalanche models, such as the classical stick--slip self-organized critical model of earthquakes, in that they are driven globally by a fully deterministic energy loading process. The model design leads to a systematic deficit of small scale avalanches. In some portions of model space, mid-size and large avalanching behavior is scale-free, being characterized by event size distributions that have the form of power-laws with index values, which, in some parameter regimes, compare favorably to those inferred from solar EUV and X-ray flare data. For models using conservative or near-conservative redistribution rules, a population of large, quasiperiodic avalanches can also appear. Although without direct counterparts in the observational global statistics of flare energy release, this latter behavior may be relevant to recurrent flaring in individual coronal loops. This class of models could provide a basis for the prediction of large solar flares.Comment: 24 pages, 11 figures, 2 tables, accepted for publication in Solar Physic

Src Dependent Pancreatic Acinar Injury Can Be Initiated Independent of an Increase in Cytosolic Calcium

Several deleterious intra-acinar phenomena are simultaneously triggered on initiating acute pancreatitis. These culminate in acinar injury or inflammatory mediator generation in vitro and parenchymal damage in vivo. Supraphysiologic caerulein is one such initiator which simultaneously activates numerous signaling pathways including non-receptor tyrosine kinases such as of the Src family. It also causes a sustained increase in cytosolic calcium- a player thought to be crucial in regulating deleterious phenomena. We have shown Src to be involved in caerulein induced actin remodeling, and caerulein induced changes in the Golgi and post-Golgi trafficking to be involved in trypsinogen activation, which initiates acinar cell injury. However, it remains unclear whether an increase in cytosolic calcium is necessary to initiate acinar injury or if injury can be initiated at basal cytosolic calcium levels by an alternate pathway. To study the interplay between tyrosine kinase signaling and calcium, we treated mouse pancreatic acinar cells with the tyrosine phosphatase inhibitor pervanadate. We studied the effect of the clinically used Src inhibitor Dasatinib (BMS-354825) on pervanadate or caerulein induced changes in Src activation, trypsinogen activation, cell injury, upstream cytosolic calcium, actin and Golgi morphology. Pervanadate, like supraphysiologic caerulein, induced Src activation, redistribution of the F-actin from its normal location in the sub-apical area to the basolateral areas, and caused antegrade fragmentation of the Golgi. These changes, like those induced by supraphysiologic caerulein, were associated with trypsinogen activation and acinar injury, all of which were prevented by Dasatinib. Interestingly, however, pervanadate did not cause an increase in cytosolic calcium, and the caerulein induced increase in cytosolic calcium was not affected by Dasatinib. These findings suggest that intra-acinar deleterious phenomena may be initiated independent of an increase in cytosolic calcium. Other players resulting in acinar injury along with the Src family of tyrosine kinases remain to be explored. © 2013 Mishra et al

Stability of gene contributions and identification of outliers in multivariate analysis of microarray data

BACKGROUND: Multivariate ordination methods are powerful tools for the exploration of complex data structures present in microarray data. These methods have several advantages compared to common gene-by-gene approaches. However, due to their exploratory nature, multivariate ordination methods do not allow direct statistical testing of the stability of genes. RESULTS: In this study, we developed a computationally efficient algorithm for: i) the assessment of the significance of gene contributions and ii) the identification of sample outliers in multivariate analysis of microarray data. The approach is based on the use of resampling methods including bootstrapping and jackknifing. A statistical package of R functions was developed. This package includes tools for both inferring the statistical significance of gene contributions and identifying outliers among samples. CONCLUSION: The methodology was successfully applied to three published data sets with varying levels of signal intensities. Its relevance was compared with alternative methods. Overall, it proved to be particularly effective for the evaluation of the stability of microarray data