SCIENTIFIC COMPUTING AND IMAGING INSTITUTE
at the University of Utah

It is well known that the violence of slow cook-off explosions can greatly exceed the comparatively mild case burst events typically observed for rapid heating. However, there have been few studies that examine the reaction violence as a function of applied heat flux that explore the dependence on heating geometry and device size. Here we report progress on a study using the Uintah Computation Framework, a high-performance computer model capable of modeling deflagration, material damage, deflagration to detonation transition and detonation for PBX9501 and similar explosives. Our results suggests the existence of a sharp threshold for increased reaction violence with decreasing heat flux. The critical heat flux was seen to increase with increasing device size and decrease with the heating of multiple surfaces, suggesting that the temperature gradient in the heated energetic material plays an important role the violence of reactions.

Keywords: DDT, cook-off, deflagration, detonation, violence of reaction, c-safe

We propose a novel difference metric, called the graph diffusion distance (GDD), for quantifying the difference between two weighted graphs with the same number of vertices. Our approach is based on measuring the average similarity of heat diffusion on each graph. We compute the graph Laplacian exponential kernel matrices, corresponding to repeatedly solving the heat diffusion problem with initial conditions localized to single vertices. The GDD is then given by the Frobenius norm of the difference of the kernels, at the diffusion time yielding the maximum difference. We study properties of the proposed distance on both synthetic examples, and on real-data graphs representing human anatomical brain connectivity.

In this paper we present a novel Bayesian approach for fractional segmentation of white matter tracts and simultaneous estimation of a multi-tensor diffusion model. Our model consists of several white matter tracts, each with a corresponding weight and tensor compartment in each voxel. By incorporating a prior that assumes the tensor fields inside each tract are spatially correlated, we are able to reliably estimate multiple tensor compartments in fiber crossing regions, even with low angular diffusion-weighted imaging (DWI). Our model distinguishes the diffusion compartment associated with each tract, which reduces the effects of partial voluming and achieves more reliable statistics of diffusion measurements.We test our method on synthetic data with known ground truth and show that we can recover the correct volume fractions and tensor compartments. We also demonstrate that the proposed method results in improved segmentation and diffusion measurement statistics on real data in the presence of crossing tracts and partial voluming.

An objective measurement technique to quantify 3D femoral head shape was developed and applied to normal subjects and patients with cam-type femoroacetabular impingement (FAI). 3D reconstructions were made from high-resolution CT images of 15 cam and 15 control femurs. Femoral heads were fit to ideal geometries consisting of rotational conchoids and spheres. Geometric similarity between native femoral heads and ideal shapes was quantified. The maximum distance native femoral heads protruded above ideal shapes and the protrusion area were measured. Conchoids provided a significantly better fit to native femoral head geometry than spheres for both groups. Cam-type FAI femurs had significantly greater maximum deviations (4.99 ± 0.39 mm and 4.08 ± 0.37 mm) than controls (2.41 ± 0.31 mm and 1.75 ± 0.30 mm) when fit to spheres or conchoids, respectively. The area of native femoral heads protruding above ideal shapes was significantly larger in controls when a lower threshold of 0.1 mm (for spheres) and 0.01 mm (for conchoids) was used to define a protrusion. The 3D measurement technique described herein could supplement measurements of radiographs in the diagnosis of cam-type FAI. Deviations up to 2.5 mm from ideal shapes can be expected in normal femurs while deviations of 4–5 mm are characteristic of cam-type FAI.

Statistical shape modeling (SSM) was used to quantify 3D variation and morphologic differences between femurs with and without cam femoroacetabular impingement (FAI). 3D surfaces were generated from CT scans of femurs from 41 controls and 30 cam FAI patients. SSM correspondence particles were optimally positioned on each surface using a gradient descent energy function. Mean shapes for groups were defined. Morphological differences between group mean shapes and between the control mean and individual patients were calculated. Principal component analysis described anatomical variation. Among all femurs, the first six modes (or principal components) captured significant variations, which comprised 84% of cumulative variation. The first two modes, which described trochanteric height and femoral neck width, were significantly different between groups. The mean cam femur shape protruded above the control mean by a maximum of 3.3 mm with sustained protrusions of 2.5–3.0 mm along the anterolateral head-neck junction/distal anterior neck. SSM described variations in femoral morphology that corresponded well with areas prone to damage. Shape variation described by the first two modes may facilitate objective characterization of cam FAI deformities; variation beyond may be inherent population variance. SSM could characterize disease severity and guide surgical resection of bone.

Advances in computational mechanics, constitutive modeling, and techniques for subject-specific modeling have opened the door to patient-specific simulation of the relationships between joint mechanics and osteoarthritis (OA), as well as patient-specific preoperative planning. This article reviews the application of computational biomechanics to the simulation of joint contact mechanics as relevant to the study of OA. This review begins with background regarding OA and the mechanical causes of OA in the context of simulations of joint mechanics. The broad range of technical considerations in creating validated subject-specific whole joint models is discussed. The types of computational models available for the study of joint mechanics are reviewed. The types of constitutive models that are available for articular cartilage are reviewed, with special attention to choosing an appropriate constitutive model for the application at hand. Issues related to model generation are discussed, including acquisition of model geometry from volumetric image data and specific considerations for acquisition of computed tomography and magnetic resonance imaging data. Approaches to model validation are reviewed. The areas of parametric analysis, factorial design, and probabilistic analysis are reviewed in the context of simulations of joint contact mechanics. Following the review of technical considerations, the article details insights that have been obtained from computational models of joint mechanics for normal joints; patient populations; the study of specific aspects of joint mechanics relevant to OA, such as congruency and instability; and preoperative planning. Finally, future directions for research and application are summarized.

Elastin is a structural protein that provides resilience to biological tissues. We examined the contributions of elastin to the quasi-static tensile response of porcine medial collateral ligament through targeted disruption of the elastin network with pancreatic elastase. Elastase concentration and treatment time were varied to determine a dose response. Whereas elastin content decreased with increasing elastase concentration and treatment time, the change in peak stress after cyclic loading reached a plateau above 1 U/ml elastase and 6 h treatment. For specimens treated with 2 U/ml elastase for 6 h, elastin content decreased approximately 35%. Mean peak tissue strain after cyclic loading (4.8%, p ≥ 0.300), modulus (275 MPa, p ≥ 0.114) and hysteresis (20%, p ≥ 0.553) were unaffected by elastase digestion, but stress decreased significantly after treatment (up to 2 MPa, p ≤ 0.049). Elastin degradation had no effect on failure properties, but tissue lengthened under the same pre-stress. Stiffness in the linear region was unaffected by elastase digestion, suggesting that enzyme treatment did not disrupt collagen. These results demonstrate that elastin primarily functions in the toe region of the stress–strain curve, yet contributes load support in the linear region. The increase in length after elastase digestion suggests that elastin may pre-stress and stabilize collagen crimp in ligaments

Background
A contributory factor to hip osteoarthritis (OA) is abnormal cartilage mechanics. Acetabular retroversion, a version deformity of the acetabulum, has been postulated to cause OA via decreased posterior contact area and increased posterior contact stress. Although cartilage mechanics cannot be measured directly in vivo to evaluate the causes of OA, they can be predicted using finite element (FE) modeling.

Objective
The objective of this study was to compare cartilage contact mechanics between hips with normal and retroverted acetabula using subject-specific FE modeling.

Methods
Twenty subjects were recruited and imaged: 10 with normal acetabula and 10 with retroverted acetabula. FE models were constructed using a validated protocol. Walking, stair ascent, stair descent and rising from a chair were simulated. Acetabular cartilage contact stress and contact area were compared between groups.

Results
Retroverted acetabula had superomedial cartilage contact patterns, while normal acetabula had widely distributed cartilage contact patterns. In the posterolateral acetabulum, average contact stress and contact area during walking and stair descent were 2.6–7.6 times larger in normal than retroverted acetabula (P ≤ 0.017). Conversely, in the superomedial acetabulum, peak contact stress during walking was 1.2–1.6 times larger in retroverted than normal acetabula (P ≤ 0.044). Further differences varied by region and activity.

Conclusions
This study demonstrated superomedial contact patterns in retroverted acetabula vs widely distributed contact patterns in normal acetabula. Smaller posterolateral contact stress in retroverted acetabula than in normal acetabula suggests that increased posterior contact stress alone may not be the link between retroversion and OA.

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.

The AGENT simulation system is used for detailed three-dimensional modeling of neutron transport and corresponding properties of nuclear reactors of any design. Numerical solution to the neutron transport equation in the AGENT system is based on the Method of Characteristics (MOCs) and the theory of R-functions. The latter of which is used for accurately describing current and future heterogeneous lattices of reactor core configurations. The AGENT code has been extensively verified to assure a high degree of accuracy for predicting neutron three-dimensional point-wise flux spatial distributions, power peaking factors, reaction rates, and eigenvalues. In this paper, a new AGENT code feature, a computational steering, is presented. This new feature provides a novel way for using deterministic codes for fast evaluation of reactor core parameters, at no loss to accuracy. The computational steering framework as developed at the Technische Universität München is smoothly integrated into the AGENT solver. This framework allows for an arbitrary interruption of AGENT simulation, allowing the solver to restart with updated parameters. One possible use of this is to accelerate the convergence of the final values resulting in significantly reduced simulation times. Using this computational steering in the AGENT system, coarse MOC resolution parameters can initially be selected and later update them – while the simulation is actively running – into fine resolution parameters. The utility of the steering framework is demonstrated using the geometry of a research reactor at the University of Utah: this new approach provides a savings in CPU time on the order of 50%.

Keywords: Numerical steering, AGENT code, Deterministic neutron transport codes, Method of Characteristics, R-functions, Numerical visualizations

Background: Magnetic resonance imaging (MRI) can visualize locations of both the ablation scar on the left atrium (LA) after atrial fibrillation (AF) ablation and epicardial fat pads (FPs) containing ganglionated plexi (GP).

Methods: We investigated 60 patients who underwent pulmonary vein antrum (PVA) isolation along with LA posterior wall and septal debulking for AF. FPs around the LA surface in well-known GP areas (which were considered as the substitution of GP areas around the LA) were segmented from the dark-blood MRI. Then the FP and the ablation scar image visualized by late gadolinium enhancement (LGE)-MRI on the LA were merged together. Overlapping areas of FP and the ablation scar image were considered as the ablated FP areas containing GP. Patients underwent 24-hour Holter monitoring after ablation for the analysis of heart rate variability.

Results: Ablated FP area was significantly wider in patients without AF recurrence than those in patients with recurrence (5.6 ± 3.1 cm² vs 4.2 ± 2.7 cm² ,P = 0.03). The mean values of both percentage of differences greater than 50 ms in the RR intervals (pRR > 50) and standard deviation of RR intervals over the entire analyzed period (SDNN), which were obtained from 24-hour Holter monitoring 1-day post-AF ablation, were significantly lower in patients without recurrence than those in patients with recurrence (5.8 ± 6.0% vs 14.0 ± 10.1%; P = 0.0005, 78.7 ± 32.4 ms vs 109.2 ± 43.5 ms; P = 0.005). There was a significant negative correlation between SDNN and the percentage of ablated FP area (Y =- 1.3168X + 118.96, R² = 0.1576, P = 0.003).

Conclusion: Extensively ablating LA covering GP areas along with PVA isolation enhanced the denervation of autonomic nerve system and seemed to improve procedural outcome in patients with AF.

Keywords: ganglionated plexi, fat pad, atrial fibrillation, catheter ablation, LGE-MRI

We present a novel integrated visualization system that enables interactive visual analysis of ensemble simulations used in ocean forecasting, i.e, simulations of sea surface elevation. Our system enables the interactive planning of both the placement and operation of off-shore structures. We illustrate this using a real-world simulation of the Gulf of Mexico. Off-shore structures, such as those used for oil exploration, are vulnerable to hazards caused by strong loop currents. The oil and gas industry therefore relies on accurate ocean forecasting systems for planning their operations. Nowadays, these forecasts are based on multiple spatio-temporal simulations resulting in multidimensional, multivariate and multivalued data, so-called ensemble data. Changes in sea surface elevation are a good indicator for the movement of loop current eddies, and our visualization approach enables their interactive exploration and analysis. We enable analysis of the spatial domain, for planning the placement of structures, as well as detailed exploration of the temporal evolution at any chosen position, for the prediction of critical ocean states that require the shutdown of rig operations.

Keywords: Uncertainty, Ensemble Simulation, Risk Estimate

In this paper, we propose a new and accurate technique for uncertainty analysis and uncertainty visualization based on fiber orientation distribution function (ODF) glyphs, associated with high angular resolution diffusion imaging (HARDI). Our visualization applies volume rendering techniques to an ensemble of 3D ODF glyphs, which we call SIP functions of diffusion shapes, to capture their variability due to underlying uncertainty. This rendering elucidates the complex heteroscedastic structural variation in these shapes. Furthermore, we quantify the extent of this variation by measuring the fraction of the volume of these shapes, which is consistent across all noise levels, the certain volume ratio. Our uncertainty analysis and visualization framework is then applied to synthetic data, as well as to HARDI human-brain data, to study the impact of various image acquisition parameters and background noise levels on the diffusion shapes.

We introduce a novel method for utilizing user input to sparsely label membranes in electron microscopy images. Using gridlines as guides, the user marks where the guides cross the membrane to generate a sparsely labeled image. We use a best path algorithm to connect each of the sparse membrane labels. The resulting segmentation has a significantly better Rand error than automatic methods while requiring as little as 2\% of the image to be labeled.

In connectomics, neuroscientists seek to identify the synaptic connections between neurons. Segmentation of cell membranes using supervised learning algorithms on electron microscopy images of brain tissue is often done to assist in this effort. Here we present a partial differential equation with a novel growth term to improve the results of a supervised learning algorithm. We also introduce a new method for representing the resulting image that allows for a more dynamic thresholding to further improve the result. Using these two processes we are able to close small to medium sized gaps in the cell membrane detection and improve the Rand error by as much as 9\% over the initial supervised segmentation.